Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,864,795 G→A 6.2% A29A (GCG→GCA hda → DnaA regulatory inactivator Hda

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,864,7950GA 6.2% 241.2 / 4.4 81A29A (GCG→GCAhdaDnaA regulatory inactivator Hda
Reads supporting (aligned to +/- strand):  ref base G (36/40);  new base A (2/3);  total (38/43)
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.74e-01

GGGTATAGACTTTGCAGATTGTTCACCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCAGCAACGTGGCGAGCCCGCTGT  >  NC_002947/1864646‑1864941
                                                                                                                                                     |                                                                                                                                                  
gggTATAGACTTTGCAGATTGTTCACCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCg                                                                                                                                                    <  1:1268272/150‑1 (MQ=255)
  gTATAGACTTTGCAGATTGTTCACCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGc                                                                                                                                                  >  1:931481/1‑150 (MQ=255)
     taGACTTTGCCGATTGTTCACCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCAtt                                                                                                                                               <  1:291991/150‑1 (MQ=255)
      aGACTTTGCAGATTGTTCACCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTg                                                                                                                                              <  2:931481/150‑1 (MQ=255)
      aGACTTTGCAGATTGTTCACCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCCGCAGCATTg                                                                                                                                              <  1:1389934/150‑1 (MQ=255)
        aCTTTGCAGATTGTTCACCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTggg                                                                                                                                            <  2:1719424/150‑1 (MQ=255)
          tttGCAGATTGTTCACCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGtt                                                                                                                                          >  1:1323977/1‑150 (MQ=255)
           ttGCAGATTGTTCACCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTa                                                                                                                                         <  2:1080984/150‑1 (MQ=255)
            tGCAGATTGTTCACCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCCGCAGCATTggcgccc                                                                                                                                        <  2:1879791/150‑6 (MQ=255)
                gATTGTTCACCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCg                                                                                                                                    >  1:419009/1‑150 (MQ=255)
                        aCCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTa                                                                                                                            >  1:1284097/1‑150 (MQ=255)
                        aCCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTa                                                                                                                            >  2:1771638/1‑150 (MQ=255)
                        aCCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCATGCGGGGTGATGCAACCTTCATCAACTACTATCCGGGCGCCCATGCTGCGGCATTGGGTTACGTCGAGCGGCTa                                                                                                                            >  1:1468006/1‑150 (MQ=255)
                         ccAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTAt                                                                                                                           >  1:1215631/1‑150 (MQ=255)
                          cAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCg                                                                                                                                                    <  1:1248665/124‑1 (MQ=255)
                          cAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCg                                                                                                                                                    >  2:1248665/1‑124 (MQ=255)
                          cAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATg                                                                                                                          >  1:223279/1‑150 (MQ=255)
                          cAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATg                                                                                                                          >  2:411617/1‑150 (MQ=255)
                          cAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATg                                                                                                                          >  1:427139/1‑150 (MQ=255)
                           aaaGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTggg                                                                                                                                            <  2:1191479/131‑1 (MQ=255)
                           aaaGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTggg                                                                                                                                            >  1:1191479/1‑131 (MQ=255)
                              ggCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGTGTTACGTCGAGCGGCTATGCGaa                                                                                                                      <  1:258422/150‑1 (MQ=255)
                                cccccAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGc                                                                                                                    <  2:707915/150‑1 (MQ=255)
                                 ccccAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGaa                                                                                                                      >  2:1200556/1‑147 (MQ=255)
                                 ccccAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGaa                                                                                                                      <  1:1200556/147‑1 (MQ=255)
                                 ccccAGCGGTCCCATCGACCGCCTGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGcc                                                                                                                   >  2:941045/1‑150 (MQ=255)
                                  cccAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTgg                                                                                                                                             >  1:308597/1‑123 (MQ=255)
                                  cccAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTgg                                                                                                                                             <  2:308597/123‑1 (MQ=255)
                                          tCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCACCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGc                                                                                                          <  1:502551/150‑1 (MQ=255)
                                             cGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTgg                                                                                                       >  2:1362447/1‑150 (MQ=255)
                                                  aCCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCga                                                                                                  <  2:1563055/150‑1 (MQ=255)
                                                    cgcccgcGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCgaga                                                                                                >  2:621582/1‑150 (MQ=255)
                                                     gcccgcGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCgaga                                                                                                <  1:1331615/149‑1 (MQ=255)
                                                     gcccgcGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCgaga                                                                                                >  2:1331615/1‑149 (MQ=255)
                                                      cccgcGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGc                                                                                              <  2:419009/150‑1 (MQ=255)
                                                       ccgcGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATGCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGcc                                                                                             >  2:1313362/1‑150 (MQ=255)
                                                         gcgAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCtt                                                                                           <  1:941045/150‑1 (MQ=255)
                                                         gcgAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCtt                                                                                           <  2:1760627/150‑1 (MQ=255)
                                                           gAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTAt                                                                                         <  1:555505/150‑1 (MQ=255)
                                                            aGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCATGGGCGGTCTCACCACCTTTATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATc                                                                                        <  2:50560/150‑1 (MQ=255)
                                                             gCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCt                                                                                       >  2:768455/1‑150 (MQ=255)
                                                               aTGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGATGCCGTCGCCGGCTGGACCGAGAGCCTTATCTAc                                                                                     <  1:250422/150‑1 (MQ=255)
                                                               aTGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTAc                                                                                     <  2:1736910/150‑1 (MQ=255)
                                                                   aaCCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCtgt                                                                                 >  1:1424350/1‑150 (MQ=255)
                                                                     ccaccGATCCTGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTgg                                                                               <  1:1612313/150‑1 (MQ=255)
                                                                     ccaccGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGGGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTgg                                                                               <  1:1149209/150‑1 (MQ=255)
                                                                         cGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCa                                                                           <  2:1323977/150‑1 (MQ=255)
                                                                           aTCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAg                                                                         <  1:416855/150‑1 (MQ=255)
                                                                             ccAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCa                                                                       <  1:488465/150‑1 (MQ=255)
                                                                                    ccccTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGtt                                                                <  2:1470243/150‑1 (MQ=255)
                                                                                      ccTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTgg                                                              >  2:1909952/1‑150 (MQ=255)
                                                                                       cTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCgg                                                                                                                               <  1:1232301/84‑1 (MQ=255)
                                                                                       cTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCgg                                                                                                                               >  2:1232301/1‑84 (MQ=255)
                                                                                        tGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCt                                                                                   >  1:1612442/1‑127 (MQ=255)
                                                                                        tGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCt                                                                                   <  2:1612442/127‑1 (MQ=255)
                                                                                              tgCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCACGCAGGGTGTTGGCCGTTccc                                                      <  1:963637/150‑1 (MQ=255)
                                                                                                       gcgATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGtt                                                                <  1:356240/131‑1 (MQ=255)
                                                                                                       gcgATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGtt                                                                >  2:356240/1‑131 (MQ=255)
                                                                                                        cgATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCa                                            >  2:1008378/1‑150 (MQ=255)
                                                                                                         gATGACGCCAGCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGATAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAg                                           <  1:769553/150‑1 (MQ=255)
                                                                                                               gCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCGAGCAGTGTGTTGGCCGTTTCCACCTGCTGCAGgccgcc                                     >  1:927796/1‑150 (MQ=255)
                                                                                                                         tcaACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCt                           >  2:1676303/1‑150 (MQ=255)
                                                                                                                          caACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCtt                          >  1:925668/1‑150 (MQ=255)
                                                                                                                             ctactaTCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCa                       <  1:621582/150‑1 (MQ=255)
                                                                                                                             ctactaTCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCa                       <  2:723990/150‑1 (MQ=255)
                                                                                                                             ctaattaccGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCa                       <  1:465085/143‑1 (MQ=255)
                                                                                                                              tactaTCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCAg                      <  1:1909952/150‑1 (MQ=255)
                                                                                                                               actaTCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCAGc                     <  1:768455/150‑1 (MQ=255)
                                                                                                                                 taTCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCCACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCAGCaa                   >  2:643515/1‑150 (MQ=255)
                                                                                                                                   tCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGAACATACGTTGCGGACGCCGTCTGTACAGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCAGCAACg                 >  2:66002/1‑150 (MQ=255)
                                                                                                                                     cGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCGGATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCAGCAACGTg               <  1:1771638/150‑1 (MQ=255)
                                                                                                                                       ggCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCAGCAACGTGGc             <  2:1561958/150‑1 (MQ=255)
                                                                                                                                          gcCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCAGCAACGTGGCGAg          >  1:1014292/1‑150 (MQ=255)
                                                                                                                                          cccAATGCCGCAGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCtgtg                                                                                <  1:774548/79‑1 (MQ=255)
                                                                                                                                          cccAATGCCGCAGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCtgtg                                                                                >  2:774548/2‑80 (MQ=255)
                                                                                                                                           ccAATGCCGCAGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCAGCAACGTGGCGAGc         >  2:659683/1‑150 (MQ=255)
                                                                                                                                            cAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCAGCAACGTGGCGAGcc        >  2:256746/1‑150 (MQ=255)
                                                                                                                                             aaTGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCAGCAACGTGGCGAGccc       >  1:421577/1‑150 (MQ=255)
                                                                                                                                             aaTGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAAAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCAGCAACGTGGCGAGccc       <  1:411617/150‑1 (MQ=255)
                                                                                                                                                  tgcGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCtt                          >  1:957077/1‑126 (MQ=255)
                                                                                                                                                  tgcGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCtt                          <  2:957077/126‑1 (MQ=255)
                                                                                                                                                  tgcGGCATTGGGTTACGTCGAGCGGATATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCa                                    <  2:1810229/116‑1 (MQ=255)
                                                                                                                                                  tgcGGCATTGGGTTACGTCGAGCGGATATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCa                                    >  1:1810229/1‑116 (MQ=255)
                                                                                                                                                  tgcGGCATTGGGCTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCAGCAACGTGGCGAGCCCGCTGt  >  1:1327174/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GGGTATAGACTTTGCAGATTGTTCACCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCAGCAACGTGGCGAGCCCGCTGT  >  NC_002947/1864646‑1864941

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