Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 551,295 C→T T286M (ACG→ATG)  ybcF → putative carbonate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913551,2950CT99.1% 418.2 / NA 110T286M (ACG→ATG) ybcFputative carbonate kinase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (30/79);  minor base A (1/0);  total (31/79)
Fisher's exact test for biased strand distribution p-value = 2.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 110 total aligned reads displayed.

CGTGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTT  >  NC_000913/551147‑551439
                                                                                                                                                    |                                                                                                                                                
 gTGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCt                                                                                                                                               <  1:788481/151‑1 (MQ=255)
   gCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTgg                                                                                                                                             <  1:1000657/151‑1 (MQ=255)
   gCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTg                                                                                                                                              <  1:1549238/150‑1 (MQ=255)
      aTTCGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTggcgg                                                                                                                                          <  2:504484/151‑1 (MQ=255)
        tCGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTggcggg                                                                                                                                         <  1:1118037/150‑1 (MQ=255)
         cGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTggcgggcg                                                                                                                                       <  2:682820/151‑1 (MQ=255)
         cGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTggcgggcg                                                                                                                                       <  2:1334202/151‑1 (MQ=255)
         cGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTggcgggc                                                                                                                                        <  1:953009/150‑1 (MQ=255)
            cATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAg                                                                                                                                    <  2:2033641/151‑1 (MQ=255)
             aTGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGc                                                                                                                                   <  1:1965218/151‑1 (MQ=255)
             aTGCCACACCGGAGGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGc                                                                                                                                   <  2:393011/151‑1 (MQ=255)
                ccACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGGGGGGGTATGTCAGAAGCCGTGGTAAACCCGCGGGGGTTGGGGGGGTATCGCGCATTGAAGAGATGCTtgcgggggggggggg                                                                                                                                >  2:397881/1‑142 (MQ=255)
                 cacaCCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGGGGGGGTATGTTAGAAACCGTGGGAAAACCGCGGGGGTTGGGGGGGTATCGCGAATTTAAGAGATGGTggcgggggggaggggg                                                                                                                               >  2:194563/1‑141 (MQ=255)
                 cacaCCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCgggg                                                                                                                               <  2:22266/151‑1 (MQ=255)
                  acacCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGGGGGGGTATGTCAGAAACCGTGGGAAAACCGCGGGGGTTGGGGGGGTATCGCGGATTGAAGAGATGCTggcggggggaagggggg                                                                                                                              >  1:834391/1‑140 (MQ=255)
                   cacCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGAc                                                                                                                             <  2:570469/151‑1 (MQ=255)
                     ccGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGGGGGGGTATGTTAGAAACCGGGGGTAACCCCCGGGGGTTGGGGCGGTATCGCGGATTGAAGAGATGTTtgcggggggaagggggg                                                                                                                              >  2:160619/1‑137 (MQ=255)
                      cGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTg                                                                                                                          <  2:593000/151‑1 (MQ=255)
                      cGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTg                                                                                                                          <  2:2132037/151‑1 (MQ=255)
                       ggATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGt                                                                                                                         <  1:538809/151‑1 (MQ=255)
                         aTGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTAt                                                                                                                       <  2:1173547/151‑1 (MQ=255)
                         aTGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTAt                                                                                                                       <  2:1777819/151‑1 (MQ=255)
                          tGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTAt                                                                                                                       <  1:1141277/150‑1 (MQ=255)
                           gAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTAttt                                                                                                                     <  1:446379/151‑1 (MQ=255)
                           gAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTAttt                                                                                                                     <  1:1944896/151‑1 (MQ=255)
                            aGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGGGTGGTTATGTCCGAAACCGGGGGAAAACCGCGTGGGTTGGGGGGTTATTGCGGATTGAAAAGAAGCTGGGGGGGGGAggggggggccggttttt                                                                                                                    >  2:638764/1‑140 (MQ=255)
                               tAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGc                                                                                                                  <  1:991670/150‑1 (MQ=255)
                                aGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAGCCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCt                                                                                                                 <  1:90785/150‑1 (MQ=255)
                                aGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCt                                                                                                                 <  2:1533095/150‑1 (MQ=255)
                                 gcgcCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGt                                                                                                               <  1:905275/151‑1 (MQ=255)
                                 gcgcCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGt                                                                                                               <  1:1990999/151‑1 (MQ=255)
                                   gcCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAg                                                                                                             <  2:2044271/151‑1 (MQ=255)
                                   gcCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAg                                                                                                             <  2:67024/151‑1 (MQ=255)
                                    ccATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGt                                                                                                            <  1:1504098/151‑1 (MQ=255)
                                      aTTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCg                                                                                                          <  2:142799/151‑1 (MQ=255)
                                         tGCCAAAGCCAATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGAGCTGTATTTCGCTGTAGTCGTa                                                                                                        <  1:801740/150‑1 (MQ=255)
                                           ccAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGc                                                                                                     <  2:1820292/151‑1 (MQ=255)
                                            cAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGGGGGTTTGTCAGAAGCCGTGGGAAAACCCCGGGGATTGGGGGGTTATCGCGAATTGAAGAGATGCTGGCGGGGGGAgggggggggcgtgttttt                                                                                                                   >  1:1355851/1‑124 (MQ=255)
                                               aGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTa                                                                                                 <  1:1649820/151‑1 (MQ=255)
                                                 ccGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGa                                                                                               <  1:18963/151‑1 (MQ=255)
                                                   gATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACa                                                                                             <  2:397421/151‑1 (MQ=255)
                                                      ggTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACAttt                                                                                          <  1:160043/151‑1 (MQ=255)
                                                      ggTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACAttt                                                                                          <  1:1287314/151‑1 (MQ=255)
                                                       gTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTtg                                                                                         <  1:1157899/151‑1 (MQ=255)
                                                         tCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTtgt                                                                                        <  1:1967362/150‑1 (MQ=255)
                                                         tCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTtgt                                                                                        <  1:85833/150‑1 (MQ=255)
                                                          cGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTtgtg                                                                                       <  2:656216/150‑1 (MQ=255)
                                                          cGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTtgtg                                                                                       <  2:1453778/150‑1 (MQ=255)
                                                            aTGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCt                                                                                    <  1:2136490/151‑1 (MQ=255)
                                                              ggggCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGa                                                                                  <  2:1714119/151‑1 (MQ=255)
                                                                  ccGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCg                                                                              <  1:1060942/151‑1 (MQ=255)
                                                                    gAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGAc                                                                            <  1:428245/151‑1 (MQ=255)
                                                                          aaCGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGa                                                                      <  1:1684759/151‑1 (MQ=255)
                                                                            cggcggTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGAcg                                                                    <  2:1781736/151‑1 (MQ=255)
                                                                                  tGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGGGGAGCGGGGGACTGTTTTTTGCTGTAGGCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGGCtt                                                              >  2:1157349/1‑151 (MQ=255)
                                                                                  tGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCtt                                                              >  1:1803097/1‑151 (MQ=255)
                                                                                   gAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGGCCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTTTGCCTGATGCGACGCTTGACGCGTCtt                                                              >  2:1745790/1‑150 (MQ=255)
                                                                                   gAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTa                                                             <  1:1242258/151‑1 (MQ=255)
                                                                                   gAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTCTAGTCGCAGGCATTAGACATTTTTGCCTGATGCGACGCTTGACGCGTCtt                                                              >  1:2078800/1‑150 (MQ=255)
                                                                                    aGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTa                                                             <  2:931569/150‑1 (MQ=255)
                                                                                      tGGTTATGTCAGAAACCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTTGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATc                                                           >  1:1425319/1‑150 (MQ=255)
                                                                                         ttATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGGAGCGGGGGACTGTATTTCGCTGTAGTCGGAGGCATTAGAAATTTGTGCCTGATGCGACGCTTGACGCGTCTTTTCAgg                                                        >  2:2116217/1‑150 (MQ=255)
                                                                                          tATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGcc                                                      >  1:370533/1‑151 (MQ=255)
                                                                                           aTGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGcc                                                      >  1:895434/1‑150 (MQ=255)
                                                                                             gTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTa                                                    <  2:1631241/150‑1 (MQ=255)
                                                                                             gTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTa                                                    >  1:1902054/1‑150 (MQ=255)
                                                                                             gTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTAc                                                   <  1:380748/151‑1 (MQ=255)
                                                                                              tCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCTTGATGCGACGCTTGACGCGTCTTATCAGGCCTAc                                                   >  1:897468/1‑150 (MQ=255)
                                                                                              tCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTAc                                                   >  1:1858573/1‑150 (MQ=255)
                                                                                               cAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACaa                                                 <  1:746628/151‑1 (MQ=255)
                                                                                               cAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGGATTTTTCTGTAGTCGTAGGCATTTGACATTTTTGCCTGAGGCGACGCTTGACGCGTCTTCTCCGGCCt                                                     >  2:356027/1‑147 (MQ=255)
                                                                                                  aaGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCg                                              <  1:1928466/151‑1 (MQ=255)
                                                                                                       gTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGcc                                         <  2:1503863/151‑1 (MQ=255)
                                                                                                       gTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGcc                                         <  1:131070/151‑1 (MQ=255)
                                                                                                       gTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGGGcc                                         >  1:561280/1‑151 (MQ=255)
                                                                                                            aaaCCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCCGTATTTACCTGTAGTCGGAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAAGCCTACAACCCGTGCCGCATc                                    >  1:159612/1‑151 (MQ=39)
                                                                                                             aaCCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATcc                                   <  1:575022/151‑1 (MQ=255)
                                                                                                                ccGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGGCCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGGCTGATGCGACGCTTGACGCGGCTTATCAGGCCTACAACCGGGGCCGCATCCgg                                 >  2:1357553/1‑150 (MQ=255)
                                                                                                                  gcgTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGGGGGGGCCTGTATTTCGCCGTAAGCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCa                               >  1:1252883/1‑150 (MQ=38)
                                                                                                                  gcgTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCaa                              <  1:571216/151‑1 (MQ=255)
                                                                                                                  gcgTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCaa                              <  1:935674/151‑1 (MQ=255)
                                                                                                                   cgTGGATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAAt                             <  1:1651420/151‑1 (MQ=255)
                                                                                                                      ggATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTgg                          <  1:990648/151‑1 (MQ=255)
                                                                                                                       gATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGGCCTGTATTTTGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAACGGTGCCGCATCCGGCAAtttg                          >  2:1763983/1‑148 (MQ=255)
                                                                                                                       gATTGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGt                         <  1:201355/151‑1 (MQ=255)
                                                                                                                          tGGGGCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGc                       <  1:821784/150‑1 (MQ=255)
                                                                                                                           ggggCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGcac                     <  1:313870/151‑1 (MQ=255)
                                                                                                                             ggCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACaa                   <  1:651714/151‑1 (MQ=255)
                                                                                                                              gCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAAt                  <  1:1995114/151‑1 (MQ=255)
                                                                                                                              gCGTTATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAAt                  <  1:474120/151‑1 (MQ=255)
                                                                                                                                  tATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGcc               <  2:1899451/150‑1 (MQ=255)
                                                                                                                                  tATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGcc               <  2:1022465/150‑1 (MQ=255)
                                                                                                                                  tATCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGCATTTTTCTGTAGTCGTAGGCATTAGACAATTTTGCCTGATGCGGCGCTTTACGCGTCTTATAAAGCCCACAAACAGTGCCGCATCCCGCAATTGGGGCACAAAGcc               >  2:73163/1‑150 (MQ=255)
                                                                                                                                   aTCGCGAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTTGCTGTAGTCGTAGGCATTAGACATTTGTGGCTGATGCGACGCTTGACGCGTCTTATTAGGCCTACAACCGGGGGCGCATCCGGCAAATTGTGGACAATGcc               >  2:1855454/1‑149 (MQ=255)
                                                                                                                                     cgcgAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTgat           <  1:416139/151‑1 (MQ=255)
                                                                                                                                     cgcgAATTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTga            <  1:1225838/150‑1 (MQ=255)
                                                                                                                                         aaTTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTTTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACCATGCCTgattcg        >  1:1400059/1‑150 (MQ=255)
                                                                                                                                          aTTGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTgatgcgat      <  2:1778133/151‑1 (MQ=255)
                                                                                                                                            tGAAGAGATGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTTTCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTgctgcgatg     >  2:2004872/1‑150 (MQ=255)
                                                                                                                                                gagaTGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATTCCTGATGCGATGCtt  >  2:898406/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CGTGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTT  >  NC_000913/551147‑551439

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: