Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,974,081 T→A 7.5% K523M (AAG→ATG) ‡ flgK ← flagellar hook‑associated protein FlgK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,974,0810TA 7.5% 235.8 / 6.2 78K523M (AAG→ATG) ‡flgKflagellar hook‑associated protein FlgK
Reads supporting (aligned to +/- strand):  ref base T (38/34);  new base A (3/3);  total (41/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATCGAATCGCGAATCGACTGCTGGATCAGTGCGGT  >  NC_002947/4973936‑4974229
                                                                                                                                                 |                                                                                                                                                    
gTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACttgtt                                                                                                                                                  <  2:87161/150‑1 (MQ=255)
gTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACttgtt                                                                                                                                                  <  2:842905/150‑1 (MQ=255)
 tCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACttgttg                                                                                                                                                 >  1:752818/1‑150 (MQ=255)
 tCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACttgttg                                                                                                                                                 <  2:2111508/150‑1 (MQ=255)
   ggCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGTTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGTCATGGCGATGGACAACTTGTTGTc                                                                                                                                               <  2:827872/150‑1 (MQ=255)
   ggCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACAAGTTGTc                                                                                                                                               <  1:1193847/150‑1 (MQ=255)
      aTTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTg                                                                                                                                            >  2:1674499/1‑150 (MQ=255)
      aTTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTg                                                                                                                                            >  2:1412873/1‑150 (MQ=255)
      aTTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTg                                                                                                                                            <  2:660771/150‑1 (MQ=255)
         gAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACttg                                                                                                                                                    >  1:1422204/1‑139 (MQ=255)
         gAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACttg                                                                                                                                                    <  2:1422204/139‑1 (MQ=255)
           aCGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAg                                                                                                                                       <  2:390668/150‑1 (MQ=255)
            cGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAgg                                                                                                                                      <  2:1163620/150‑1 (MQ=255)
                   gTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTTTTGTCCTGCCCAGGCGCCACg                                                                                                                               >  2:78692/1‑150 (MQ=255)
                    tAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCATAACTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGc                                                                                                                              >  1:1870805/1‑150 (MQ=255)
                    tAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGc                                                                                                                              >  1:705993/1‑150 (MQ=255)
                     aaaGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCt                                                                                                                             >  1:1547156/1‑150 (MQ=255)
                        gCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGcc                                                                                                                          >  2:2065417/1‑150 (MQ=255)
                          tGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTg                                                                                                                            <  2:2379642/146‑1 (MQ=255)
                          tGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTg                                                                                                                            >  1:2379642/1‑146 (MQ=255)
                          tGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCgg                                                                                                                        >  1:988996/1‑150 (MQ=255)
                               ttGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCa                                                                                                                   >  2:2085438/1‑150 (MQ=255)
                               ttGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCa                                                                                                                   >  1:1670250/1‑150 (MQ=255)
                                    cGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTg                                                                                                              >  2:1449869/1‑150 (MQ=255)
                                     gCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTgg                                                                                                             <  2:1852661/150‑1 (MQ=255)
                                            tGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCt                                                                                                      <  2:988996/150‑1 (MQ=255)
                                              ccACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTg                                                                                                    <  1:1399819/150‑1 (MQ=255)
                                                aCCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGc                                                                                                  >  1:2338761/1‑150 (MQ=255)
                                                   gATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTc                                                                                               <  1:132414/150‑1 (MQ=255)
                                                   gATTGTGGTCTCGACGGCAAAGGAAAGCGGGATACCGCTACCATCAAGCCTAGGGTCACCGCCCGCACCCCCAATAGGCCCAGGAACGCACAACTGGCATTCACGCCCCGGGACCAGCCTGCCGGTACCAATTTGGGTACCCATAGCGTc                                                                                               >  2:644246/1‑150 (MQ=255)
                                                          gTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATAATGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGc                                                                                        <  1:262807/150‑1 (MQ=255)
                                                           tctcGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCt                                                                                       >  2:1739538/1‑150 (MQ=255)
                                                           tctcGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCt                                                                                       >  2:114621/1‑150 (MQ=255)
                                                             tcGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTg                                                                                     >  1:1796896/1‑150 (MQ=255)
                                                               gACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGc                                                                                                  <  2:1186653/135‑1 (MQ=255)
                                                               gACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGc                                                                                                  >  1:1186653/1‑135 (MQ=255)
                                                                aCGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAg                                                                                  <  2:1670250/150‑1 (MQ=255)
                                                                   gCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGccc                                                                               <  2:123437/150‑1 (MQ=255)
                                                                     aaaGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCCGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCt                                                                              >  1:1836416/1‑149 (MQ=255)
                                                                     aaaGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCt                                                                              <  2:1836416/149‑1 (MQ=255)
                                                                      aaGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTgg                                                                            <  1:1346728/150‑1 (MQ=255)
                                                                      aaGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTgg                                                                            <  2:2338761/150‑1 (MQ=255)
                                                                        ggTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGa                                                                          <  2:705993/150‑1 (MQ=255)
                                                                        ggTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGa                                                                          <  1:2085438/150‑1 (MQ=255)
                                                                                  ttACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTCGTCATGGGCAGGCACCACGCTGTTGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTg                                                                <  1:899635/150‑1 (MQ=255)
                                                                                     ccGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCa                                                             >  2:257161/1‑150 (MQ=255)
                                                                                      cGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAg                                                            >  1:1720481/1‑150 (MQ=255)
                                                                                        cTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGc                                                                                        >  1:1787599/1‑120 (MQ=255)
                                                                                        cTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGc                                                                                        <  2:1787599/120‑1 (MQ=255)
                                                                                        cTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTgg                                                                            >  2:639595/1‑132 (MQ=255)
                                                                                        cTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTgg                                                                            <  1:639595/132‑1 (MQ=255)
                                                                                         tACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGt                                                         <  1:1674499/150‑1 (MQ=255)
                                                                                           ccATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTAc                                                                                                         >  2:799529/1‑100 (MQ=255)
                                                                                           ccATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTAc                                                                                                         <  1:799529/100‑1 (MQ=255)
                                                                                           ccATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACAAGTTGTccatc                                                                                                                                           <  2:1727830/66‑4 (MQ=39)
                                                                                           ccATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACAAGTTGTccatc                                                                                                                                           >  1:1727830/1‑63 (MQ=38)
                                                                                             aTCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTAc                                                                                                                     <  1:495585/86‑1 (MQ=255)
                                                                                             aTCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTAc                                                                                                                     >  2:495585/1‑86 (MQ=255)
                                                                                              tCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAac                                                    >  1:724695/1‑150 (MQ=255)
                                                                                              tCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAac                                                    >  1:21251/1‑150 (MQ=255)
                                                                                               ccAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAAAAACCAGTTTCCCTGGGAGCCACCGCTGGCAGCGTCGAaca                                                   <  1:1171949/150‑1 (MQ=255)
                                                                                                cAGAATCGGTGCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAacac                                                  >  2:705940/1‑150 (MQ=255)
                                                                                                  gAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGCAcccca                                                <  2:2211172/150‑1 (MQ=255)
                                                                                                         gTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTAc                                         <  2:21251/150‑1 (MQ=255)
                                                                                                          tCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTAcc                                        <  2:1547156/150‑1 (MQ=255)
                                                                                                           cGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCCCGCTGCCGGCACTATAATGTGTACCCTTTTCGGCTAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCg                                       <  1:1497611/150‑1 (MQ=255)
                                                                                                               ggccgCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCAt                                   >  1:137543/1‑150 (MQ=255)
                                                                                                                   gCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATcgaa                               <  1:705940/150‑1 (MQ=255)
                                                                                                                     aTCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGTATcgaatc                             >  1:1278317/1‑150 (MQ=255)
                                                                                                                     aTCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATcgaatc                             >  1:1439589/1‑150 (MQ=255)
                                                                                                                     aTCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTGGTCGTCGAACAGCTTGTAGCCATGGGAGCCAACGCTGGCAGCGTCGAACACCACGCGTACCGGCATcgaatc                             >  1:1738194/1‑150 (MQ=255)
                                                                                                                          gCATCGGCACGGCGATGGACAACTTGTTGCCCTGCCCAGGCACCACGCTGCCGGTACCCATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATcgaatcgcgaa                        >  1:1123777/1‑150 (MQ=255)
                                                                                                                                   cGGCGATGGCCAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGGCGAACGCCACGCGTACCGGCATCGAATCGCGAATCGActgct               <  2:1796896/150‑1 (MQ=255)
                                                                                                                                     gCGATGGACAACTTGTTGTCCTGCCCAGGCGCCATGCTTCCGGTACTAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATCGAATCGCGAATCGACTGCTgg             <  2:1123777/150‑1 (MQ=255)
                                                                                                                                         aggACAACAAGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGa                                                                          >  2:982377/2‑85 (MQ=255)
                                                                                                                                         aggACAACAAGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGa                                                                          <  1:982377/84‑1 (MQ=255)
                                                                                                                                         aggACAACAAGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATCGAATCGCGAATCGACTGCTGGATCa         >  2:2062656/2‑150 (MQ=255)
                                                                                                                                            aCAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATCGAATCGCGAATCGACTGCTGGATCAGTg      >  2:1349681/1‑150 (MQ=255)
                                                                                                                                            aCAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATCGAATCGCGAATCGACTGCTGGATCAGTg      >  1:1127192/1‑150 (MQ=255)
                                                                                                                                              aaCTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCt                                                                              >  1:1668686/1‑76 (MQ=255)
                                                                                                                                              aaCTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCt                                                                              <  2:1668686/76‑1 (MQ=255)
                                                                                                                                               aCTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCTTACCGGCATCGAATCGCGAATCGACTGCTGGATCAGTGCgg   >  2:1374036/1‑150 (MQ=255)
                                                                                                                                               aCTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATCGAATCGCGAATCGACTGCTGGATCAGTGCgg   <  2:724695/150‑1 (MQ=255)
                                                                                                                                                cTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATCGAATCGCGAATCGACTGCTGGATCAGTGCGgt  >  2:29047/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATCGAATCGCGAATCGACTGCTGGATCAGTGCGGT  >  NC_002947/4973936‑4974229

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: