Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 629,579 T→C S94P (TCT→CCT)  entH → proofreading thioesterase in enterobactin biosynthesis

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913629,5790TC97.7% 289.3 / NA 87S94P (TCT→CCT) entHproofreading thioesterase in enterobactin biosynthesis
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (56/29);  minor base G (2/0);  total (58/29)
Fisher's exact test for biased strand distribution p-value = 5.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAA  >  NC_000913/629443‑629714
                                                                                                                                        |                                                                                                                                       
gCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACACTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTg                                                                                                                                      <  2:647232/140‑1 (MQ=255)
  cGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAg                                                                                                                                    <  2:821289/140‑1 (MQ=255)
    ttCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAggg                                                                                                                                  <  2:1092273/140‑1 (MQ=255)
    ttCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAggg                                                                                                                                  >  1:1092273/1‑140 (MQ=255)
    ttCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAggg                                                                                                                                  >  1:1328558/1‑140 (MQ=255)
     tCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGa                                                                                                                                 <  1:57530/140‑1 (MQ=255)
                   ggcggcGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACATACCATCGACCGATGCCTGAGGGAAAGGTACGCGGCGt                                                                                                                   >  1:542441/1‑140 (MQ=255)
                   ggcggcGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGTCCGGTGGCTGAGGGAAAGGTACGCGGCGt                                                                                                                   >  1:598416/1‑140 (MQ=255)
                   ggcggcGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGt                                                                                                                   >  1:42980/1‑140 (MQ=255)
                   ggcggcGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGt                                                                                                                   >  1:1113289/1‑140 (MQ=255)
                     cggcgCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCt                                                                                                                 >  2:1432039/1‑140 (MQ=255)
                     cggcgCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCt                                                                                                                 >  1:1387693/1‑140 (MQ=255)
                        cgcgTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATTCAACACACCATCGCACTGTGCCTGAGGGAATCGTACGCGGCGTCTGcc                                                                                                              >  2:728828/1‑140 (MQ=255)
                        cgcgTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTAcg                                                                                                                         >  2:356630/1‑129 (MQ=255)
                        cgcgTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTAcg                                                                                                                         <  1:356630/129‑1 (MQ=255)
                        cgcgTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGcc                                                                                                              >  2:271155/1‑140 (MQ=255)
                        cgcgTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATACAACACACGATCGCCCAGTGCCTGAGGGAAAGGTACGCGGCGTCTGcc                                                                                                              >  2:422617/1‑140 (MQ=255)
                             cggcggCGCTGGCGGAAACGCTGGGGTCGATGGCCGGTTTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCg                                                                                                         <  1:586256/140‑1 (MQ=255)
                             cggcggCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCg                                                                                                         <  1:140505/140‑1 (MQ=255)
                               gcggcgCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCt                                                                                                       <  2:1387693/140‑1 (MQ=255)
                                    gcTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATAGCCCGGTGGCTGGGGGTAAGGTACGCGGCGACTGCCAGCCGCTGCATc                                                                                                  >  2:586256/1‑140 (MQ=255)
                                             aaCGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGc                                                                                         >  2:1203247/1‑140 (MQ=255)
                                             aaCGCTGGGATCGATAGCCGGACTTCTGACGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGc                                                                                         <  1:728828/140‑1 (MQ=255)
                                                    ggATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCa                                                                                  <  1:271155/140‑1 (MQ=255)
                                                     gATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCag                                                                                 >  1:252299/1‑140 (MQ=255)
                                                      aTCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCaga                                                                                >  2:340343/1‑140 (MQ=255)
                                                       tCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCagag                                                                               <  2:360574/140‑1 (MQ=255)
                                                          aTGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCa                                                                                  >  2:872206/1‑134 (MQ=255)
                                                          aTGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCa                                                                                  <  1:872206/134‑1 (MQ=255)
                                                          aTGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTg                                                                            >  2:913767/1‑140 (MQ=255)
                                                          aTGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTg                                                                            >  1:215254/1‑140 (MQ=255)
                                                          aTGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGCACGCGGCGTCTGACAGACGCTGCATTTTGGTCGGCAAAATCAGAGCTg                                                                            >  2:706862/1‑140 (MQ=255)
                                                              ccGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGaa                                                                        >  1:162474/1‑140 (MQ=255)
                                                              ccGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGaa                                                                        >  1:1053904/1‑140 (MQ=255)
                                                                   tttATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATTCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAAtcgt                                                                   <  2:633574/140‑1 (MQ=255)
                                                                      atgatgACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATcgtcgt                                                                <  2:1427907/140‑1 (MQ=255)
                                                                           gaCCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCg                                                           <  1:1432039/140‑1 (MQ=255)
                                                                            aCCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGa                                                          >  1:854590/1‑140 (MQ=255)
                                                                            aCCCGCGACGGACAGTGTGTGGAAGGCACAGAACTAAATCCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGa                                                          <  2:598416/140‑1 (MQ=255)
                                                                                gcgACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAGAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGaa                                                      <  2:1050429/140‑1 (MQ=255)
                                                                                           tgtgtgGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCTGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGgcggcg                                           >  2:719420/1‑140 (MQ=255)
                                                                                           tgtgtgGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGgcggcg                                           <  2:252299/140‑1 (MQ=255)
                                                                                                 gTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTtgctg                                     <  1:1203247/140‑1 (MQ=255)
                                                                                                    ggCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCAc                                  >  1:1484817/1‑140 (MQ=255)
                                                                                                     gCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGCTGTACAGGGACGGCGTTGCTGCACt                                 >  1:715815/1‑140 (MQ=255)
                                                                                                     gCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGggcg                                              <  1:762202/127‑1 (MQ=255)
                                                                                                     gCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGggcg                                              >  2:762202/1‑127 (MQ=255)
                                                                                                     gCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCa                                   >  1:880377/1‑138 (MQ=255)
                                                                                                     gCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCa                                   <  2:880377/138‑1 (MQ=255)
                                                                                                     gCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACt                                 >  2:1467002/1‑140 (MQ=255)
                                                                                                     gCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACt                                 >  2:903857/1‑140 (MQ=255)
                                                                                                     gCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACt                                 >  2:1408150/1‑140 (MQ=255)
                                                                                                     gCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACt                                 >  1:915513/1‑140 (MQ=255)
                                                                                                     gCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACt                                 >  2:1037438/1‑140 (MQ=255)
                                                                                                     gCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACt                                 >  1:198852/1‑140 (MQ=255)
                                                                                                     gCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACt                                 >  1:841324/1‑140 (MQ=255)
                                                                                                     gCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACt                                 >  1:64657/1‑140 (MQ=255)
                                                                                                     gCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAAACGTCGTTGTCGATGCACAGGGGCGGCGTTGCTGCACt                                 >  1:445431/1‑140 (MQ=255)
                                                                                                     gCACAGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCATCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACt                                 >  1:202044/1‑140 (MQ=255)
                                                                                                       acaGAACTTAATGCAACACACCATCGCCCGGAGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAACTGGAAATTGTAGTTTACGATGAAATGGGGAGGCGTTGCTGCACTTGt                              >  2:588068/1‑140 (MQ=255)
                                                                                                         aGAACTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTc                             >  1:1273339/1‑140 (MQ=255)
                                                                                                           aaCTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCgg                           <  2:1375596/140‑1 (MQ=255)
                                                                                                           aaCTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCgg                           >  1:376162/1‑140 (MQ=255)
                                                                                                           aaCTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCCGCATCTCGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTATcag                           >  2:1387106/1‑138 (MQ=255)
                                                                                                           aaCTTAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGAACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCgg                           >  2:1279188/1‑140 (MQ=255)
                                                                                                               tAATGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTgg                       >  1:1025342/1‑140 (MQ=255)
                                                                                                                aaTGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTggg                      >  1:950093/1‑140 (MQ=255)
                                                                                                                aaTGCAACACACCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCAATGAACACGGGAGGAGTTGCTGCACTTGTCGGCTggg                      >  1:435496/1‑140 (MQ=255)
                                                                                                                      acacacCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGc                <  2:1328558/140‑1 (MQ=255)
                                                                                                                      acacacCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGc                >  1:239079/1‑140 (MQ=255)
                                                                                                                       cacacCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCa               >  1:1095362/1‑140 (MQ=255)
                                                                                                                       cacacCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCa               >  2:235890/1‑140 (MQ=255)
                                                                                                                       cacacCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGAGACCGCGTTGCTTCACTGGTCGGCTGGGTACGGCa               >  2:672008/1‑140 (MQ=255)
                                                                                                                          acCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGtt            >  1:1415877/1‑140 (MQ=255)
                                                                                                                          acCATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGtt            >  2:624499/1‑140 (MQ=255)
                                                                                                                           ccATCGCCCGGTGCCGGAGGGAAAGTTACGCGGCGTCGGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGttt           <  1:672008/140‑1 (MQ=255)
                                                                                                                            cATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGTTGAACAGGGGCGGCGTGGCTGCACTCGTCGACTGGGTACGGCAGtttt          >  2:681169/1‑140 (MQ=255)
                                                                                                                            cATCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGtttt          >  1:1255382/1‑140 (MQ=255)
                                                                                                                              tCGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTACACGTGACGGCTGAGTTCGGCAGTTTTgg        >  1:713018/1‑140 (MQ=255)
                                                                                                                               cGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTggg       <  2:239079/140‑1 (MQ=255)
                                                                                                                               cGCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTggg       <  2:1415877/140‑1 (MQ=255)
                                                                                                                                gCCCGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATAGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGa      >  1:1022165/1‑140 (MQ=255)
                                                                                                                                 cccGGTGCTTGAGGGAGAGGTATGCGGCGTGTGTCGGCCTCTGCATCTGGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGAt     <  2:715815/140‑1 (MQ=255)
                                                                                                                                 cccGGTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGAt     >  1:147465/1‑140 (MQ=255)
                                                                                                                                    ggTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGttt           >  1:382527/1‑131 (MQ=255)
                                                                                                                                    ggTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGttt           <  2:382527/131‑1 (MQ=255)
                                                                                                                                    ggTGCCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGaa  <  2:841324/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                       
GCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAA  >  NC_000913/629443‑629714

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