Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 272,330:1 +G 5.7% intergenic (+141/‑250) insN → / → eyeA pseudogene, IS911 transposase A;IS, phage, Tn; Transposon‑related functions; extrachromosomal; transposon related/novel sRNA, function unknown, CP4‑6 propahge

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913272,3301.G 5.7% 768.8 / 39.8 174intergenic (+141/‑250)insN/eyeApseudogene, IS911 transposase A;IS, phage, Tn; Transposon‑related functions; extrachromosomal; transposon related/novel sRNA, function unknown, CP4‑6 propahge
Reads supporting (aligned to +/- strand):  ref base . (75/89);  new base G (8/2);  total (83/91)
Fisher's exact test for biased strand distribution p-value = 4.93e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.28e-01

Only 100 of 180 total aligned reads displayed.

TGACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGTTGAAAACGTCACCACG  >  NC_000913/272182‑272475
                                                                                                                                                     |                                                                                                                                                 
 gACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑g                                                                                                                                                  <  4:1014781/149‑1 (MQ=255)
    cTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCa                                                                                                                                               >  3:846545/1‑149 (MQ=255)
    cTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCa                                                                                                                                               >  1:509379/1‑149 (MQ=255)
     tataGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCAt                                                                                                                                              >  1:951336/1‑149 (MQ=255)
            tgCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCa                                                                                                                                       <  2:255845/149‑1 (MQ=255)
             gCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAg                                                                                                                                      <  3:989094/149‑1 (MQ=255)
              cccGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAgg                                                                                                                                     <  4:435485/149‑1 (MQ=255)
              cccGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAgg                                                                                                                                     <  4:863318/149‑1 (MQ=255)
               ccGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCAt                                                                                                                                              >  4:641989/1‑139 (MQ=255)
               ccGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAggg                                                                                                                                    >  4:749726/1‑149 (MQ=255)
                 ggAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTg                                                                                                                                  <  2:823582/149‑1 (MQ=255)
                   aGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGAt                                                                                                                                <  4:1067478/149‑1 (MQ=255)
                   aGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTATGGAATATTGAACCGGTAACTCACGATGAGAAAACCAACAATCCCACCGCGTGTGACGGTGGAGAACCTGAGCGGCAGTCACCTGC‑GGCATGCCCGCAGGGTGAt                                                                                                                                <  4:672151/149‑1 (MQ=255)
                      tCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTa                                                                                                                             >  1:842391/1‑149 (MQ=255)
                        aGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAAc                                                                                                                           >  4:115738/1‑149 (MQ=255)
                        aGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAAc                                                                                                                           >  1:862926/1‑149 (MQ=255)
                            cgggcATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGc                                                                                                                       <  2:1083852/149‑1 (MQ=255)
                             gggcATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCt                                                                                                                      >  4:814871/1‑149 (MQ=255)
                             gggcATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCt                                                                                                                      <  2:49155/149‑1 (MQ=255)
                                   ggATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAAc                                                                                                                >  1:283624/1‑149 (MQ=255)
                                     aTGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACgg                                                                                                              <  3:577136/149‑1 (MQ=255)
                                     aTGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACgg                                                                                                              <  1:1047485/149‑1 (MQ=255)
                                         ttAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑g                                                                                                                                                  >  2:270452/1‑109 (MQ=255)
                                         ttAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGAt                                                                                                          >  3:849638/1‑149 (MQ=255)
                                         ttAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGAAGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGTAGGGTGATGTAACCCGCTGACAACGGGGAt                                                                                                          >  1:713467/1‑149 (MQ=255)
                                            aaTGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGa                                                                                                       <  1:477408/149‑1 (MQ=255)
                                            aaTGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGa                                                                                                       >  3:652188/1‑149 (MQ=255)
                                             aTGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAg                                                                                                      <  3:115738/149‑1 (MQ=255)
                                             aTGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAg                                                                                                      <  1:341781/149‑1 (MQ=255)
                                             aTGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAg                                                                                                      <  3:724717/149‑1 (MQ=255)
                                               gAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGc                                                                                                    >  2:658524/1‑149 (MQ=255)
                                               gAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGc                                                                                                    <  1:222944/149‑1 (MQ=255)
                                                  ccGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCgag                                                                                                 <  4:199436/149‑1 (MQ=255)
                                                    gcgAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGAt                                                                                               <  3:748282/149‑1 (MQ=255)
                                                       aGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCAc                                                                                            <  1:749785/149‑1 (MQ=255)
                                                       aGTCTGTCTGGAATATTGAACCGGTAACTCACGACGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCAc                                                                                            <  1:50115/149‑1 (MQ=255)
                                                        gtctgtctGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACt                                                                                           <  2:283624/149‑1 (MQ=255)
                                                         tctgtctgGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTa                                                                                          >  3:281246/1‑149 (MQ=255)
                                                               tgGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGAt                                                                                                          >  3:1021963/1‑127 (MQ=255)
                                                               tgGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCga                                                                                    >  2:694362/1‑149 (MQ=255)
                                                               tgGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCga                                                                                    <  3:815627/149‑1 (MQ=255)
                                                               tgGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCga                                                                                    <  3:237847/149‑1 (MQ=255)
                                                                ggAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCgag                                                                                   >  4:952154/1‑149 (MQ=255)
                                                                 gAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGAt                                                                                                          <  2:397154/125‑1 (MQ=255)
                                                                 gAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGAt                                                                                                          >  1:397154/1‑125 (MQ=255)
                                                                    tatTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAgatga                                                                               <  1:455715/149‑1 (MQ=255)
                                                                      ttGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATc                                                                             >  2:614771/1‑149 (MQ=255)
                                                                      ttGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATc                                                                             <  3:752386/149‑1 (MQ=255)
                                                                        gAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATc                                                                             <  1:1048167/147‑1 (MQ=255)
                                                                         aaCCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTc                                                                          >  1:681130/1‑149 (MQ=255)
                                                                                  cTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAg                                                                 <  3:198981/149‑1 (MQ=255)
                                                                                   tCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCgaga                                                                                                <  4:468060/117‑1 (MQ=255)
                                                                                   tCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCgaga                                                                                                >  3:468060/1‑117 (MQ=255)
                                                                                   tCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGt                                                                >  2:818348/1‑149 (MQ=255)
                                                                                          gagaAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGaaa                                                         <  2:456257/149‑1 (MQ=255)
                                                                                          gagaAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGaaa                                                         <  3:391770/149‑1 (MQ=255)
                                                                                           agaAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAg                                                        <  3:634364/149‑1 (MQ=255)
                                                                                             aaaCCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGt                                                      <  1:694362/149‑1 (MQ=255)
                                                                                             aaaCCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGt                                                      <  2:879603/149‑1 (MQ=255)
                                                                                             aaaCACAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGt                                                      <  1:543114/149‑1 (MQ=255)
                                                                                              aaCCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGtt                                                     <  2:951336/149‑1 (MQ=255)
                                                                                               aCCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGaa                                                          <  1:620400/143‑1 (MQ=255)
                                                                                               aCCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGaa                                                          >  2:620400/1‑143 (MQ=255)
                                                                                               aCCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTg                                                    >  4:807210/1‑149 (MQ=255)
                                                                                                cccAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTga                                                   >  3:274211/1‑149 (MQ=255)
                                                                                                 ccAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCga                                                                                    >  1:929988/1‑115 (MQ=255)
                                                                                                 ccAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCga                                                                                    <  2:929988/115‑1 (MQ=255)
                                                                                                 ccAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATc                                                                             >  2:919805/1‑122 (MQ=255)
                                                                                                 ccAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATc                                                                             <  1:919805/122‑1 (MQ=255)
                                                                                                 ccAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAg                                                        <  4:497332/143‑1 (MQ=255)
                                                                                                 ccAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTgaa                                                  >  2:559879/1‑149 (MQ=255)
                                                                                                  caacaaTCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTgaag                                                 >  1:951285/1‑149 (MQ=255)
                                                                                                   aacaaTCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTgaaga                                                >  2:384183/1‑149 (MQ=255)
                                                                                                   aacaaTCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTgaaga                                                >  3:263555/1‑149 (MQ=255)
                                                                                                     caaTCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAAc                                              <  3:805247/149‑1 (MQ=255)
                                                                                                       aTCCCACCGGGTGTGACGGTGGAGAACCTGATCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCCCTGACAACGGGGAATGAGGCGAGATCACTAAGCCGACATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACAt                                            <  3:667396/149‑1 (MQ=255)
                                                                                                          ccACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGaa                                         >  4:1029981/1‑149 (MQ=255)
                                                                                                           cACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAAc                                        >  2:543644/1‑149 (MQ=255)
                                                                                                              cGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCg                                     >  1:966035/1‑149 (MQ=255)
                                                                                                               gggTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAgg                                                                      >  3:728858/1‑115 (MQ=255)
                                                                                                               gggTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAgg                                                                      <  4:728858/115‑1 (MQ=255)
                                                                                                               gggTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGt                                    <  1:256654/149‑1 (MQ=255)
                                                                                                               gggTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGt                                    <  1:906048/149‑1 (MQ=255)
                                                                                                               gggTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGt                                    <  1:551982/149‑1 (MQ=255)
                                                                                                                   gtgACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAAt                                <  4:419095/149‑1 (MQ=255)
                                                                                                                     gACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATcc                              >  3:637829/1‑149 (MQ=255)
                                                                                                                             agaACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCtgt                      <  4:184289/149‑1 (MQ=255)
                                                                                                                              gaACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCtgtg                     >  1:343408/1‑149 (MQ=255)
                                                                                                                                 ccTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGt                  <  2:745991/149‑1 (MQ=255)
                                                                                                                                  cTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTgaa                                                  <  2:359720/116‑1 (MQ=255)
                                                                                                                                  cTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGtt                 >  3:594204/1‑149 (MQ=255)
                                                                                                                                   tGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGTTg                >  3:460143/1‑149 (MQ=255)
                                                                                                                                            gTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGTTGAAAACGTca       <  4:846545/149‑1 (MQ=255)
                                                                                                                                             cacccTGCGGGCATG‑CCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGTTGAAAACGTcac      >  2:429713/4‑149 (MQ=255)
                                                                                                                                              acccTGCGGGCATG‑CCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTgg                    >  1:658576/3‑134 (MQ=255)
                                                                                                                                                ccTGCGGGCATG‑CCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGAt                                                                                               <  4:407093/57‑1 (MQ=37)
                                                                                                                                                ccTGCGGGCATG‑CCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGAt                                                                                               >  3:407093/1‑57 (MQ=37)
                                                                                                                                                ccTGCGGGCATG‑CCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGTTGAAAACGTcaccac   >  4:254667/1‑149 (MQ=255)
                                                                                                                                                 cTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATc                               <  3:523395/120‑1 (MQ=255)
                                                                                                                                                 cTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGTTGAAAACGTCACCACg  >  4:604656/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
TGACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGTTGAAAACGTCACCACG  >  NC_000913/272182‑272475

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑