Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 272,330:1 +G 7.4% 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 7.4% 337.7 / 23.0 81intergenic (+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 . (40/35);  new base G (3/3);  total (43/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

ACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGTTGAAAACGTCACCACGGC  >  NC_000913/272184‑272477
                                                                                                                                                   |                                                                                                                                                   
aCCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑gg                                                                                                                                                   >  4:310185/1‑149 (MQ=255)
  cTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCa                                                                                                                                                 >  2:187580/1‑149 (MQ=255)
  cTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAAACCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCa                                                                                                                                                 <  1:474857/149‑1 (MQ=255)
     taGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGc                                                                                                                                              <  3:165238/149‑1 (MQ=255)
           gCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAg                                                                                                                                        <  1:187580/149‑1 (MQ=255)
           gCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAg                                                                                                                                        <  3:564161/149‑1 (MQ=255)
            cccGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGGCATG‑CCGCAgg                                                                                                                                       >  1:354053/1‑149 (MQ=255)
            cccGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAgg                                                                                                                                       >  3:502715/1‑149 (MQ=255)
             ccGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGGCATG‑CCGCAggt                                                                                                                                      <  3:179631/149‑2 (MQ=255)
              cGGAGTTCGGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCATCAATCCCTCCGGTTGTGACGGTGGAGATCCTGAGCGGCATTGACCTGC‑GGCATGCCCCCAGGGt                                                                                                                                     <  4:283358/149‑1 (MQ=255)
              cGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGt                                                                                                                                     <  3:131449/149‑1 (MQ=255)
               ggAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGGCATG‑CCGCAggtca                                                                                                                                    <  2:354053/149‑4 (MQ=255)
                                      cTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGa                                                                                                             >  1:504091/1‑149 (MQ=255)
                                       ttAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGAt                                                                                                            >  3:22711/1‑149 (MQ=255)
                                           aTGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAg                                                                                                        >  4:350426/1‑149 (MQ=255)
                                            tGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAgg                                                                                                       >  1:441501/1‑149 (MQ=255)
                                             gAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGAGATGTAACCCGCTGACAACGGGGATTGAGGc                                                                                                      >  2:294594/1‑149 (MQ=255)
                                               aCCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCga                                                                                                    >  3:13275/1‑149 (MQ=255)
                                                ccGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCgag                                                                                                   >  2:584724/1‑149 (MQ=255)
                                                ccGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCga                                                                                                    >  3:66912/1‑148 (MQ=255)
                                                ccGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCga                                                                                                    <  4:66912/148‑1 (MQ=255)
                                                    gAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGATCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCa                                                                                               >  4:347658/1‑149 (MQ=255)
                                                    gAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCa                                                                                               >  2:229017/1‑149 (MQ=255)
                                                    gAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCa                                                                                               <  3:45146/149‑1 (MQ=255)
                                                               gAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCgaga                                                                                    >  2:303328/1‑149 (MQ=255)
                                                                        aCCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCa                                                                           >  3:337942/1‑149 (MQ=255)
                                                                        aCCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCa                                                                           >  3:107278/1‑149 (MQ=255)
                                                                          cGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAg                                                                         <  4:495063/149‑1 (MQ=255)
                                                                            gTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACTGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGt                                                                       <  3:488353/149‑1 (MQ=255)
                                                                            gTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGt                                                                       <  1:78151/149‑1 (MQ=255)
                                                                               aCTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTaa                                                                    <  3:257125/149‑1 (MQ=255)
                                                                                 tCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCgag                                                                                                   <  3:92031/116‑1 (MQ=255)
                                                                                 tCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCgag                                                                                                   >  4:92031/1‑116 (MQ=255)
                                                                                 tCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGt                                                                  >  1:388194/1‑149 (MQ=255)
                                                                                 tCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGt                                                                  >  3:336437/1‑149 (MQ=255)
                                                                                  cACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTa                                                                 <  4:270912/149‑1 (MQ=255)
                                                                                   aCGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTAc                                                                >  2:356656/1‑149 (MQ=255)
                                                                                   aCGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTAc                                                                >  3:265252/1‑149 (MQ=255)
                                                                                     gATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTAc                                                                >  4:248141/1‑147 (MQ=255)
                                                                                     gATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTAc                                                                <  3:248141/147‑1 (MQ=255)
                                                                                     gATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTg                                                              >  2:182359/1‑149 (MQ=255)
                                                                                         agaAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAg                                                          <  2:47687/149‑1 (MQ=255)
                                                                                           aaaCCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGt                                                        >  1:493911/1‑149 (MQ=255)
                                                                                           aaaCCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGt                                                        >  1:75584/1‑149 (MQ=255)
                                                                                            accccAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGtt                                                       <  2:75584/147‑1 (MQ=255)
                                                                                             aCCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTg                                                      <  1:184784/149‑1 (MQ=255)
                                                                                             aCCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTg                                                      >  4:365842/1‑149 (MQ=255)
                                                                                               ccAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGtt                                                       <  3:292828/146‑1 (MQ=255)
                                                                                               ccAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGtt                                                       >  4:292828/1‑146 (MQ=255)
                                                                                               ccAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTgaa                                                    >  1:435838/1‑149 (MQ=255)
                                                                                               ccAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTgaa                                                    <  2:242527/149‑1 (MQ=255)
                                                                                               ccAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTgaa                                                    <  1:51054/149‑1 (MQ=255)
                                                                                               ccAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCACAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTgaa                                                    <  1:41791/149‑1 (MQ=255)
                                                                                                caacaaTCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTgaag                                                   >  4:331530/1‑149 (MQ=255)
                                                                                                caacaaTCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTgaag                                                   >  1:56203/1‑149 (MQ=255)
                                                                                                     aTCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGTGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACAt                                              >  4:130603/1‑149 (MQ=255)
                                                                                                      tCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATg                                             <  4:32159/149‑1 (MQ=255)
                                                                                                        ccACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGaa                                           <  4:13275/149‑1 (MQ=255)
                                                                                                            cGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCg                                       <  2:435838/149‑1 (MQ=255)
                                                                                                                      ggTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGcc                             >  3:400945/1‑149 (MQ=255)
                                                                                                                       gTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAgatga                                                                                 <  3:509505/96‑1 (MQ=255)
                                                                                                                       gTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAgatga                                                                                 >  4:509505/1‑96 (MQ=255)
                                                                                                                          gagaACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCtg                         >  1:237138/1‑149 (MQ=255)
                                                                                                                          gagaACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCtg                         >  3:586508/1‑149 (MQ=255)
                                                                                                                          gagaACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGCGGATTGAGGCGAGATCACTAAGCCGAGATGATACTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCtg                         >  4:557920/1‑149 (MQ=255)
                                                                                                                            gaACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCtgt                        >  2:284007/1‑148 (MQ=255)
                                                                                                                            gaACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCtgt                        <  1:284007/148‑1 (MQ=255)
                                                                                                                             aaCCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTgg                      <  3:350426/149‑1 (MQ=255)
                                                                                                                                cTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGtt                   <  3:365842/149‑1 (MQ=255)
                                                                                                                                   aGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGTTGaa                <  2:484522/149‑1 (MQ=255)
                                                                                                                                   aGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGTTGaa                <  4:400945/149‑1 (MQ=255)
                                                                                                                                     cGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGTTGaaaa              <  1:229017/149‑1 (MQ=255)
                                                                                                                                          gTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGTTGAAAACGTca         <  4:22711/149‑1 (MQ=255)
                                                                                                                                              ccTGCGGGCATG‑CCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGTTGaa                <  1:463665/138‑1 (MQ=255)
                                                                                                                                              ccTGCGGGCATG‑CCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGTTGaa                >  2:463665/1‑138 (MQ=255)
                                                                                                                                              ccTGCGGGCATG‑CCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGTTGAAAACGTcaccac     >  4:286823/1‑149 (MQ=255)
                                                                                                                                              ccTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGtt                                                       >  3:523312/1‑99 (MQ=255)
                                                                                                                                              ccTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGtt                                                       <  4:523312/99‑1 (MQ=255)
                                                                                                                                              ccTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATcc                                >  3:135101/1‑122 (MQ=255)
                                                                                                                                              ccTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATcc                                <  4:135101/122‑1 (MQ=255)
                                                                                                                                                tGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGTTGAAAACGTCACCACgg   >  1:207164/1‑149 (MQ=255)
                                                                                                                                                 gC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGTTGAAAACGTCACCACGGc  >  3:396159/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGTTGAAAACGTCACCACGGC  >  NC_000913/272184‑272477

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