Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,493,494 A→T 5.7% W462R (TGG→AGG)  PP_1306 ← Pyocin S‑type killer domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,493,4940AT 5.7% 313.5 / 4.6 103W462R (TGG→AGG) PP_1306Pyocin S‑type killer domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (42/55);  new base T (4/2);  total (46/57)
Fisher's exact test for biased strand distribution p-value = 4.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 120 total aligned reads displayed.

TTCTCAACCCCCGCCAGAAGGTAGGGCAGGTCGATGCTGGTCGAGGATGCCGCCAAGGCAGCAAGGTCGGGCCCACCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGCCAAGGCAGAGTAAGCCGCTTCGGCTACGGCAAAAGTGGTAGCCCCAATTGGGATTACGGCAGGTGT  >  NC_002947/1493346‑1493641
                                                                                                                                                    |                                                                                                                                                   
ttCTCAACCCCCGCCAGAAGGTAGGGCAGGTCGATGCTGGTCGAGGATGCCGCCAAGGCAGCAAGGTCGGGCCCACCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATggccag                                                                                                                                                    <  1:895118/150‑1 (MQ=255)
ttCTCAACCCCCGCCAGAAGGTAGGGCAGGTCGATGCTGGTCGAGGATGCCGCCAAGGCAGCAAGGTCGGGCCCACCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATggccag                                                                                                                                                    >  1:1200621/1‑150 (MQ=255)
ttCTCAACCCCCGCCAGAAGGTAGGGCAGGTCGATGCTGGTCGAGGATGCCGCCAAGGCAGCAAGGTCGGGCCCACCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATggccag                                                                                                                                                    >  2:186591/1‑150 (MQ=255)
  ctcAACCCCCGCCAGAAGGTAGGGCAGGTCGATGCTGGTCGAGGATGCCGCCAAGGCAGCAAGGTCGGGCCCACCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATggccaggc                                                                                                                                                  >  2:1428412/1‑150 (MQ=255)
    cAACCCCCGCCAGAAGGTAGGGCAGGTCGATGCTGGTCGAGGATGCCGCCAAGGCAGCAAGGTCGGGCCCACCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATggccaggcca                                                                                                                                                <  2:295317/150‑1 (MQ=255)
         cccGCCAGAAGGTAGGGCAGGTCGATGCTGGTCGAGGATGCCGCCAAGGCAGCAAGGTCGGGCCCACCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCTGGCCATcaac                                                                                                                                           >  1:1169801/1‑147 (MQ=255)
         cccGCCAGAAGGTAGGGCAGGTCGATGCTGGTCGAGGATGCCGCCAAGGCAGCAAGGTCGGGCCCACCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCa                                                                                                                                           >  1:267138/1‑150 (MQ=255)
         cccGCCAGAAGGTAGGGCAGGTCGATGCTGGTCGAGGATGCCGCCAAGGCAGCAAGGTCGGGCCCACCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCa                                                                                                                                           >  2:1917003/1‑150 (MQ=255)
         cccGCCAGAAGGTAGGGCAGGTCGATGCTGGTCGAGGATGCCGCCAAGGCAGCAAGGTCGGGCCCACCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCa                                                                                                                                           >  1:1525210/1‑150 (MQ=255)
           cGCCAGAAGGTAGGGCAGGTCGATGCTGGTCGAGGATGCCGCCAAGGCAGCAAGGTCGGGCCCACCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAg                                                                                                                                         <  2:457256/150‑1 (MQ=255)
           atctaGAAGGTAGGGCAGGTCGATGCTGGTCGAGGATGAAACCAAGGCAGCAAGGTCGGGCCCACCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAg                                                                                                                                         <  1:625135/146‑1 (MQ=255)
            gCCAGAAGGTAGGGCAGGTCGATGCTGGTCGAGGATGCCGCCAAGGCAGCAAGGTCGGGCCCACCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAgg                                                                                                                                        <  2:2430097/150‑1 (MQ=255)
             ccAGAAGGTAGGGCAGGTCGATGCTGGTCGAGGATGCCGCCAAGGCAGCAAGGTCGGGACCACCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTCGGTGCCCAAGCCATCGCAAggg                                                                                                                                       >  2:204419/1‑150 (MQ=255)
                gAAGGTAGGGCAGGTCGATGCTGGTCGAGGATGCCGCCAAGGCAGCAAGGTCGGGCCCACCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACACGACATCACACTGGCCCAACGATGATGGCCAGGTCTTCGCAAggcgat                                                                                                                                    >  2:625135/1‑146 (MQ=255)
                gAAGGTAGGGCAGGTCGATGCTGGTCGAGGATGCCGCCAAGGCAGAAAAGTCGGTGCCGTCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGc                                                                                                                                    <  2:2276501/150‑1 (MQ=255)
                   ggTAGGGCAGGTCGATGCTGGTCGAGGATGCCGCCAAGGCAGCAAGGTCGGGCCCACCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTa                                                                                                                                 >  1:394195/1‑150 (MQ=255)
                       gggCAGGTCGATGCTGGTCGAGGATGCCGCCAAGGCAGCAAGGTCGGGCCCACCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGc                                                                                                                             <  1:381227/150‑1 (MQ=255)
                          cAGGTCGATGCTGGTCGAGGATGCCGCCAAGGCAGCAAGGTCGGGCCCACCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGAc                                                                                                                          <  2:1525210/150‑1 (MQ=255)
                             gTCGATGCTGGTCGAGGATGCCGCCAAGGCAGCAAGGTCGGGCCCACCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCTGGCCATcaacg                                                                                                                                          <  1:137348/131‑5 (MQ=255)
                             gTCGATGCTGGTCGAGGATGCCGCCAAGGCAGCAAGGTCGGGCCCACCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCTGGCCATcaacg                                                                                                                                          >  2:137348/1‑127 (MQ=255)
                             gTCGATGCTGGTCGAGGATGCCGCCAAGGCAGCAAGGTCGGGCCCACCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGa                                                                                                                       >  2:278717/1‑150 (MQ=255)
                               cGATGCTGGTCGAGGATGCCGCCAAGGCAGCAAGGTCGGGCCCACCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAg                                                                                                                     <  1:981746/150‑1 (MQ=255)
                                gATGCTGGTCGAGGATGCCGCCAAGGCAGCAAGGTCGGGCCCACCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAgg                                                                                                                    <  2:394195/150‑1 (MQ=255)
                                gATGCTGGTCGAGGATGCCGCCAAGGCAGCAAGGTCGGGCCCACCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAgg                                                                                                                    >  1:1870894/1‑150 (MQ=255)
                                  tGCTGGTCGAGGATGCCGCCAAGGCAGCAAGGTCGGGCCCACCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTGCATGTTCAGCTACGTTTATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGa                                                                                                                  >  2:63702/1‑150 (MQ=255)
                                     tGGTCGAGGATGCCGCCAAGGCAGCAAGGTCGGGCCCACCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAg                                                                                                               <  2:1119174/150‑1 (MQ=255)
                                       gTCGAGGATGCCGCCAAGGCAGCAAGGTCGGGCCCACCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGAc                                                                                                                          <  1:2229406/137‑1 (MQ=255)
                                       gTCGAGGATGCCGCCAAGGCAGCAAGGTCGGGCCCACCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGAc                                                                                                                          >  2:2229406/1‑137 (MQ=255)
                                           aGGATGCCGCCAAGGCAGCAAGGTCGGGCCCACCGGCAGGGCTAAGAAAACGGTGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAg                                                                                                         <  1:145805/150‑1 (MQ=255)
                                               tGCCGCCAAGGCAGCAAGGTCGGGCCCACCGGCAGTGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGcaac                                                                                                     <  1:295501/150‑1 (MQ=255)
                                               tGCCGCCAAGGCAGCAAGGTCGGGCCCACCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGATCAACACGACCGAAGGGACAGTCGTAGcaac                                                                                                     <  1:663530/150‑1 (MQ=255)
                                                  cgccAAGGCAGCAAGGTCGGGCCCACCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAg                                                                                                  <  1:1428412/150‑1 (MQ=255)
                                                  cgccAAGGCAGCAAGGTCGGGCCCACCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAg                                                                                                  >  1:1806807/1‑150 (MQ=255)
                                                   gccAAGGCAGCAAGGTCGGGCCCACCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATc                                                                                                                                              >  2:921808/1‑105 (MQ=255)
                                                   gccAAGGCAGCAAGGTCGGGCCCACCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATc                                                                                                                                              <  1:921808/105‑1 (MQ=255)
                                                   accAAGCCAACAACGTCGGGCTCACCCCCAGTGCTAAAGCAGTCGAAGAGGGTGATCAGGAGGTACAGGAGTTCACTGTTACCCAGCGCCAATGGCCAGGCCATCGCAAGCGTCCCTACTGCCACCCAAGGGACAGTCGCCGcaacaaac                                                                                                 >  2:1031685/2‑148 (MQ=255)
                                                    ccAAGGCAGCAAGGTCGGGCCCACCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGcc                                                                                                >  1:208348/1‑150 (MQ=255)
                                                         gcagcaAGGTCGGGCCCACCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTa                                                                                           >  2:2316130/1‑150 (MQ=255)
                                                           agcaAGGTCGGGCCCACCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAg                                                                                         >  2:2363539/1‑150 (MQ=255)
                                                            gaaaGGTAGGGACCACCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGt                                                                                        >  2:699618/3‑150 (MQ=255)
                                                                  tCGGGGCCGCCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGATCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCt                                                                                  <  2:134738/150‑1 (MQ=255)
                                                                   cGGGCCCACCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTg                                                                                 <  2:60122/150‑1 (MQ=255)
                                                                    gggCCCACCGGCAGGGCTAAGGCTGTCGAGTTGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGc                                                                                >  1:1025572/1‑150 (MQ=255)
                                                                    gggCCCACCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGc                                                                                <  2:2147394/150‑1 (MQ=255)
                                                                    gggCCCACCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGc                                                                                <  1:173919/150‑1 (MQ=255)
                                                                        ccACCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGATAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGcgatc                                                                            <  2:619498/150‑1 (MQ=255)
                                                                        ccACCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGcgatc                                                                            >  2:2197070/1‑150 (MQ=255)
                                                                          aCCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGcgatcga                                                                          <  1:1347782/150‑1 (MQ=255)
                                                                          aCCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGcgatcga                                                                          <  2:1870894/150‑1 (MQ=255)
                                                                                   ttccgatcTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGAATGGGTGCCTACTGCGACCTAATATTCCTGCATAACAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGCCaa                                                                 <  1:414520/143‑1 (MQ=255)
                                                                                   gCTAAGGCTGTCGATTTGGGTGCTGATCAGGTACTGGCTTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGCCaa                                                                 <  1:2374815/150‑1 (MQ=255)
                                                                                    cTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGCCAAg                                                                >  1:1876813/1‑150 (MQ=255)
                                                                                    cTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGCCAAg                                                                <  1:1895153/150‑1 (MQ=255)
                                                                                       aGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGa                                                                                                                  <  2:1975095/97‑1 (MQ=255)
                                                                                       aGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGa                                                                                                                  >  1:1975095/1‑97 (MQ=255)
                                                                                       aGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGATCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGCCAAGGCa                                                             <  1:524080/150‑1 (MQ=255)
                                                                                          cTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAAGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGCCAAGGCagag                                                          <  1:1217730/150‑1 (MQ=255)
                                                                                              cGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGCCAAGGCAGAGTAAg                                                      <  1:1917003/150‑1 (MQ=255)
                                                                                                     ggTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGCCAAGGCAGAGTAAGCCGCTTc                                               <  1:2197070/150‑1 (MQ=255)
                                                                                                     ggTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGCCAAGGCAGAGTAAGCCGCTTc                                               <  2:1656601/150‑1 (MQ=255)
                                                                                                       tGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGCCaa                                                                 >  1:1030667/1‑130 (MQ=255)
                                                                                                          tgatgaGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCTGGTAAGTGCCGCTGCGATCGATTCAGCCAAGGCAGAGTAAGCCGCTTCGGCTa                                          <  1:278717/150‑1 (MQ=255)
                                                                                                           gatgaGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGCCAAGGCAGAGTAAGCCGCTTCGGCTAc                                         <  2:1025572/150‑1 (MQ=255)
                                                                                                                     tGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGCCAAGGCAGAGTAAGCCGCTTCGGCTACGGCAAAAGTg                               >  2:349778/1‑150 (MQ=255)
                                                                                                                       gCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGCCAAGGCAGAGTAAGCCGCTTCGGCTACGGCAAAAGTGGt                             <  2:1353465/150‑1 (MQ=255)
                                                                                                                          ttCACTGTTGCGCAGCGTTGATGGCTAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGCCAAGGCAGAGTAAGCCGCTTCGGCTACGGCAAAAGTGGTAGc                          <  2:1136443/150‑1 (MQ=255)
                                                                                                                          ttCACTGTTGCCCAGCGTTGATGGGCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGCCAAGGCAGAGTAAGCCGCTTCGGCTACGGCAAAAGTGGTAGc                          <  2:1136311/150‑1 (MQ=255)
                                                                                                                          ttCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGCCAAGGCAGAGTAAGCCGCTTCGGCTACGGCAAAAGTGGt                             >  1:491813/1‑147 (MQ=255)
                                                                                                                          ttCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGCCAAGGCAGAGTAAGCCGCTTCGGCTACGGCAAAAGTGGt                             <  2:491813/147‑1 (MQ=255)
                                                                                                                          ttCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGCCAAGGCAGAGTAAGCCGCTTCGGCTACGGCAAAAGTGGTAGc                          >  1:175993/1‑150 (MQ=255)
                                                                                                                            cACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGCCAAGGCAGAGTAAGCCGCTTCGGCTACGGCaaaa                                  >  1:1702793/1‑140 (MQ=255)
                                                                                                                            cACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGCCAAGGCAGAGTAAGCCGCTTCGGCTACGGCaaaa                                  <  2:1702793/140‑1 (MQ=255)
                                                                                                                               tGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGCCAAGGCAGAGTAAGc                                                     >  2:564669/1‑118 (MQ=255)
                                                                                                                               tGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGCCAAGGCAGAGTAAGc                                                     <  1:564669/118‑1 (MQ=255)
                                                                                                                               tGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGGCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGCCAAGGCAGAGTAAGCCGCTTCGGCTCCGGCAAAAGTGGTAGCCCCaa                     >  2:921485/1‑150 (MQ=255)
                                                                                                                                 ttGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGc                                                                    >  1:2385749/1‑101 (MQ=255)
                                                                                                                                 ttGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGc                                                                    <  2:2385749/101‑1 (MQ=255)
                                                                                                                                 ttGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGCCAAGGCAGAGTAAGCCGCTTCGGCTACGGCAAAAGTGGTAGCCCCAAtt                   <  1:283624/150‑1 (MQ=255)
                                                                                                                                    cccAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGCCAAGGCAGAGTAAGCCGCTTCgg                                             <  1:946038/121‑1 (MQ=255)
                                                                                                                                    cccAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGCCAAGGCAGAGTAAGCCGCTTCgg                                             >  2:946038/1‑121 (MQ=255)
                                                                                                                                    cccAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGCCAAGGCAGAGTAAGCCGCTTCGGCTACGGCAAAAGTGGTAGCCCCAATTggg                <  2:1806807/150‑1 (MQ=255)
                                                                                                                                      cAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGCCAAGGCAGAGTAAGCCGCTTCGGCTACGGCAAAAGTGGTAGCCCCAATTGGGAt              >  2:433787/1‑150 (MQ=255)
                                                                                                                                      cAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGCCAAGGCAGAGTAAGCCGCTTCGGCTACGGCAAAAGTGGTAGCCCCAATTGGGAt              >  2:159938/1‑150 (MQ=255)
                                                                                                                                      cAGCGTTGATGGCCAGGACGTGTCGAGGGGGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGCCAAGGCAGAGTAAGCCGCTTCGGCTACGGCAAAAGTGGTAGCCCCAATTGGGAt              >  2:1402156/1‑150 (MQ=255)
                                                                                                                                       aGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGCCAAGGCAGAGTAAGCCGCTTCGGCTACGGCAAAAGTGGTAGc                          >  1:769828/1‑137 (MQ=255)
                                                                                                                                       aGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGCCAAGGCAGAGTAAGCCGCTTCGGCTACGGCAAAAGTGGTAGc                          <  2:769828/137‑1 (MQ=255)
                                                                                                                                         ttgcgATGGCCTGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGCCAAGGCAGAGTAAGCCGCTTCGGCTACGGCAAAAGTGGTAGCCCCAATTGGGATTAc           >  2:2389771/5‑150 (MQ=255)
                                                                                                                                         ttgcgATGGCCTGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGCCAAGGCAGAGTAAGCCGCTTCGGCTACGGCAAAAGTGGTAGCCCCAATTGGGATTAc           >  2:1159295/5‑150 (MQ=255)
                                                                                                                                         ctttGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGCCAAGGCAGAGTAAGCCGCTTCGGCTACGGCAAAAGTGGTAGCCCCAATTGGGATTAc           <  1:346191/148‑1 (MQ=255)
                                                                                                                                         cGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGCCAAGGCAGAGTAAGCCGCTTCGGCTACGGCAAAAGTGGTAGCCCCAATTGGGATTAc           <  1:389579/150‑1 (MQ=255)
                                                                                                                                              aTGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGCCAAGGCAGAGTaa                                                       <  1:312038/101‑1 (MQ=255)
                                                                                                                                                ggccaggccaTCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGCCAAGGCAGAGTAAGCCGCTTCGGCTACGGCAACAGTGCTAGCCCCAATTGGGATTACGGCAGgt    <  1:63702/150‑1 (MQ=255)
                                                                                                                                                ggccaggccaTCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGCCAAGGCAGAGTAAGCCGCTTCGGCTACGGCAAAAGTGGTAGCCCCAATTGGGATTACGGCAGgt    >  1:840533/1‑150 (MQ=255)
                                                                                                                                                 gccaggccaTCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGCCAAGGCAGAGTAAGCCGCTTCGGCTACGGCAAAAGTGGTAGCCCCAATTGGGATTACGGCAGgtg   >  1:1807929/1‑150 (MQ=255)
                                                                                                                                                  ccaggccaTCGCAAGGGTGCCTACTGCGACCGAAg                                                                                                                     >  2:252857/1‑35 (MQ=255)
                                                                                                                                                  ccaggccaTCGCAAGGGTGCCTACTGCGACCGAAg                                                                                                                     >  2:251519/1‑35 (MQ=255)
                                                                                                                                                  ccaggccaTCGCAAGGGTGCCTACTGCGACCGAAg                                                                                                                     <  1:252857/35‑1 (MQ=255)
                                                                                                                                                  ccaggccaTCGCAAGGGTGCCTACTGCGACCGAAg                                                                                                                     <  1:251519/35‑1 (MQ=255)
                                                                                                                                                  ccaggccaTCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGCCAAGGCAGAGTAAGCCGCTTCGGCTACGGCAAAAGTGGTAGCCCCAATTGGGATTACGGCAGgtgt  >  2:2126978/1‑150 (MQ=255)
                                                                                                                                                  ccaggccaTCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGCCAAGGCAGAGTAAGCCGCTTCGGCTACGGCAAAAGTGGTAGCCCCAATTGGGATTACGGCAGgtgt  >  2:2125869/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TTCTCAACCCCCGCCAGAAGGTAGGGCAGGTCGATGCTGGTCGAGGATGCCGCCAAGGCAGCAAGGTCGGGCCCACCGGCAGGGCTAAGGCTGTCGAGTGGGGTGCTGATGAGGTACTGGCGTTCACTGTTGCCCAGCGTTGATGGCCAGGCCATCGCAAGGGTGCCTACTGCGACCGAAGGGACAGTCGTAGCAACAAGCCGGGTAAGTGCCGCTGCGATCGATTCAGCCAAGGCAGAGTAAGCCGCTTCGGCTACGGCAAAAGTGGTAGCCCCAATTGGGATTACGGCAGGTGT  >  NC_002947/1493346‑1493641

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: