Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,484,207 C→A 5.8% T165N (ACC→AAC)  PP_3976 → isochorismatase family hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,484,2070CA 5.8% 251.8 / 2.3 80T165N (ACC→AAC) PP_3976isochorismatase family hydrolase
Reads supporting (aligned to +/- strand):  ref base C (38/37);  new base A (2/3);  total (40/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

AGAGGTCTGCAACACCACCTTCGCGCTGTCCGCCATGCACGATGCAGGCTATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCAAAGAGCATTCGGGCGCCTACGGCATGGGCGTCGACTACGCCTATACCATGGTGCACAAAGCGCCAGAGC  >  NC_002947/4484058‑4484353
                                                                                                                                                     |                                                                                                                                                  
agagGTCTGCAACACCACCTTCGCGCTGTCCGCCATGCACGATGCAGGCTATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTAc                                                                                                                                                    <  1:712132/150‑1 (MQ=255)
agagGTCTGCAACACCACCTTCGCGCTGTCCGCCATGCACGATGCAGGCTATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTAc                                                                                                                                                    >  1:269106/1‑150 (MQ=255)
       tGCAACACCACCTTCGCGCTGTCCGCCATGCACGATGCAGGCTATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGgcag                                                                                                                                             >  2:191652/1‑150 (MQ=255)
        gCAACACCACCTTCGCGCTGTCCGCCATGCACGATGCAGGCTATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGgcagc                                                                                                                                            >  2:238965/1‑150 (MQ=255)
         cAACACCACCTTCGCGCTGTCCGCCATGCACGATGCAGGCTATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGCTGCCAGGTAACTGgcacca                                                                                                                                           <  1:901586/150‑1 (MQ=255)
          aaCACCACCTTCGCGCTGTCCGCCATGCACGATGCAGGCTATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGcagcag                                                                                                                                          <  2:1509610/150‑1 (MQ=255)
              ccaccTTCGCGCTGTCCGCCATGCACGATGCAGGCTATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGc                                                                                                                                      >  1:707515/1‑150 (MQ=255)
                  cTTCGCGCTGTCCGCCATGCACGATGCAGGCTATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACt                                                                                                                                  <  2:269106/150‑1 (MQ=255)
                     cgcgCTGTCCGCCATGCACGATGCAGGCTATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGa                                                                                                                               <  2:707515/150‑1 (MQ=255)
                       cgcTGTCCGCCATGCACGATGCAGGCTATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGt                                                                                                                             >  2:337084/1‑150 (MQ=255)
                       cgcTGTCCGCCATGCACGATGCAGGCTATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGt                                                                                                                             >  2:139705/1‑150 (MQ=255)
                        gcTGTCCGCCATGCACGATGCAGGCTATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTg                                                                                                                            >  1:1769074/1‑150 (MQ=255)
                            tCCGCCATGCACGATGCAGGCTATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAg                                                                                                                        >  2:1925604/1‑150 (MQ=255)
                                ccATGCACGATGCAGGCTATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGcgcg                                                                                                                    >  1:248356/1‑150 (MQ=255)
                                     cACGATGCAGGCTATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGt                                                                                                                             >  1:578167/1‑136 (MQ=255)
                                     cACGATGCAGGCTATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGt                                                                                                                             <  2:578167/136‑1 (MQ=255)
                                      aCGATGCAGGCTATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTg                                                                                                                                       <  2:943376/125‑1 (MQ=255)
                                      aCGATGCAGGCTATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTg                                                                                                                                       >  1:943376/1‑125 (MQ=255)
                                          tGCAGGCTATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGaaaaa                                                                                                          <  2:1769074/150‑1 (MQ=255)
                                            cAGGCTATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCGGCAGGTGATACTGGAGTGGCAGAGCGACTGGAAAAAcc                                                                                                        >  2:371289/1‑150 (MQ=255)
                                             aGGCTATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACcg                                                                                                       >  1:1767364/1‑150 (MQ=255)
                                             aGGCTATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACcg                                                                                                       >  1:232389/1‑150 (MQ=255)
                                             aGGCTATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACcg                                                                                                       <  2:232389/150‑1 (MQ=255)
                                               gCTATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCAAGACCAAGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGGGCTACTGGAGTGGCAGCGCGACTGGAAAAACcgcg                                                                                                     >  1:1711859/1‑150 (MQ=255)
                                                       aTCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTAc                                                                                             <  1:618029/150‑1 (MQ=255)
                                                           atatGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATg                                                                                         >  2:1242604/1‑150 (MQ=255)
                                                           atatGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATg                                                                                         >  2:1754933/1‑150 (MQ=255)
                                                             atGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTCACGACTCAGCTATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCt                                                                                       <  1:1573427/150‑1 (MQ=255)
                                                              tGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTg                                                                                      >  1:183355/1‑150 (MQ=255)
                                                                      gACGCTTCCGGGGGCACGACCAAGGAGGCTCACGAATACGCCATGCAGCGCATGATCCAGGCGGGAGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGAGCGAATCGAAAAAACGCGACACCTACGATGCTGTCATGgcg                                                                              >  2:1852473/1‑150 (MQ=255)
                                                                          cTTCCGGGGGCACGACCAAGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGTAAAAACCGCGACACCTACGATGCTGTCATGGCGCTgg                                                                          <  1:1568948/150‑1 (MQ=255)
                                                                          cTTCCGGGGGCACGACCAAGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTgg                                                                          <  1:1925604/150‑1 (MQ=255)
                                                                           ttCCGGGGGCACGACCAAGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTg                                                                           <  2:1125744/148‑1 (MQ=255)
                                                                           ttCCGGGGGCACGACCAAGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTg                                                                           >  1:1125744/1‑148 (MQ=255)
                                                                           ttCCGGGGGCACGACCAAGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGt                                                                         >  1:546353/1‑150 (MQ=255)
                                                                           ttCCGGGGGCACGACAAAAGAGGATCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACCGGCGTGCGAGCCCTACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGt                                                                         <  2:759542/150‑1 (MQ=255)
                                                                             ccGGGGGCACGACCAAGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCa                                                                       >  2:1353474/1‑150 (MQ=255)
                                                                               gggggCACGACCAAGGAGGCTAAAGACAACGCGAAGGAGCGCAGTAGGCAGGCGGGGGTGGGGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCaaa                                                                     <  2:211964/150‑1 (MQ=255)
                                                                                 gggCACGACCAAGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCAAaga                                                                   <  1:1310503/150‑1 (MQ=255)
                                                                                   gCAAGACCAAGGAGGCTCACGACTACGCCATGCAGCGCATTATCCAGGCGGGCGTGGGGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCAAAGAGc                                                                 <  2:459080/150‑1 (MQ=255)
                                                                                     aCGACCAAGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCAAAGAGCAt                                                               >  1:552381/1‑150 (MQ=255)
                                                                                     aCGACCAAGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTCCCAGTTACCTGGCAGCAGGTGCGACTGGATTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCAAAGAGCAt                                                               >  1:1691119/1‑150 (MQ=255)
                                                                                       gACCAAGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTg                                                                                                                                       >  1:15075/1‑76 (MQ=255)
                                                                                       gACCAAGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTg                                                                                                                                       <  2:15075/76‑1 (MQ=255)
                                                                                       gACCAAGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGAGTGGTAAAACCGCGGCACATATGATGCTGTCATGGCGCTGGTAAAAGAGCATTc                                                             >  2:543985/1‑150 (MQ=255)
                                                                                       gACCAAGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCAAAGAGCATTc                                                             >  1:696480/1‑150 (MQ=255)
                                                                                       gACCAAGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGCGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCCACACCTACGATGCTGTCATGGCGCTGGTCAAAGAGCATTc                                                             <  1:1852473/150‑1 (MQ=255)
                                                                                           aaGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCAAAGAGCATTCGGgc                                                         <  2:1570660/150‑1 (MQ=255)
                                                                                           aaGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCAAAGAGCATTCGGgc                                                         <  2:1767364/150‑1 (MQ=255)
                                                                                                 gCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCAAAGAGCATTCGGGCGCCTAc                                                   <  1:371289/150‑1 (MQ=255)
                                                                                                   tCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCt                                                                                       >  2:1352423/1‑112 (MQ=255)
                                                                                                   tCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCt                                                                                       <  1:1352423/112‑1 (MQ=255)
                                                                                                   tCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCAAAGAGCATTCGGGCGCCTACgg                                                 >  2:439758/1‑150 (MQ=255)
                                                                                                    cACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCAAAGAGCATTCGGGCGCCTACGGc                                                <  1:1353474/150‑1 (MQ=255)
                                                                                                        aCTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCAAAGAGCATTCGGGCGCCTACGGCATgg                                            <  1:1242604/150‑1 (MQ=255)
                                                                                                          tACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCAAAGAGCATTCGGGCGCCTACGGCATGGGc                                          >  2:739450/1‑150 (MQ=255)
                                                                                                               cATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGTCGCCTACGATGCTGTCATGGCGCTGGTCAAAGAGCATTCGGGCGCCTACGGCATGGGCGTCGa                                     <  2:469527/150‑1 (MQ=255)
                                                                                                               cATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTgg                                                                          <  1:677655/113‑1 (MQ=255)
                                                                                                               cATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTgg                                                                          >  2:677655/1‑113 (MQ=255)
                                                                                                               cATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCAAAGAGCATTCGGGCGCCTACGGCATGGGCGTCGa                                     >  2:435125/1‑150 (MQ=255)
                                                                                                               cATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCAAAGAGCATTCGGGCGCCTACGGCATGGGCGTCGa                                     >  1:1166940/1‑150 (MQ=255)
                                                                                                                aTGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTg                                                                                                                >  2:1545492/1‑74 (MQ=255)
                                                                                                                aTGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTg                                                                                                                <  1:1545492/74‑1 (MQ=255)
                                                                                                                   cAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCAAAGAGCATTCGGGCGCCTACGGCATGGGCGTCGACTAc                                 <  1:439758/150‑1 (MQ=255)
                                                                                                                   cAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCAAAGAGCATTCGGGCGCCTACGGCATGGGCGTCGACTAc                                 <  1:1754933/150‑1 (MQ=255)
                                                                                                                    aGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCaaa                                                                     >  2:1688274/1‑113 (MQ=255)
                                                                                                                    aGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCaaa                                                                     <  1:1688274/113‑1 (MQ=255)
                                                                                                                     gcgcAGGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATCGCGCTGGTCAAAGAGCATTCGGGCACCTACGGCATGGGCGTCGACTACGc                               <  1:446467/150‑1 (MQ=255)
                                                                                                                               cAGGCGGGCGTGGTGCCAGTTCCCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCAAANAGCATTCGGNCGCCTACGNGATGGGCGTCGACTACGCCTATACCATg                     <  1:543985/150‑1 (MQ=255)
                                                                                                                                  gcgggcgTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCAAAGAGCATTCGGGCGCCTACGGCATGGGCGTTGACTACGCCTATACCATGGTg                  <  1:238965/150‑1 (MQ=255)
                                                                                                                                      acctgctgCCAGGTAACTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCAAAGAGCATTCGGGCGCCTACGGCATGGGCGTCGACTACGCCTATACCATGGTGcaca              >  1:942538/7‑150 (MQ=255)
                                                                                                                                         tgctgCCAGGTAACTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCAAAGAGCATTCGGGCGCCTACGGCATGGGCGTc                                       >  1:706825/4‑122 (MQ=255)
                                                                                                                                         tgctgCCAGGTAACTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCAAAGAGCATTCGGGCGCCTACGGCATGGGCGTc                                       <  2:706825/119‑1 (MQ=255)
                                                                                                                                          gctgCCAGGTAACTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCAAAGAGCATTCGGGCGCCTACGGCATGGGCGTCGACTACGCCTATACCATGGTGCACAAAgc          <  2:942538/148‑1 (MQ=255)
                                                                                                                                            tgCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGCAAACCCGCGACACCTACGATGCTGTCATGGCGCTGGTCAAAGAGCATTCGGGCGCCTACGGCATGGGCGTCGACTACGCCTATACCATGGTGCACAAAgcgc        <  2:52291/150‑1 (MQ=255)
                                                                                                                                             gCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCAAAGAGCATTCGGGCGCCTACGGCATGGGCGTCGACTACGCCTATACCATGGTGCACAAAGCGcc       <  1:139705/150‑1 (MQ=255)
                                                                                                                                                aGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCAAAGAGCATTCGGGCGCCTACGGCATGGGCGTCGACTACGCCTATACCATGGTGCACAAAGCGCCaga    >  2:1138605/1‑150 (MQ=255)
                                                                                                                                                aGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCAAAGAGCATTCGGGCGCCTACGGCATAAGCCTCGACTACGCCTATACCATGGTGCACAAAGCGCCaga    >  1:990997/1‑150 (MQ=255)
                                                                                                                                                 gTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCAAag                                                                    >  1:1310243/1‑85 (MQ=255)
                                                                                                                                                 gTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCAAag                                                                    <  2:1310243/85‑1 (MQ=255)
                                                                                                                                                  ttACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCAAAGAGCATTCGGGCGCCTACGGCATGGGCGTCGACTACGCCTATACCATGGTGCACAAAGCGCCAGAgc  >  1:423775/1‑150 (MQ=255)
                                                                                                                                                  ttACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGATACGTTTTTGGTCAAAGAGCATTCGGGCGCCTACGGCATGGGCGTCGACTACGCCTATACCATGGTGCACAAAGCGCCAGAgc  <  1:337084/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
AGAGGTCTGCAACACCACCTTCGCGCTGTCCGCCATGCACGATGCAGGCTATGAGATCTATATGGTTGCCGACGCTTCCGGGGGCACGACCAAGGAGGCTCACGACTACGCCATGCAGCGCATGATCCAGGCGGGCGTGGTGCCAGTTACCTGGCAGCAGGTGCTACTGGAGTGGCAGCGCGACTGGAAAAACCGCGACACCTACGATGCTGTCATGGCGCTGGTCAAAGAGCATTCGGGCGCCTACGGCATGGGCGTCGACTACGCCTATACCATGGTGCACAAAGCGCCAGAGC  >  NC_002947/4484058‑4484353

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: