Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,781,103 A→C 12.5% intergenic (+86/‑112) lldD → / → trmL L‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,781,1030AC12.5% 270.7 / 4.1 121intergenic (+86/‑112)lldD/trmLL‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (51/51);  new base C (7/10);  total (60/61)
Fisher's exact test for biased strand distribution p-value = 6.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 123 total aligned reads displayed.

TAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCG  >  NC_000913/3780966‑3781250
                                                                                                                                         |                                                                                                                                                   
tAAAGAGTCGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCTATTGCGCCCCCTCACCATAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATTTACTCCCTCACTCTGAAACGACACCGCACTCtt                                                                                                                                           <  5:107467/148‑1 (MQ=255)
 aaaGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCCCCCTAACCCTCTACCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCttt                                                                                                                                          >  1:211274/1‑148 (MQ=255)
 aaaGAGTTGCCTGCGGCACTGGCTCCCAAGGCGAAAGCGAATTCGGCATAGTCGTTTTCGCCCCCTAACCCTAAACCAATCCCTCAGGGCGAAGGGACCCGTCGGCGCAGTAAGGACACCCTCACTCTGTAACGAAACCCCACTAttt                                                                                                                                          >  6:18952/1‑148 (MQ=255)
      gTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACACttttttt                                                                                                                                      >  2:319797/1‑147 (MQ=255)
       ttGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtct                                                                                                                                    <  5:192476/148‑1 (MQ=255)
        tGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCCCGCTGAAACGCCCCCGCCCTCTTTTTTtctc                                                                                                                                   >  4:139757/1‑148 (MQ=255)
        tGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACGCCCTCACTCTCCAACGACCCCCCCCTCGTTTTTtctc                                                                                                                                   >  7:312355/1‑148 (MQ=255)
         gCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtctc                                                                                                                                   <  2:287677/147‑1 (MQ=255)
          ccTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCCCCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGGATGTACTCCCTCACGCTGAAACGACACCGCACTCtttttttatacc                                                                                                                                 >  3:143771/1‑143 (MQ=255)
          ccTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTccc                                                                                                                                 <  6:354019/148‑1 (MQ=255)
              cGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGc                                                                                                                             >  5:373116/1‑148 (MQ=255)
              cGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTAACCCTAACCCACTCCCACAGGGAGAGGGGACCGTTCGGCGATGGATGTACACCCTCACTCTGAAACGACACCGCACTCGTTTTTTCTCCCTCGc                                                                                                                             >  4:126087/1‑148 (MQ=255)
                  aCTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTAGCCCCt                                                                                                                         <  7:180252/148‑1 (MQ=255)
                    tGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCACCCTAACCCTCTCCCGCAGGGAGAGGGGCCCGTTCGGGGCTGTATGTACTCCCTCACTCTGAAACGCCCCCGCCCTCGTTTTTTCTCCCTCGCCCCGcc                                                                                                                       >  7:367961/1‑148 (MQ=255)
                      gCTCCCCAGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCgg                                                                                                                     <  5:282352/148‑1 (MQ=255)
                      gCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCgg                                                                                                                     <  8:329560/148‑1 (MQ=255)
                       cTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCCCTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCggg                                                                                                                    >  1:308287/1‑148 (MQ=255)
                         cccATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCggg                                                                                                                    >  4:205310/1‑146 (MQ=255)
                          ccATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGgag                                                                                                                 <  3:348792/148‑1 (MQ=255)
                          ccATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGgag                                                                                                                 <  1:302267/148‑1 (MQ=255)
                           cATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCACCCTAAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGCCACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGga                                                                                                                  >  2:111812/1‑146 (MQ=255)
                              ggCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTGAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAggg                                                                                                             >  5:75866/1‑148 (MQ=255)
                                 gAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCg                                                                                                          <  4:345511/148‑1 (MQ=255)
                                  aaaGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACGCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCgg                                                                                                         >  1:117682/1‑148 (MQ=255)
                                   aaGGGAATGCGGCATAGTCGTTTGCGCCCCCTCGCCCTAACCCTCTCCCTCAGGGAGAGTGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACCACACAGCACACATATTATCTCCCTCCCCCCTCCGCGGACAGGccccgc                                                                                                        >  5:118236/1‑145 (MQ=255)
                                      ggAATGCGGCATAGTCGTTTGCCCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTg                                                                                                     <  7:63199/148‑1 (MQ=255)
                                           gCGGCATAGTCGTTTGCGCCCCCTAACCCTAACCCTCACCCTCAGGGAGAGGGGACCGTTAGGCGCTGTATGTACTCCCTCCCACTGAAACGCAACCGCCCTATTTTTTTCACCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAggg                                                                                                 >  8:263898/1‑147 (MQ=255)
                                            cGGCAATGTAGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCCGGGAGAGGGGAGCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCTGTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGATGGCCGGGGTGAGGGGa                                                                                               <  6:241557/148‑1 (MQ=255)
                                            aggCATAGTCTTTTGCTCCCCCTCACCCTAACCCTCTCCCTCAGTGATAGGGGACCTTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTATCTCCATCGCCCCTCCGGGGAGAGGGCAGGGGTGAGGGGa                                                                                               <  6:141444/147‑1 (MQ=255)
                                              gCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCCGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaa                                                                                             <  5:117234/148‑1 (MQ=255)
                                              gCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTGGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaa                                                                                             <  5:219313/148‑1 (MQ=255)
                                              gCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGgag                                                                                                                 >  1:72205/1‑128 (MQ=255)
                                              gCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGgag                                                                                                                 <  2:72205/128‑1 (MQ=255)
                                              gCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaa                                                                                             <  2:325236/148‑1 (MQ=255)
                                               cATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaa                                                                                              >  8:317447/1‑146 (MQ=255)
                                                 tAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAgg                                                                                          <  7:198945/148‑1 (MQ=255)
                                                  aGTCGTTTGCGCCCCCTCACCCTAACCCTATCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTACCTCACTCTTAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGTAAAGGc                                                                                         >  3:247002/1‑148 (MQ=255)
                                                  aGTCGTTTGCGCCCCCGCACCCAAACCCACTCCCGCAGGGGGAGGGGGCCGTTCGGGGCGGTTTGTTCTACCACCCACTGAAACGCCCCCGCACTCTTTTTTTCTCCCACGCCCCACCGGGGCGAGGGCCGGGGGGGGGGGAAAATGc                                                                                         >  6:17912/1‑148 (MQ=255)
                                                   gTCGTTTGCGCCCCCTCAACCTAACCCACTCCCACAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCAAACTCTGAAACGACAACGCACTCTTTTTTTATACCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGcc                                                                                        >  8:231234/1‑148 (MQ=255)
                                                      gTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGTGTGAGGGGAAAAGGCCGc                                                                                      <  5:313561/147‑1 (MQ=255)
                                                       tttGCGCCCCCTCACCCAAACCCTCACCCTCAGGGAGAGGGGACCGATCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCACCGGGGAGAGGGCCGGGGTCAGGGGAAAAGGCCGCa                                                                                     >  4:331274/1‑147 (MQ=255)
                                                                 cTCACCCTAACCCTTTCCCTCAGGGAGTGGGGACCGTGGGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGTGTGAGGGGAAAAGGCCGCACTGCTCACCaa                                                                          <  3:325116/148‑1 (MQ=255)
                                                                 cTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAgg                                                                                          >  6:336042/1‑132 (MQ=255)
                                                                 cTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAgg                                                                                          <  5:336042/132‑1 (MQ=255)
                                                                  tCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAt                                                                         >  6:167049/1‑148 (MQ=255)
                                                                   cACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGCGGGCCGGGGTGCGGGGAAAAGGCCGCACTGCGAACCAAtt                                                                        >  6:66798/1‑148 (MQ=255)
                                                                    acccacaCCCTCTCCGTCAGGGTGAGGGGACCGTGCGGCGCTGTATGTACTCCCCCACTCTCAAACCACCCCGCACTATTTTTTTCTCCCTCCCCCCTCGGGGGATTGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAt                                                                         <  6:37226/140‑1 (MQ=255)
                                                                    aCCCTAACCCTCTCCCTCTGGGTGTGGGGACCGTTCGGCGCTGTAGGTACTCCCTCACTCTCAAACGACCCCGCACTCTTTTTTTCCCCCTCGCCCCTCGGGGGTGCGGGCGGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAttt                                                                       <  2:178140/148‑1 (MQ=255)
                                                                     cccTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAtttt                                                                      >  3:178439/1‑148 (MQ=255)
                                                                      ccTAACCCGCGCCCTCAGGGCGAGGGGACCGTGCGGCGCGGTTTGTACTCCCCCCCCGTGAAACCCCCCCGCACGCTTTTTTTCTCCCTCGCCCCGCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTc                                                                     <  2:270998/148‑1 (MQ=255)
                                                                      caaaccccTCGCCCCCAGGGCGCGGGGACCGTGCGGCGCTGTATGAACCCCCCCACTCTGAACCGACCCCGCACTCTTTTTTTCTCCCTCCCCCCGCCGGGGAGAGGGCGGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTc                                                                     <  5:17912/143‑1 (MQ=255)
                                                                       cTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCt                                                                    >  4:306172/1‑148 (MQ=255)
                                                                         aaCCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCCGTCTGAAACGACCCCGCACACTTTTTTTCTCCCTCGCCCCTCCGGGGAGGGGCCCGGGGTGAGGGGAAAAGGCCGCACTGGTCCCAAATTTTCTGc                                                                  >  4:4562/1‑148 (MQ=255)
                                                                           cccTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGgag                                                                                                                 >  6:62401/1‑99 (MQ=255)
                                                                           cccTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGgag                                                                                                                 <  5:62401/99‑1 (MQ=255)
                                                                             ctctCCCTGAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaa                                                                                             <  1:154130/117‑1 (MQ=255)
                                                                             ctctCCCTGAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaa                                                                                             >  2:154130/1‑117 (MQ=255)
                                                                             ctctCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaa                                                                                             <  6:55883/117‑1 (MQ=255)
                                                                             ctctCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaa                                                                                             >  5:55883/1‑117 (MQ=255)
                                                                             ctctCCCTCAGGGAGAGGGGACCGGTCGGCGCTGTATGTTCGCCCTCACTCTGAAACGCCACCGGACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATc                                                              >  4:70152/1‑148 (MQ=255)
                                                                              tctcCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGCCGCCGCCCTCTTTTTTTCTCCCTCGCCCCTCCGGGGCGAGGGCGGGGGGGCGGGGAAAAGGCCGCACTGCGCACCAATTTTCTGCTATcc                                                             >  2:193528/1‑148 (MQ=255)
                                                                              tctcCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATcc                                                             >  1:261041/1‑148 (MQ=255)
                                                                                  ccTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTg                                                           <  2:261041/146‑1 (MQ=255)
                                                                                   cTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACCCCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGccc                                                        >  4:128460/1‑148 (MQ=255)
                                                                                    tCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGcccc                                                       >  8:224243/1‑148 (MQ=255)
                                                                                     cAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGccccc                                                      >  1:63939/1‑148 (MQ=255)
                                                                                      aGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGcccc                                                       <  6:155366/146‑1 (MQ=255)
                                                                                       gggAGAGGCGACGGTTCGCCTCAGTATCTCCTCCCTCACGCTGACACTCCCCCGCCCTATTTTTTTCTCCCTCGCCCCTCCGGGGCGCGGGCCGTGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGc                                                    <  5:305901/148‑1 (MQ=255)
                                                                                         gagaGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCAc                                                  >  4:214634/1‑148 (MQ=255)
                                                                                           gtggggACCGTTGGGCGCTGTATGTACTCCCCCACTCGGAACCGCCCCCGCACTTTTTTTTTCTCCCTCGCCCCTCGGGGGAGCGGGCGGGGGGGCGGGGACAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTa                                                <  1:362123/146‑1 (MQ=255)
                                                                                           gatgggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCTCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTa                                                <  5:18952/145‑1 (MQ=255)
                                                                                           gagGGGACCGTTCGGCGCTGTATGTACGCCCTCACTCTGAAACGACCCCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTa                                                <  2:308287/148‑1 (MQ=255)
                                                                                            tggggACCGTTGGGCGCTGTATTTACTCCCTCACTGTGAAAGGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCGGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTa                                                <  4:168142/146‑1 (MQ=255)
                                                                                            agGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGTGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTaa                                               <  3:331274/148‑1 (MQ=255)
                                                                                                 aCCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggggg                                          >  4:199578/1‑148 (MQ=255)
                                                                                                 aCCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggggg                                          <  7:317447/148‑1 (MQ=255)
                                                                                                  ccGTTCGGCGCTGTATGTACTCCCTCACTCTGAAATGAAACCGCACTCTTTTTTTCTCCCTCGCCCGTCCGGGGAGAGGGCCGGGGTTAGGGGAAAAGACCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGc                                         <  6:114412/148‑1 (MQ=255)
                                                                                                  ccGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGc                                         <  3:128460/148‑1 (MQ=255)
                                                                                                  ccGTGCTGCGCTGTATTTCCGCCCTCACTCTGAACCCCCCCCGCACTTTTTTTTTCTCCCTCGCCCCTCCGGGGAGTGGGCCGGGGTGAGGGGACAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGc                                         <  6:75866/148‑1 (MQ=255)
                                                                                                    gTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAg                                       >  1:307521/1‑148 (MQ=255)
                                                                                                       cGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGCGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTa                                     <  3:306172/147‑1 (MQ=255)
                                                                                                       cGGCGCTGTATATCCTCCCTCCCCCTGAAACGACACCGCACGCTTTTTTTCTCCCTCGCCCCTCCGGGGCGAGGGCCGGGGGGCGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTa                                     <  3:186501/147‑1 (MQ=255)
                                                                                                              gtatgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaaa                              >  7:52441/1‑147 (MQ=255)
                                                                                                                 tgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGt                                      <  5:125383/136‑1 (MQ=255)
                                                                                                                 tgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGt                                      >  6:125383/1‑136 (MQ=255)
                                                                                                                  gtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCg                         >  6:182967/1‑148 (MQ=255)
                                                                                                                     ctccctccCTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCCCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTAc                      >  4:202198/1‑148 (MQ=255)
                                                                                                                     cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTAc                      >  7:265965/1‑148 (MQ=255)
                                                                                                                         cTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTa                  >  7:338429/1‑148 (MQ=255)
                                                                                                                         cTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTa                  >  6:229196/1‑148 (MQ=255)
                                                                                                                         cTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTa                  >  7:225087/1‑148 (MQ=255)
                                                                                                                            aCTCTTAACCGACCCCGCACTATTTTTTGAACCGGCGCCCCTCGGGGTCGTGGGCGGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTAc                 <  7:231234/146‑1 (MQ=255)
                                                                                                                            aCTCTGAAACGACCCCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGa               <  4:143771/148‑1 (MQ=255)
                                                                                                                             ctctGAAACGACACCGCACTCTTTTTTTCTCCCTTGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGaa              >  6:153098/1‑148 (MQ=255)
                                                                                                                             ctctGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGaa              <  2:117682/148‑1 (MQ=255)
                                                                                                                              tctGAACCGACACCGCACTCTTTTTTTCTCACTCGCCCCTCCGGGGAGAGTGCCGGGCAGAGCGGAAAAGGCCGCACTGCTCACCAATTTTCTGATATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCgaa                       >  1:199920/1‑136 (MQ=255)
                                                                                                                               ctGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaa                                                                                             <  1:206881/67‑1 (MQ=255)
                                                                                                                               ctGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaa                                                                                             >  2:206881/1‑67 (MQ=255)
                                                                                                                                       acacCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCCCCAATTTTCTGCTATCCTGCCCCCGCACTATGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTc    >  2:28963/1‑148 (MQ=255)
                                                                                                                                         acCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCg  <  8:265965/148‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                   
TAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCG  >  NC_000913/3780966‑3781250

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