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

ACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGTTGAAAACGTCACCA  >  NC_000913/272184‑272473
                                                                                                                                                   |                                                                                                                                               
aCCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCgg                                                                                                                                               >  1:15222/1‑149 (MQ=255)
  cTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCa                                                                                                                                             >  4:398319/1‑149 (MQ=255)
   tataGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCAt                                                                                                                                            >  1:176200/1‑149 (MQ=255)
     taGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGACATGc                                                                                                                                          >  2:316662/1‑149 (MQ=255)
       gtgtGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGccc                                                                                                                                        >  1:40567/1‑149 (MQ=255)
             ccGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAggg                                                                                                                                  >  3:26081/1‑149 (MQ=255)
               ggAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTg                                                                                                                                >  3:133168/1‑149 (MQ=255)
                       gggcgggcATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAAcc                                                                                                                        >  1:451452/1‑149 (MQ=255)
                             gcATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGa                                                                                                                  <  1:368033/149‑1 (MQ=255)
                                  gATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACg                                                                                                             >  3:258982/1‑149 (MQ=255)
                                    tGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGACATGCCCGCAGGGTGATGTAACCCGCTGACAACggg                                                                                                           <  1:316662/149‑1 (MQ=255)
                                     gCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACgggg                                                                                                          >  1:103669/1‑149 (MQ=255)
                                       ttAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAgg                                                                                                                                   <  1:240265/122‑1 (MQ=255)
                                       ttAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAgg                                                                                                                                   >  2:240265/1‑122 (MQ=255)
                                       ttAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACgg                                                                                                            <  4:86904/145‑1 (MQ=255)
                                       ttAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACgg                                                                                                            >  3:86904/1‑145 (MQ=255)
                                         aaaTGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTg                                                                                                      >  3:994/1‑149 (MQ=255)
                                               aCCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCga                                                                                                >  3:314566/1‑149 (MQ=255)
                                                 cgcgAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAgg                                                                                                   <  2:85756/144‑1 (MQ=255)
                                                 cgcgAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAgg                                                                                                   >  1:85756/1‑144 (MQ=255)
                                                  gcgAGTCTGTCTGGAACATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGAt                                                                                             <  4:26081/149‑1 (MQ=255)
                                                   cgAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATc                                                                                            <  2:103669/149‑1 (MQ=255)
                                                     aGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCAc                                                                                          <  4:273372/149‑1 (MQ=255)
                                                      gtctgtctGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACt                                                                                         <  2:107692/149‑1 (MQ=255)
                                                        ctgtctgGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTaa                                                                                       <  3:63359/149‑1 (MQ=255)
                                                        ctgtctgGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTaa                                                                                       <  3:34223/149‑1 (MQ=255)
                                                         tgtctgGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAg                                                                                      >  1:91996/1‑149 (MQ=255)
                                                             tgGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGa                                                                                                         <  1:19980/126‑1 (MQ=255)
                                                             tgGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGa                                                                                                         >  2:19980/1‑126 (MQ=255)
                                                             tgGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCga                                                                                  >  2:415395/1‑149 (MQ=255)
                                                             tgGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCga                                                                                  <  2:423441/149‑1 (MQ=255)
                                                             tgGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCga                                                                                  >  1:287627/1‑149 (MQ=255)
                                                             tgGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCga                                                                                  <  3:395630/149‑1 (MQ=255)
                                                             tgGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCga                                                                                  >  1:66938/1‑149 (MQ=255)
                                                             tgGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCga                                                                                  >  1:50812/1‑149 (MQ=255)
                                                              ggAATATTGAACCGGTAACTCACGATGAGAAACCCAACAACCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCgag                                                                                 >  4:276711/1‑149 (MQ=255)
                                                               gAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCgaga                                                                                >  3:150132/1‑149 (MQ=255)
                                                                  tatTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAgatga                                                                             >  1:330769/1‑149 (MQ=255)
                                                                      gAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCt                                                                         >  3:286583/1‑149 (MQ=255)
                                                                       aaCCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGAt                                                                                                        >  4:274044/1‑117 (MQ=255)
                                                                       aaCCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGAt                                                                                                        <  3:274044/117‑1 (MQ=255)
                                                                         ccGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCaa                                                                      <  2:343098/149‑1 (MQ=255)
                                                                         ccGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCaa                                                                      <  3:405314/149‑1 (MQ=255)
                                                                          cGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTCACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAg                                                                     <  4:286583/149‑1 (MQ=255)
                                                                            gTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGt                                                                   >  4:147920/1‑149 (MQ=255)
                                                                               aCTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTaa                                                                <  3:132832/149‑1 (MQ=255)
                                                                                cTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAg                                                               <  3:6956/149‑1 (MQ=255)
                                                                                cTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAg                                                               <  4:150132/149‑1 (MQ=255)
                                                                                 tCACGATTAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGt                                                              >  2:259923/1‑149 (MQ=255)
                                                                                 tCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGt                                                              >  3:331463/1‑149 (MQ=255)
                                                                                 tCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGt                                                              >  3:307077/1‑149 (MQ=255)
                                                                                 tCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGt                                                              >  2:327157/1‑149 (MQ=255)
                                                                                 tCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGt                                                              >  1:221283/1‑149 (MQ=255)
                                                                                 tCACGATGAGAAACCCAACAATCCAACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCAGTAAGACGAGATGATCCTCGAGGTTAAGt                                                              >  3:298406/1‑149 (MQ=255)
                                                                                     gATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGc                                                                                                                                                  <  1:213915/62‑1 (MQ=255)
                                                                                     gATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGc                                                                                                                                                  >  2:213915/1‑62 (MQ=255)
                                                                                         agaAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAg                                                      >  2:384940/1‑149 (MQ=255)
                                                                                          gaAACCCAACACTCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATACTCAAGGTTAAGTACTGAAAgg                                                     >  3:271779/1‑149 (MQ=255)
                                                                                          gaAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAgg                                                     <  2:20367/149‑1 (MQ=255)
                                                                                          gaAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAgg                                                     <  4:258982/149‑1 (MQ=255)
                                                                                          gaAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAgg                                                     >  2:161142/1‑149 (MQ=255)
                                                                                           aaaCCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGt                                                    >  4:193302/1‑149 (MQ=255)
                                                                                            aaCCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGACCCTCAAGGTTAAGTACTGAAAGGtt                                                   <  4:447744/149‑1 (MQ=255)
                                                                                             aCCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTg                                                  >  4:357847/1‑149 (MQ=255)
                                                                                             aCCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTg                                                  <  4:472251/149‑1 (MQ=255)
                                                                                             aCCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTg                                                  <  2:287627/149‑1 (MQ=255)
                                                                                             aCCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTg                                                  <  4:236299/149‑1 (MQ=255)
                                                                                              cccAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTga                                                 <  1:283700/149‑1 (MQ=255)
                                                                                               ccAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTgaa                                                <  3:339502/149‑1 (MQ=255)
                                                                                               ccAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTgaa                                                <  2:176200/149‑1 (MQ=255)
                                                                                               ccAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTgaa                                                >  4:16658/1‑149 (MQ=255)
                                                                                               ccAACAATCCCACCGGGCGTGACGGTGGAGAACCTGAGCGGCAGTGACCAGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCTGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTgaa                                                <  4:298406/149‑1 (MQ=255)
                                                                                                  acaaTCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTgaagaa                                             <  3:357847/149‑1 (MQ=255)
                                                                                                    aaTCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACa                                           <  4:133168/149‑1 (MQ=255)
                                                                                                    aaTCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACa                                           <  2:91996/149‑1 (MQ=255)
                                                                                                        ccACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGaa                                       <  1:415395/149‑1 (MQ=255)
                                                                                                         cACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGt                                                    >  2:248059/1‑135 (MQ=255)
                                                                                                         cACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGt                                                    <  1:248059/135‑1 (MQ=255)
                                                                                                           ccGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAAccc                                    <  3:398319/149‑1 (MQ=255)
                                                                                                           ccGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAAccc                                    <  2:66938/149‑1 (MQ=255)
                                                                                                                     cGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGc                          <  2:40567/149‑1 (MQ=255)
                                                                                                                         ggAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCctct                      <  4:331463/149‑1 (MQ=255)
                                                                                                                            gaACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCtgtg                   >  3:166454/1‑149 (MQ=255)
                                                                                                                             aaCCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTgg                  <  3:87602/149‑1 (MQ=255)
                                                                                                                               ccTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGt                >  2:450427/1‑149 (MQ=255)
                                                                                                                               ccTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGt                <  3:276711/149‑1 (MQ=255)
                                                                                                                               ccTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGt                <  2:330769/149‑1 (MQ=255)
                                                                                                                                cTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACgg                                                                                                            <  3:347039/56‑1 (MQ=255)
                                                                                                                                cTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACgg                                                                                                            >  4:347039/1‑56 (MQ=255)
                                                                                                                                      ggCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGTTGAAAAc         <  3:16658/149‑1 (MQ=255)
                                                                                                                                        cAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGTTGAAAACGt       <  2:451452/149‑1 (MQ=255)
                                                                                                                                         aGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGTTGAAAACGTc      >  1:55914/1‑149 (MQ=255)
                                                                                                                                           cacccTGCGGGCATGCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATc                             <  1:193104/121‑1 (MQ=255)
                                                                                                                                           cacccTGCGGGCATGCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATc                             >  2:193104/4‑124 (MQ=255)
                                                                                                                                             cccTGCGGGCATGCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGTTGAAAACGTcacca  >  3:121693/2‑149 (MQ=255)
                                                                                                                                             cccTGCGGGCATGCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGTTGAAAACGTcacca  >  1:70236/2‑149 (MQ=255)
                                                                                                                                              ccTGCGGGCATGCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGcc                         <  3:133114/125‑1 (MQ=255)
                                                                                                                                              ccTGCGGGCATGCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGcc                         >  4:133114/1‑125 (MQ=255)
                                                                                                                                               cTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATcc                                                                          <  4:467738/75‑1 (MQ=255)
                                                                                                                                               cTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATcc                                                                          >  3:467738/1‑75 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
ACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGTTGAAAACGTCACCA  >  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: