Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,855,065 G→T 5.8% A27A (GCC→GCA 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,0650GT 5.8% 297.3 / 3.9 95A27A (GCC→GCAdnaXDNA polymerase III subunit gamma
Reads supporting (aligned to +/- strand):  ref base G (46/43);  new base T (2/4);  total (48/47)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 101 total aligned reads displayed.

GAACAGGTGCCACAGGGCGTGGAGGTGATGCCGGTTTCACAGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTAGCGGAGCGGGGGCAAAATTGCACC  >  NC_002947/4854916‑4855214
                                                                                                                                                     |                                                                                                                                                     
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                                                                                                                                                       >  2:1550464/1‑143 (MQ=255)
       tGCCACAGGGCGTGGAGGTGATGCCGGTTTCACAGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAg                                                                                                                                                       <  1:1550464/143‑1 (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                                                                                                                                        <  2:2184422/150‑1 (MQ=255)
               ggCGTGGAGGTGATGCCGGTTTCACAGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATgg                                                                                                                                        >  1:841577/1‑150 (MQ=255)
               ggCGTGGAGGTGATGCCGGTTTCACAGTTCAGGCACTTGGCGATAATACAGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATgg                                                                                                                                        <  2:1760378/150‑1 (MQ=255)
                     gAGGTGATGCCGGTTTCACAGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTgg                                                                                                                                  >  2:298908/1‑150 (MQ=255)
                      aGGTGATGCCGGTTTCACAGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGc                                                                                                                                 <  2:1417343/150‑1 (MQ=255)
                             gCCGGTTTCACAGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCAt                                                                                                                          >  1:1914918/1‑150 (MQ=255)
                                      acaGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATgg                                                                                                                                        >  2:1034045/1‑127 (MQ=255)
                                      acaGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATgg                                                                                                                                        <  1:1034045/127‑1 (MQ=255)
                                          ttCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAg                                                                                                             >  2:429058/1‑150 (MQ=255)
                                            cAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCtt                                                                                                                                                  <  2:257207/111‑1 (MQ=255)
                                            cAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCtt                                                                                                                                                  >  1:257207/1‑111 (MQ=255)
                                                cACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGa                                                                                                       <  2:2149322/150‑1 (MQ=255)
                                                   ttGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGc                                                                                                    >  1:896572/1‑150 (MQ=255)
                                                   ttGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGc                                                                                                    >  2:2136597/1‑150 (MQ=255)
                                                         ataataCGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGAc                                                                                                      <  1:1687751/142‑1 (MQ=255)
                                                         ataataCGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGAc                                                                                                      >  2:1687751/1‑142 (MQ=255)
                                                          taataCGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGTTGGTTTTCCAAGGCATTGATCAATGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTc                                                                                             <  1:2121360/150‑1 (MQ=255)
                                                                  gCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTcgcg                                                                                                                    >  2:941167/1‑119 (MQ=255)
                                                                  gCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTcgcg                                                                                                                    <  1:941167/119‑1 (MQ=255)
                                                                    gATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAAc                                                                                   <  1:2136597/150‑1 (MQ=255)
                                                                         tcgtcTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCCTTGATCAATGCCtt                                                                                                                                                  <  2:2306446/82‑1 (MQ=255)
                                                                         tcgtcTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCCTTGATCAATGCCtt                                                                                                                                                  >  1:2306446/1‑82 (MQ=255)
                                                                         tcgtcTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGAt                                                                              <  2:460743/150‑1 (MQ=255)
                                                                         tcgtcTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGAt                                                                              >  1:648768/1‑150 (MQ=255)
                                                                           gtcTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATaa                                                                            >  1:1303539/1‑150 (MQ=255)
                                                                            tcTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCCTTGATCAATGCCtt                                                                                                                                                  <  2:1023814/79‑1 (MQ=39)
                                                                            tcTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCCTTGATCAATGCCtt                                                                                                                                                  >  1:1023814/1‑79 (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                                                                                                                                  >  2:2017243/1‑85 (MQ=255)
                                                                                      cccGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTgg                                                                                                                                  <  1:2017243/85‑1 (MQ=255)
                                                                                       ccGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCtt                                                                                                                                                  <  1:282728/68‑1 (MQ=255)
                                                                                       ccGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCtt                                                                                                                                                  >  2:282728/1‑68 (MQ=255)
                                                                                       ccGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGaaaaa                                                                >  2:1537671/1‑150 (MQ=255)
                                                                                       ccGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGaaaaa                                                                >  2:1470555/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                                                                                          <  2:292739/104‑1 (MQ=255)
                                                                                                           gTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTg                                                                                          >  1:292739/1‑104 (MQ=255)
                                                                                                           gTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGc                                            >  2:1035486/1‑150 (MQ=255)
                                                                                                           gTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGc                                            >  1:1564340/1‑150 (MQ=255)
                                                                                                           gTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGc                                            >  1:688265/1‑150 (MQ=255)
                                                                                                           gTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGc                                            >  2:1758618/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                                     >  2:612468/1‑147 (MQ=255)
                                                                                                                     tGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGAt                                     <  1:612468/147‑1 (MQ=255)
                                                                                                                          gcgcTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATg                             <  2:841577/150‑1 (MQ=255)
                                                                                                                             cTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTa                          <  1:1537671/150‑1 (MQ=255)
                                                                                                                               ggTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTagc                        <  1:504588/150‑1 (MQ=255)
                                                                                                                                gTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTagcg                       <  1:1835805/150‑1 (MQ=255)
                                                                                                                                gTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTagcg                       >  2:1275642/1‑150 (MQ=255)
                                                                                                                                      cAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGc                                                                                         >  1:2022816/1‑78 (MQ=255)
                                                                                                                                      cAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGc                                                                                         <  2:2022816/78‑1 (MQ=255)
                                                                                                                                      cAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTagcggagcgg                 <  1:980929/150‑1 (MQ=255)
                                                                                                                                       aaGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCAt                                                            <  1:219382/106‑1 (MQ=255)
                                                                                                                                       aaGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCAt                                                            >  2:219382/1‑106 (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           >  2:1343484/1‑150 (MQ=255)
                                                                                                                                            aTTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTAGCGGAGCGGGGGCaa           >  1:2335722/1‑150 (MQ=255)
                                                                                                                                            aTTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTAGCGGAGCGGGGGCaa           >  1:550289/1‑150 (MQ=255)
                                                                                                                                                aTCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTAGCGGAGCGGGGGCAAAAtt       >  1:186296/1‑150 (MQ=255)
                                                                                                                                                 tCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTAGCGGAGCGGGGGCAAAAtt       <  2:58893/149‑1 (MQ=255)
                                                                                                                                                 tCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTAGCGGAGCGGGGGCAAAAtt       >  1:58893/1‑149 (MQ=255)
                                                                                                                                                  cAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTAGCGGAGCGGGGGCAAAATTGc     >  2:105485/1‑150 (MQ=255)
                                                                                                                                                   aaGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTAGCGGAGCGGGGGCAAAATTGCa    >  1:673316/1‑150 (MQ=255)
                                                                                                                                                   aaGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTAGCGGAGCGGGGGCAAAATTGCa    >  1:1469187/1‑150 (MQ=255)
                                                                                                                                                   aaGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTAGCGGAGCGGGGGCAAAATTGCa    <  1:2012087/150‑1 (MQ=255)
                                                                                                                                                   aaGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTAGCGGAGCGGGGGCAAAATTGCa    <  1:1470555/150‑1 (MQ=255)
                                                                                                                                                    aGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTAGCGGAGCGGGGGCAAAATTGCAc   <  1:429058/150‑1 (MQ=255)
                                                                                                                                                     ggCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTAGCGGAGCGGGGGCAAAATTGCAcc  >  1:1671233/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAACAGGTGCCACAGGGCGTGGAGGTGATGCCGGTTTCACAGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTAGCGGAGCGGGGGCAAAATTGCACC  >  NC_002947/4854916‑4855214

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: