Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 272,330:1 +G 6.2% 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 6.2% 426.2 / 21.8 97intergenic (+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 . (44/47);  new base G (4/2);  total (48/49)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

ACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGTTGAAAACGTCACCA  >  NC_000913/272184‑272473
                                                                                                                                                   |                                                                                                                                               
aCCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑gg                                                                                                                                               >  1:15222/1‑149 (MQ=255)
  cTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCa                                                                                                                                             >  4:398319/1‑149 (MQ=255)
   tataGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCAt                                                                                                                                            >  1:176200/1‑149 (MQ=255)
     taGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GACATGc                                                                                                                                          >  2:316662/1‑149 (MQ=255)
       gtgtGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGccc                                                                                                                                        >  1:40567/1‑149 (MQ=255)
             ccGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAggg                                                                                                                                  >  3:26081/1‑149 (MQ=255)
               ggAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTg                                                                                                                                >  3:133168/1‑149 (MQ=255)
                       gggcgggcATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAAcc                                                                                                                        >  1:451452/1‑149 (MQ=255)
                             gcATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGa                                                                                                                  <  1:368033/149‑1 (MQ=255)
                                  gATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACg                                                                                                             >  3:258982/1‑149 (MQ=255)
                                    tGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GACATGCCCGCAGGGTGATGTAACCCGCTGACAACggg                                                                                                           <  1:316662/149‑1 (MQ=255)
                                     gCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACgggg                                                                                                          >  1:103669/1‑149 (MQ=255)
                                       ttAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAgg                                                                                                                                   <  1:240265/122‑1 (MQ=255)
                                       ttAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAgg                                                                                                                                   >  2:240265/1‑122 (MQ=255)
                                       ttAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACgg                                                                                                            <  4:86904/145‑1 (MQ=255)
                                       ttAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACgg                                                                                                            >  3:86904/1‑145 (MQ=255)
                                         aaaTGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTg                                                                                                      >  3:994/1‑149 (MQ=255)
                                               aCCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCga                                                                                                >  3:314566/1‑149 (MQ=255)
                                                 cgcgAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAgg                                                                                                   <  2:85756/144‑1 (MQ=255)
                                                 cgcgAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAgg                                                                                                   >  1:85756/1‑144 (MQ=255)
                                                  gcgAGTCTGTCTGGAACATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGAt                                                                                             <  4:26081/149‑1 (MQ=255)
                                                   cgAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATc                                                                                            <  2:103669/149‑1 (MQ=255)
                                                     aGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCAc                                                                                          <  4:273372/149‑1 (MQ=255)
                                                      gtctgtctGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACt                                                                                         <  2:107692/149‑1 (MQ=255)
                                                        ctgtctgGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTaa                                                                                       <  3:63359/149‑1 (MQ=255)
                                                        ctgtctgGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTaa                                                                                       <  3:34223/149‑1 (MQ=255)
                                                         tgtctgGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAg                                                                                      >  1:91996/1‑149 (MQ=255)
                                                             tgGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGa                                                                                                         <  1:19980/126‑1 (MQ=255)
                                                             tgGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGa                                                                                                         >  2:19980/1‑126 (MQ=255)
                                                             tgGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCga                                                                                  >  2:415395/1‑149 (MQ=255)
                                                             tgGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCga                                                                                  <  2:423441/149‑1 (MQ=255)
                                                             tgGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCga                                                                                  >  1:287627/1‑149 (MQ=255)
                                                             tgGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCga                                                                                  <  3:395630/149‑1 (MQ=255)
                                                             tgGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCga                                                                                  >  1:66938/1‑149 (MQ=255)
                                                             tgGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCga                                                                                  >  1:50812/1‑149 (MQ=255)
                                                              ggAATATTGAACCGGTAACTCACGATGAGAAACCCAACAACCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCgag                                                                                 >  4:276711/1‑149 (MQ=255)
                                                               gAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCgaga                                                                                >  3:150132/1‑149 (MQ=255)
                                                                  tatTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAgatga                                                                             >  1:330769/1‑149 (MQ=255)
                                                                      gAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCt                                                                         >  3:286583/1‑149 (MQ=255)
                                                                       aaCCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGAt                                                                                                        >  4:274044/1‑117 (MQ=255)
                                                                       aaCCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGAt                                                                                                        <  3:274044/117‑1 (MQ=255)
                                                                         ccGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCaa                                                                      <  2:343098/149‑1 (MQ=255)
                                                                         ccGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCaa                                                                      <  3:405314/149‑1 (MQ=255)
                                                                          cGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTCACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAg                                                                     <  4:286583/149‑1 (MQ=255)
                                                                            gTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGt                                                                   >  4:147920/1‑149 (MQ=255)
                                                                               aCTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTaa                                                                <  3:132832/149‑1 (MQ=255)
                                                                                cTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAg                                                               <  3:6956/149‑1 (MQ=255)
                                                                                cTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAg                                                               <  4:150132/149‑1 (MQ=255)
                                                                                 tCACGATTAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGt                                                              >  2:259923/1‑149 (MQ=255)
                                                                                 tCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGt                                                              >  3:331463/1‑149 (MQ=255)
                                                                                 tCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGt                                                              >  3:307077/1‑149 (MQ=255)
                                                                                 tCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGt                                                              >  2:327157/1‑149 (MQ=255)
                                                                                 tCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGt                                                              >  1:221283/1‑149 (MQ=255)
                                                                                 tCACGATGAGAAACCCAACAATCCAACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCAGTAAGACGAGATGATCCTCGAGGTTAAGt                                                              >  3:298406/1‑149 (MQ=255)
                                                                                     gATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGc                                                                                                                                                  <  1:213915/62‑1 (MQ=255)
                                                                                     gATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGc                                                                                                                                                  >  2:213915/1‑62 (MQ=255)
                                                                                         agaAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAg                                                      >  2:384940/1‑149 (MQ=255)
                                                                                          gaAACCCAACACTCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATACTCAAGGTTAAGTACTGAAAgg                                                     >  3:271779/1‑149 (MQ=255)
                                                                                          gaAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAgg                                                     <  2:20367/149‑1 (MQ=255)
                                                                                          gaAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAgg                                                     <  4:258982/149‑1 (MQ=255)
                                                                                          gaAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAgg                                                     >  2:161142/1‑149 (MQ=255)
                                                                                           aaaCCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGt                                                    >  4:193302/1‑149 (MQ=255)
                                                                                            aaCCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGACCCTCAAGGTTAAGTACTGAAAGGtt                                                   <  4:447744/149‑1 (MQ=255)
                                                                                             aCCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTg                                                  >  4:357847/1‑149 (MQ=255)
                                                                                             aCCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTg                                                  <  4:472251/149‑1 (MQ=255)
                                                                                             aCCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTg                                                  <  2:287627/149‑1 (MQ=255)
                                                                                             aCCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTg                                                  <  4:236299/149‑1 (MQ=255)
                                                                                              cccAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTga                                                 <  1:283700/149‑1 (MQ=255)
                                                                                               ccAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTgaa                                                <  3:339502/149‑1 (MQ=255)
                                                                                               ccAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTgaa                                                <  2:176200/149‑1 (MQ=255)
                                                                                               ccAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTgaa                                                >  4:16658/1‑149 (MQ=255)
                                                                                               ccAACAATCCCACCGGGCGTGACGGTGGAGAACCTGAGCGGCAGTGACCAGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCTGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTgaa                                                <  4:298406/149‑1 (MQ=255)
                                                                                                  acaaTCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTgaagaa                                             <  3:357847/149‑1 (MQ=255)
                                                                                                    aaTCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACa                                           <  4:133168/149‑1 (MQ=255)
                                                                                                    aaTCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACa                                           <  2:91996/149‑1 (MQ=255)
                                                                                                        ccACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGaa                                       <  1:415395/149‑1 (MQ=255)
                                                                                                         cACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGt                                                    >  2:248059/1‑135 (MQ=255)
                                                                                                         cACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGt                                                    <  1:248059/135‑1 (MQ=255)
                                                                                                           ccGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAAccc                                    <  3:398319/149‑1 (MQ=255)
                                                                                                           ccGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAAccc                                    <  2:66938/149‑1 (MQ=255)
                                                                                                                     cGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGc                          <  2:40567/149‑1 (MQ=255)
                                                                                                                         ggAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCctct                      <  4:331463/149‑1 (MQ=255)
                                                                                                                            gaACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCtgtg                   >  3:166454/1‑149 (MQ=255)
                                                                                                                             aaCCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTgg                  <  3:87602/149‑1 (MQ=255)
                                                                                                                               ccTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGt                >  2:450427/1‑149 (MQ=255)
                                                                                                                               ccTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGt                <  3:276711/149‑1 (MQ=255)
                                                                                                                               ccTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGt                <  2:330769/149‑1 (MQ=255)
                                                                                                                                cTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACgg                                                                                                            <  3:347039/56‑1 (MQ=255)
                                                                                                                                cTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACgg                                                                                                            >  4:347039/1‑56 (MQ=255)
                                                                                                                                      ggCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGTTGAAAAc         <  3:16658/149‑1 (MQ=255)
                                                                                                                                        cAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGTTGAAAACGt       <  2:451452/149‑1 (MQ=255)
                                                                                                                                         aGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGTTGAAAACGTc      >  1:55914/1‑149 (MQ=255)
                                                                                                                                           cacccTGCGGGCATG‑CCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATc                             <  1:193104/121‑1 (MQ=255)
                                                                                                                                           cacccTGCGGGCATG‑CCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATc                             >  2:193104/4‑124 (MQ=255)
                                                                                                                                             cccTGCGGGCATG‑CCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGTTGAAAACGTcacca  >  3:121693/2‑149 (MQ=255)
                                                                                                                                             cccTGCGGGCATG‑CCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGTTGAAAACGTcacca  >  1:70236/2‑149 (MQ=255)
                                                                                                                                              ccTGCGGGCATG‑CCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGcc                         <  3:133114/125‑1 (MQ=255)
                                                                                                                                              ccTGCGGGCATG‑CCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGcc                         >  4:133114/1‑125 (MQ=255)
                                                                                                                                               cTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATcc                                                                          <  4:467738/75‑1 (MQ=255)
                                                                                                                                               cTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATcc                                                                          >  3:467738/1‑75 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
ACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGTTGAAAACGTCACCA  >  NC_000913/272184‑272473

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