Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,855,056 A→C 5.2% N30K (AAT→AAG dnaX ← DNA polymerase III subunit gamma

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,855,0560AC 5.2% 282.6 / 6.0 97N30K (AAT→AAGdnaXDNA polymerase III subunit gamma
Reads supporting (aligned to +/- strand):  ref base A (45/47);  new base C (2/3);  total (47/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.53e-01

GGCACACCGAACAGGTGCCACAGGGCGTGGAGGTGATGCCGGTTTCACAGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTAGCGGAGCGGGGGCAA  >  NC_002947/4854908‑4855205
                                                                                                                                                    |                                                                                                                                                     
ggCACACCGAACAGGTGCCACAGGGCGTGGAGGTGATGCCGGTTTCACAGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCAt                                                                                                                                                      <  2:1202032/150‑1 (MQ=255)
ggCACACCGAACAGGTGCCACAGGGCGTGGAGGTGATGCCGGTTTCACAGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCAt                                                                                                                                                      >  1:852507/1‑150 (MQ=255)
 gCACACCGAACAGGTGCCACAGGGCGTGGAGGTGATGCCGGTTTCACAGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCAtt                                                                                                                                                     <  2:415683/150‑1 (MQ=255)
   acacCGAACAGGTGCCACAGGGCGTGGAGGTGATGCCGGTTTCACAGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGa                                                                                                                                                   <  1:276657/150‑1 (MQ=255)
     acCGAACAGGTGCCACAGGGCGTGGAGGTGATGCCGGTTTCACAGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATc                                                                                                                                                 >  1:1417343/1‑150 (MQ=255)
        gAACAGGTGCCACAGGGCGCGGAGGTGATGCCGGTTTCACAGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAg                                                                                                                                              >  2:341813/1‑150 (MQ=255)
            aGGTGCCACAGGGCGTGGAGGTGATGCCGGTTTCACAGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCt                                                                                                                                          <  2:852507/150‑1 (MQ=255)
             ggTGCCACAGGGCGTGGAGGTGATGCCGGTTTCACAGTGCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCtt                                                                                                                                         >  1:460743/1‑150 (MQ=255)
              gTGCCACAGGGCGTGGAGGTGATGCCGGTTTCACAGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCtt                                                                                                                                         >  2:1504665/1‑149 (MQ=255)
              gTGCCACAGGGCGTGGAGGTGATGCCGGTTTCACAGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCtt                                                                                                                                         <  1:1504665/149‑1 (MQ=255)
               tGCCACAGGGCGTGGAGGTGATGCCGGTTTCACAGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAg                                                                                                                                              <  1:1550464/143‑1 (MQ=255)
               tGCCACAGGGCGTGGAGGTGATGCCGGTTTCACAGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAg                                                                                                                                              >  2:1550464/1‑143 (MQ=255)
               tGCCACAGGGCGTGGAGGTGATGCCGGTTTCACAGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGa                                                                                                                                       <  1:2128835/150‑1 (MQ=255)
               tGCCACAGGGCGTGGAGGTGATGCCGGTTTCACAGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGa                                                                                                                                       <  2:2156169/150‑1 (MQ=255)
                  cacaGGGCGTGGAGGTGATGCCGGTTTCACAGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGca                                                                                                                                    >  2:2012087/1‑150 (MQ=255)
                    caGGGCGTGGAGGTGATGCCGGTTTCACAGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCCTTGATCAATGCCTTggacaac                                                                                                                                  <  1:1586564/150‑7 (MQ=255)
                      gggCGTGGAGGTGATGCCGGTTTCACAGTTGGGGCACTGGGCGATACTACGGGCGATGGTCCTCTTGTGGCACGCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATg                                                                                                                                <  1:1223911/150‑1 (MQ=255)
                       ggCGTGGAGGTGATGCCGGTTTCACAGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATgg                                                                                                                               >  1:841577/1‑150 (MQ=255)
                       ggCGTGGAGGTGATGCCGGTTTCACAGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATgg                                                                                                                               <  2:2184422/150‑1 (MQ=255)
                       ggCGTGGAGGTGATGCCGGTTTCACAGTTCAGGCACTTGGCGATAATACAGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATgg                                                                                                                               <  2:1760378/150‑1 (MQ=255)
                             gAGGTGATGCCGGTTTCACAGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCAtt                                                                                                                                                     <  1:1622573/122‑1 (MQ=255)
                             gAGGTGATGCCGGTTTCACAGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCAtt                                                                                                                                                     >  2:1622573/1‑122 (MQ=255)
                             gAGGTGATGCCGGTTTCACAGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTgg                                                                                                                         >  2:298908/1‑150 (MQ=255)
                              aGGTGATGCCGGTTTCACAGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGc                                                                                                                        <  2:1417343/150‑1 (MQ=255)
                                     gCCGGTTTCACAGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCAt                                                                                                                 >  1:1914918/1‑150 (MQ=255)
                                              acaGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTg                                                                                                                                                    >  2:730465/1‑106 (MQ=255)
                                              acaGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTg                                                                                                                                                    <  1:730465/106‑1 (MQ=255)
                                              acaGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATgg                                                                                                                               >  2:1034045/1‑127 (MQ=255)
                                              acaGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATgg                                                                                                                               <  1:1034045/127‑1 (MQ=255)
                                                  ttCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAg                                                                                                    >  2:429058/1‑150 (MQ=255)
                                                    cAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCtt                                                                                                                                         >  1:257207/1‑111 (MQ=255)
                                                    cAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCtt                                                                                                                                         <  2:257207/111‑1 (MQ=255)
                                                        cACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGa                                                                                              <  2:2149322/150‑1 (MQ=255)
                                                           ttGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGc                                                                                           >  2:2136597/1‑150 (MQ=255)
                                                           ttGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGc                                                                                           >  1:896572/1‑150 (MQ=255)
                                                                 ataataCGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGAc                                                                                             >  2:1687751/1‑142 (MQ=255)
                                                                 ataataCGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGAc                                                                                             <  1:1687751/142‑1 (MQ=255)
                                                                  taataCGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGTTGGTTTTCCAAGGCATTGATCAATGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTc                                                                                    <  1:2121360/150‑1 (MQ=255)
                                                                     taCGGGCGATGGTAGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTg                                                                                 >  2:980929/1‑150 (MQ=255)
                                                                          gCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTcgcg                                                                                                           >  2:941167/1‑119 (MQ=255)
                                                                          gCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTcgcg                                                                                                           <  1:941167/119‑1 (MQ=255)
                                                                            gATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAAc                                                                          <  1:2136597/150‑1 (MQ=255)
                                                                                 tcgtcTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCCTTGATCAATGCCtt                                                                                                                                         >  1:2306446/1‑82 (MQ=255)
                                                                                 tcgtcTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCCTTGATCAATGCCtt                                                                                                                                         <  2:2306446/82‑1 (MQ=255)
                                                                                 tcgtcTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGAt                                                                     >  1:648768/1‑150 (MQ=255)
                                                                                 tcgtcTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGAt                                                                     <  2:460743/150‑1 (MQ=255)
                                                                                   gtcTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATaa                                                                   >  1:1303539/1‑150 (MQ=255)
                                                                                    tcTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCCTTGATCAATGCCtt                                                                                                                                         >  1:1023814/1‑79 (MQ=39)
                                                                                    tcTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCCTTGATCAATGCCtt                                                                                                                                         <  2:1023814/79‑1 (MQ=39)
                                                                                      ttGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTc                                                                <  2:1914918/150‑1 (MQ=255)
                                                                                         cccACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATc                                                             <  2:353052/150‑1 (MQ=255)
                                                                                         cccACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTGCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATc                                                             >  1:133927/1‑150 (MQ=255)
                                                                                          ccACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCg                                                            >  1:210231/1‑150 (MQ=255)
                                                                                           cACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACAt                                                                                                                                 <  1:1724409/80‑1 (MQ=255)
                                                                                           cACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACAt                                                                                                                                 >  2:1724409/1‑80 (MQ=255)
                                                                                            aCCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAg                                                                                                                                              <  1:232554/66‑1 (MQ=255)
                                                                                            aCCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAg                                                                                                                                              >  2:232554/1‑66 (MQ=255)
                                                                                              cccGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTgg                                                                                                                         <  1:2017243/85‑1 (MQ=255)
                                                                                              cccGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTgg                                                                                                                         >  2:2017243/1‑85 (MQ=255)
                                                                                               ccGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCtt                                                                                                                                         <  1:282728/68‑1 (MQ=255)
                                                                                               ccGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCtt                                                                                                                                         >  2:282728/1‑68 (MQ=255)
                                                                                               ccGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGaaaaa                                                       >  2:1470555/1‑150 (MQ=255)
                                                                                               ccGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGaaaaa                                                       >  2:1537671/1‑150 (MQ=255)
                                                                                                   gggTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGATCCTGCGCGGATGCCATTTTCGTTCAAGAACCTGATAACTCATCGAAAAACCAt                                                   <  1:298908/150‑1 (MQ=255)
                                                                                                       aCCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGACATGGAGTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTa                                               >  1:2251432/1‑150 (MQ=255)
                                                                                                         cGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGc                                             <  2:648768/150‑1 (MQ=255)
                                                                                                           gTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCt                                           >  2:1527224/1‑150 (MQ=255)
                                                                                                               aCAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATc                                       >  2:1887957/1‑150 (MQ=255)
                                                                                                                  ggTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAg                                    <  1:1887957/150‑1 (MQ=255)
                                                                                                                   gTAGTCATGGTTCAGGCGCTGGTTGTCAAACGCATTGATCAAGTCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGCAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGc                                   <  1:821115/150‑1 (MQ=255)
                                                                                                                   gTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTg                                                                                 >  1:292739/1‑104 (MQ=255)
                                                                                                                   gTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTg                                                                                 <  2:292739/104‑1 (MQ=255)
                                                                                                                   gTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGc                                   >  2:1758618/1‑150 (MQ=255)
                                                                                                                   gTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGc                                   >  1:688265/1‑150 (MQ=255)
                                                                                                                   gTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGc                                   >  1:1564340/1‑150 (MQ=255)
                                                                                                                   gTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGc                                   >  2:1035486/1‑150 (MQ=255)
                                                                                                                       gCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGTCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGaa                               <  1:341813/150‑1 (MQ=255)
                                                                                                                         aTGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTg                                                                                 <  2:288602/98‑1 (MQ=255)
                                                                                                                         aTGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTg                                                                                 >  1:288602/1‑98 (MQ=255)
                                                                                                                            gTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCg                          >  2:1865241/1‑150 (MQ=255)
                                                                                                                             tGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGAt                            <  1:612468/147‑1 (MQ=255)
                                                                                                                             tGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGAt                            >  2:612468/1‑147 (MQ=255)
                                                                                                                                  gcgcTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATg                    <  2:841577/150‑1 (MQ=255)
                                                                                                                                     cTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTa                 <  1:1537671/150‑1 (MQ=255)
                                                                                                                                       ggTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTagc               <  1:504588/150‑1 (MQ=255)
                                                                                                                                        gTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTagcg              >  2:1275642/1‑150 (MQ=255)
                                                                                                                                        gTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTagcg              <  1:1835805/150‑1 (MQ=255)
                                                                                                                                              cAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGc                                                                                >  1:2022816/1‑78 (MQ=255)
                                                                                                                                              cAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGc                                                                                <  2:2022816/78‑1 (MQ=255)
                                                                                                                                              cAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTagcggagcgg        <  1:980929/150‑1 (MQ=255)
                                                                                                                                               aaGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCAt                                                   >  2:219382/1‑106 (MQ=255)
                                                                                                                                               aaGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCAt                                                   <  1:219382/106‑1 (MQ=255)
                                                                                                                                               aaGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTAGCGGAGCggg       <  1:1275642/150‑1 (MQ=255)
                                                                                                                                                aGGCATTGATCAAGGCCTTGAGCACCTCGGTCTGTCCGACCATTTCGCGGAACGAGCGCCGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTAGCGGAGCgggg      <  2:896572/150‑1 (MQ=255)
                                                                                                                                                aGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGa                                                                      <  1:1541611/86‑1 (MQ=255)
                                                                                                                                                aGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGa                                                                      >  2:1541611/1‑86 (MQ=255)
                                                                                                                                                 ggCATGGATCAAGGCGGTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAGCTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTAGCGGAGCggggg     <  2:699044/150‑1 (MQ=255)
                                                                                                                                                    aTTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTAGCGGAGCGGGGGCaa  >  1:550289/1‑150 (MQ=255)
                                                                                                                                                    aTTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTAGCGGAGCGGGGGCaa  >  1:2335722/1‑150 (MQ=255)
                                                                                                                                                    aTTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTAGCGGAGCGGGGGCaa  >  2:1343484/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGCACACCGAACAGGTGCCACAGGGCGTGGAGGTGATGCCGGTTTCACAGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTAGCGGAGCGGGGGCAA  >  NC_002947/4854908‑4855205

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