Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 1,513,121 C→A 7.8% F135L (TTC→TTA)  mraY → Phospho‑N‑acetylmuramoyl‑pentapeptide‑transferas e

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1431,513,1210CA 7.8% 196.9 / 4.2 67F135L (TTC→TTA) mraYPhospho‑N‑acetylmuramoyl‑pentapeptide‑transferas e
Reads supporting (aligned to +/- strand):  ref base C (24/36);  new base A (4/2);  total (29/38)
Fisher's exact test for biased strand distribution p-value = 3.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

CGACCTGAGCAACCGCTATGTATGGGTTGTGCTGATTGTCACCCTGGCGTTCGGTGCCATTGGCTGGGTCGATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGGTGTTCCTGTACAAGACGGCGCCAACCAGCGTCGAGACCACGCTGATCCTGCCGTTCATCAAGGATGTCACCATTCCGTTGGGCGTCGGCTTCGTCGTACTGAC  >  ACB143/1512989‑1513255
                                                                                                                                    |                                                                                                                                      
cGACCTGAGCAACCGCTATGTATGGGTTGTGCTGATTGTCACCCTGGCGTTCGGTGCCATTGGCTGGGTCGATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTgg                                                                                                                                     >  1:229947/1‑136 (MQ=255)
 gACCTGAGCAACCGCTATGTATGGGTTGTGCTGATTGTCACCCTGGCGTTCGGTGCCATTGGCTGGGTCGATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGc                                                                                                                                    <  4:65121/136‑1 (MQ=255)
 gACCTGAGCAACCGCTATGTATGGGTTGTGCTGATTGTCACCCTGGCGTTCGGTGCCATTGGCTGGGTCGATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGc                                                                                                                                    <  2:122207/136‑1 (MQ=255)
     tGAGCAACCGCTATGTATGGGTTGTGCTGATTGTCACCCTGGCGTTCGGTGCCATTGGCTGGGTCGATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATc                                                                                                                                >  4:108652/1‑136 (MQ=255)
            ccGCTATGTATGGGTTGTGCTGATTGTCACCCTGGCGTTCGGTGCCATTGGCTGGGTCGATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTc                                                                                                                         <  8:119125/136‑1 (MQ=255)
                    tatgGGTTGTGCTGATTGTCACCCTGGCGTTCGGTGCCATTGGCTGGGTCGATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTggc                                                                                                                 >  6:25304/1‑136 (MQ=255)
                      tgGGTTGTGCTGATTGTCACCCTGGCGTTCGGTGCCATTGGCTGGGTCGATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTggcgg                                                                                                               <  1:263770/136‑1 (MQ=255)
                         ggtgtgCTGATTGTCACCCTGGCGTTCGGTGCCATTGGCTGGGTCGATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTggcggcgg                                                                                                            <  5:198668/134‑1 (MQ=255)
                         gTTGTGCTGATTGTCACCCTGGCGTTCGGTGCCATTGGCTGGGTCGATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTggcggcgg                                                                                                            <  6:99132/136‑1 (MQ=255)
                         gTTGTGCTGATTGTCACCCTGGCGTTCGGTGCCATTGGCTGGGTCGATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTggcggcgg                                                                                                            <  6:215783/136‑1 (MQ=255)
                            gtgCTGATTGTCACCCTGGCGTTCGGTGCCATTGGCTGGGTCGATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGcggcggcgg                                                                                                         >  7:61070/1‑136 (MQ=255)
                            gtgCTGATTGTCACCCTGGCGTTCGGTGCCATTGGCTGGGTCGATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTAGGCCTGGcggcggcgg                                                                                                         >  4:59901/1‑136 (MQ=255)
                             tgCTGATTGTCACCCTGGCGTTCGGTGCCATTGGCTGGGTCGATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGcggcggcgg                                                                                                         <  8:175761/135‑1 (MQ=255)
                             tgCTGATTGTCACCCTGGCGTTCGGTGCCATTGGCTGGGTCGATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGgt                                                                                                        >  3:117036/1‑136 (MQ=255)
                              gCTGATTGTCACCCTGGCGTTCGGTGCCATTGGCTGGGTCGATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGgtg                                                                                                       >  3:120179/1‑136 (MQ=255)
                               cTGATTGTCAGCCTGTCGTTCGTTGCCATTGGCTGGGTCGATGACTACCTCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGgtg                                                                                                       <  8:196470/135‑1 (MQ=255)
                               cTGATTGTCACCCTGGCGTTCGGTGCCATTGGCTGGGTCGATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGgtgt                                                                                                      >  6:173753/1‑136 (MQ=255)
                               cTGATTGTCACCCTGGCGTTCGGTGCCATTGGCTGGGTCGATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGgtgt                                                                                                      >  3:163125/1‑136 (MQ=255)
                               cTGATTGTCACCCTGGCGTTCGGTGCCATTGGCTGGGTCGATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGgtgt                                                                                                      >  6:96912/1‑136 (MQ=255)
                               cTGATTGTCACCCTGGCGTTCGGTGCCATTGGCTGGGTCGATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGgtgt                                                                                                      >  7:133458/1‑136 (MQ=255)
                               cTGATTGTCACCCTGGCGTTCGGTGCCATTGGCTGGGTCGATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGgtgt                                                                                                      >  7:57680/1‑136 (MQ=255)
                               cTGATTGTCACCCTGGCGTTCGGTGCCATTGGCTGGGTCGATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGgtg                                                                                                       <  7:97017/135‑1 (MQ=255)
                               cTGATTGTCACCCTGGCGTTCGGTGCCATTGGCTGGGTCGATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGgtg                                                                                                       <  5:29263/135‑1 (MQ=255)
                               cTGATTGTCACCCTGGCGTTCGGTGCCATTGGCTGGGTCGATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGgtg                                                                                                       <  5:143188/135‑1 (MQ=255)
                               cTGATTGTCACCCTGGCGTTCGGTGCCATTGGCTGGGTCGATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGgtg                                                                                                       <  5:122014/135‑1 (MQ=255)
                               cTGATTGTCACCCTGGCGTTCGGTGCCATTGGCTGGGTCGATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTATGGCAATCGGTGTTCGGCCTGGCGGCGGCGgggt                                                                                                      >  6:187050/1‑136 (MQ=255)
                               cTGATTGTCACCCTGGCGTTCGGTGCCATTGGCTGGGTCGATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACGTCAGGCAATCGGTGTTCGGCCTGGCGGCGGCGgtgt                                                                                                      >  5:105853/1‑136 (MQ=255)
                                tGATTGTCACCCTGGCGTTCGGTGCCATTGGCTGGGTCGATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGGTGtt                                                                                                     >  7:196470/1‑136 (MQ=255)
                                 gATTGTCACCCTGGCGTTCGGTGCCATTGGCTGGGTCGATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGTAGTACTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGGTGTTc                                                                                                    >  4:27938/1‑136 (MQ=255)
                                 gATTGTCACCCTGGCGTTCGGTGCCATTGGCTGGGTCGATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGGTGTTc                                                                                                    >  5:198533/1‑136 (MQ=255)
                                 gATTGTCACCCTGGCGTTCGGTGCCATTGGCTGGGTCGATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTATGGCAATCGGTGTTCGGCCTGGCGGCGGCGGTGTTc                                                                                                    >  1:85491/1‑136 (MQ=255)
                                 gATTGTCACCCTGGCGTTCGGTGCCATTGGCTGGGTCGATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTATGGCAATCGGTGTTCGGCCTGGCGGCGGCGGTGTTc                                                                                                    >  5:284596/1‑136 (MQ=255)
                                   ttGTCACCCTGGCGTTCGGTGCCATTGGCTGGGTCGATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGGTGTTCCt                                                                                                  <  8:133458/136‑1 (MQ=255)
                                        aCCCTGGCGTTCGGTGCCATTGGCTGGGGCGATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGGTGTTCCTGTACa                                                                                             <  6:177889/136‑1 (MQ=255)
                                          ccTGGCGTTCGGTGCCATTGGCTGGGTCGATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGGTGTTCCTGTAc                                                                                              >  6:212685/1‑133 (MQ=255)
                                              gCGTTCGGTGCCATTGGCTGGGTCGATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGGTGTTCCTGTACAAGACg                                                                                        <  4:117036/135‑1 (MQ=255)
                                               cGTTCGGTGCCATTGGCTGGGTCGATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGGTGTTCCTGTACAAGACGgc                                                                                      >  3:112484/1‑136 (MQ=255)
                                               cGTTCGGTGCCATTGGCTGGGTCGATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCGGGAAGTACTTGTGGCAATCGGTGTTCGGCCTGGCGGCGGCGGGGTTGCTTGACAAGAAGgc                                                                                      >  6:38289/1‑136 (MQ=255)
                                                    ggggCCATTGGCTGGGTCGATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGGTGTTCCTGTACAAGACGGCGCCaa                                                                                 <  7:258976/133‑1 (MQ=255)
                                                    ggTGCCATTGGCTGGGTCGATGACTACCGCAAGGTTATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGGTGTTCCTGTACAAGACGGCGCCa                                                                                  <  6:175692/135‑1 (MQ=255)
                                                      tGCCATTGGCTGGGTCGATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGGTGTTCCTGTACAAGACGGCGCCAAcc                                                                               >  2:149182/1‑136 (MQ=255)
                                                              gCTGGGTCGATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGGTGTTCCTGTACAAGACGGCGCCAACCAGCGTCga                                                                       <  8:61070/136‑1 (MQ=255)
                                                              gCTGGGTCGATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGGTGTTCCTGTACAAGACGGCGCCAACCAGCGTCga                                                                       <  2:27889/136‑1 (MQ=255)
                                                                   gTCGATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGGTGTTCCTGTACAAGACGGCGCCAACCAGCGTCGAGACCa                                                                  >  7:152150/1‑136 (MQ=255)
                                                                     agATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTCGCCTTCCGAGCCGCAGGAAGTACTTCTGGAAATCGGTGATCGGCCTGGCGGCGGAGGTGTTCCTGTACAACAAGGCGACACACATCGTCACGACCACg                                                                >  6:213628/2‑136 (MQ=255)
                                                                          aCTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGGTCTTACTGTACAAGACGGCGCCAACCAGCGTCGAGACCACGcagct                                                           >  5:198444/1‑132 (MQ=255)
                                                                              ccGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGGTGTTCCTGTACAAGACGGCGCCAACCAGCGTCGAGACCACGCTGATCCTg                                                       <  4:112484/136‑1 (MQ=255)
                                                                                gCAAGGTGATCGAAAAGAACTCGCGTTGCGTTCCGAGCCGCTGGAAGTACTGAGGGCAATCGGTGTTCGGCCTGGCGGCGGGGGTGGTCGTGTTCAAGGCGGGGCCAACCCGCGTCGAGACCCCGCTGATccgg                                                       >  5:58489/1‑132 (MQ=255)
                                                                                gCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGGTGTTCCTGTACAAGACGGCGCCAACCAGCGTCGAGACCACGCTGATCCTGcc                                                     <  6:146434/136‑1 (MQ=255)
                                                                                gCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGGTGTTCCTGTACAAGACGGCGCCAACCAGCGTCGAGACCACGCTGATCCTGcc                                                     <  1:136407/136‑1 (MQ=255)
                                                                                gCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGGTGTTCCTGTACAAGACGGCGCCAACCAGCGTCGAGACCACGCTGATCCTGcc                                                     <  4:248159/136‑1 (MQ=255)
                                                                                    ggTGATCGAATAGAAGTCGGGTGGCCTTCCGAGGCGCGGTTTGTACTTCTGGCAATCGGTGTTCGGCCTTGCGGCGGCGGTGTTCCTGTACAAGACGGCGCCAACCAGCGTCGAGACCACGCTGATCCTGCCGt                                                   <  6:155484/134‑1 (MQ=255)
                                                                                     gTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGGTGTTCCTGTACAAGACGGCGCCAACCAGCGTCGAGACCACGCTGATCCTGCCGTtca                                                <  2:229947/136‑1 (MQ=255)
                                                                                         tCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCTGTGTTCGGCCTGGCGGCGGCGGTGTTCCTGTACAAGACGGCGCCAACCAGCGTCGAGACCACGCTGATCCTGCCGTTCATCaa                                            <  3:164056/136‑1 (MQ=255)
                                                                                         tCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGGTGTTCCTGTACAAGACGGCGCCAACCAGCGTCGAGACCACGCTGATCCTGCCGTtcatc                                              >  7:152352/1‑134 (MQ=255)
                                                                                         tCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGGTGTTCCTGTACAAGACGGCGCCAACCAGCGTCGAGACCACGCTGATCCTGCCGTTCATCaa                                            <  3:189121/136‑1 (MQ=255)
                                                                                            aaaaGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGGTGTTCCTGTACAAGACGGCGCCAACCAGCGTCGAGACCACGCTGATCCTGCCGTTCATCAAGGa                                         <  4:107304/136‑1 (MQ=255)
                                                                                              aaGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGGTGTTCCTGTACAAGACGGCGCCAACCAGCGTCGAGACCACGCTGATCCTGCCGTTCATCAAGGATg                                       <  1:261802/136‑1 (MQ=255)
                                                                                                gAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGGTGTTCCTGTACAAGACGGCGCCAACCAGCGTCGAGACCACGCTGATCCTGCCGTTCATCAAGGATGTc                                     <  1:77255/136‑1 (MQ=255)
                                                                                                        gTGGCCTTCCGAGCCGCTGGAAGTACTTATGGCAATCGGTGTTCGGCCTGGCGGCGGCGGTGTTCCTGTACAAGACGGCGCCAACCAGCGTCGAGACCACGCTGATCCTGCCGTTCATCAAGGATGTCACCATTcc                             <  6:284596/136‑1 (MQ=255)
                                                                                                        gTGGCCTTCCGAGCCGCTGGAAGTACTTATGGCAATCGGTGTTCGGCCTGGCGGCGGCGGTGTTCCTGTACAAGACGGCGCCAACCAGCGTCGAGACCACGCTGATCCTGCCGTTCATCAAGGATGTCACCATTcc                             <  2:85491/136‑1 (MQ=255)
                                                                                                                 cGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGGTGTTCCTGTACAAGACGGCGCCAACCAGCGTCGAGACCACGCTGATCCTGCCGTTCATCAAGGATGTCACCATTCCGTTGGGCGt                    <  1:36470/136‑1 (MQ=255)
                                                                                                                   aGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGGTGTTCCTGTACAAGACGGCGCCAACCAGCGTCGAGACCACGCTGATCCTGCCGTTCATCAAGGATGTCACCATTCCGTTGGGCGTCg                  >  5:153133/1‑136 (MQ=255)
                                                                                                                   aGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGGTGTTCCTGTACAAGACGGCGCCAACCAGCGTCGAGACCACGCTGATCCTGCCGTTCATCAAGGATGTCACCATTCCGTTGGGCGTCg                  <  4:163125/136‑1 (MQ=255)
                                                                                                                      cGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGGTGTTCCTGTACAAGACGGCGCCAACCAGCGTCGAGACCACGCTGATCCTGCCGTTCATCAAGGATGTCACCATTCCGTTGGGCGTCGGCt               <  2:250788/136‑1 (MQ=255)
                                                                                                                      cGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGGTGTTCCTGTACAAGACGGCGCCAACCAGCGTCGAGACCACGCTGATCCTGCCGTTCATCAAGGATGTCACCATTCCGTTGGGCGTCGGCt               <  1:149182/136‑1 (MQ=255)
                                                                                                                                 cTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGGTGTTCCTGTACAAGACGGCGCCAACCAGCGTCGAGACCACGCTGATCCTGCCGTTCATCAAGGATGTCACCATTCCGTTGGGCGTCGGCTTCGTCGTACt     <  6:198533/135‑1 (MQ=255)
                                                                                                                                    cTGGCAATCGGTGTTCGGCCTGGCGGCGGCGGTGTTCCTGTACAAGACGGCGCCAACCAGCGTCGAGACCACGCTGATCCTGCCGTTCATCAAGGATGTCACCATTCCGTTGGGCGTCGGCTTCGTCGTACTGAc  <  4:201810/135‑1 (MQ=255)
                                                                                                                                    cTGGCAATCGGTGTTCGGCCTGGCGGCGGCGGTGTTCCTGTACAAGACGGCGCCAACCAGCGTCGAGACCACGCTGATCCTGCCGTTCATCAAGGATGTCACCATTCCGTTGGGCGTCGGCTTCGTCGTACTGAc  <  3:160709/135‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                      
CGACCTGAGCAACCGCTATGTATGGGTTGTGCTGATTGTCACCCTGGCGTTCGGTGCCATTGGCTGGGTCGATGACTACCGCAAGGTGATCGAAAAGAACTCGCGTGGCCTTCCGAGCCGCTGGAAGTACTTCTGGCAATCGGTGTTCGGCCTGGCGGCGGCGGTGTTCCTGTACAAGACGGCGCCAACCAGCGTCGAGACCACGCTGATCCTGCCGTTCATCAAGGATGTCACCATTCCGTTGGGCGTCGGCTTCGTCGTACTGAC  >  ACB143/1512989‑1513255

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑