Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,928,771 G→A 5.9% R112C (CGC→TGC)  fbpA ← iron ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,928,7710GA 5.9% 246.4 / 4.6 84R112C (CGC→TGC) fbpAiron ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (41/38);  new base A (3/2);  total (44/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.82e-01

CGCAGGCACAGACGGCCTTCCCACTGTTTGTCAGCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGACGGTCAGCAGCAGGTCGGCCACGCCGTTATCGCCCTCGGCCTTGATGCGCTGC  >  NC_002947/5928632‑5928916
                                                                                                                                            |                                                                                                                                                 
cgcAGGCACAGACGGCCTTCCCACTGTTTGTCAGCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGTTAATTGGAGTAGATGCTGACGCGGATCTCGTATGCCGTGTCGTGCTAGTCATGGGACGAGGCGCGGTACTGCGcc                                                                                                                                         >  1:1375310/1‑150 (MQ=255)
 gcAGGCACAGACGGCCTTCCCACTGTTTGTCAGCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCg                                                                                                                                        >  2:1557884/1‑150 (MQ=255)
    ggCACAGACGGCCTTCCCAGCGTTTGTCAGCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGa                                                                                                                                     <  2:125034/150‑1 (MQ=255)
     gCACAGACGGCCTTCCCACTGTTTGTCAGCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGAGCGGGCGAGCAGGCTGAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGAt                                                                                                                                    >  1:881813/1‑150 (MQ=255)
       acaGACGGCCTTCCCACTGTTTGTCAGCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTgc                                                                                                                                            >  1:189105/1‑140 (MQ=255)
       acaGACGGCCTTCCCACTGTTTGTCAGCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTgc                                                                                                                                            <  2:189105/140‑1 (MQ=255)
            cGGCCTTCCCACTGTTTGTCAGCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTgg                                                                                                                             <  2:1832156/150‑1 (MQ=255)
             ggCCTTCCCCTTGTTTGTCAGCCAGGGCTTCGTAGGTGCTCACCTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGt                                                                                                                            <  2:54385/150‑1 (MQ=255)
             ggCCTTCCCACTGTTTGTCAGCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCt                                                                                                                               >  1:1607297/1‑147 (MQ=255)
             ggCCTTCCCACTGTTTGTCAGCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGt                                                                                                                            >  2:777766/1‑150 (MQ=255)
             ggCCTTCCCACTGTTTGTCAGCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGt                                                                                                                            <  1:402063/150‑1 (MQ=255)
             ggCCTTCCCACTGTTTGTCAGCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGt                                                                                                                            >  2:1144946/1‑150 (MQ=255)
             ggCCTTCCCACTGTTTGTCAGCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGt                                                                                                                            >  2:402063/1‑150 (MQ=255)
             ggCCTTCCCACTGTTTGTCAGCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGGACTGCGCCGGGATGTTCt                                                                                                                               <  2:1607297/147‑1 (MQ=255)
              gCCTTCCCACTGTTTGTCAGCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTc                                                                                                                           <  2:418102/150‑1 (MQ=255)
                cTTCCCCCTGTTTGACAGCTAGGGTTTCGGAGGTGCTGAACTCCTCGGGTTTGAGCTGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGa                                                                                                                         <  2:116460/150‑1 (MQ=255)
                cTTCCCACTGTTTGTCAGCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTa                                                                                                                                                <  2:1702449/127‑1 (MQ=255)
                cTTCCCACTGTTTGTCAGCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTa                                                                                                                                                >  1:1702449/1‑127 (MQ=255)
                cTTCCCACTGTTTGTCAGCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGa                                                                                                                         >  2:350583/1‑150 (MQ=255)
                cTTCCCACTGTTTGTCAGCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGa                                                                                                                         <  2:1061502/150‑1 (MQ=255)
                 ttCCCACTGTTTGTCAGCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGAt                                                                                                                        <  1:665382/150‑1 (MQ=255)
                    ccACTGTTTGTCAGCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGAc                                                                                                                     <  1:1144946/150‑1 (MQ=255)
                    ccACTGTTTGTCAGCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGAc                                                                                                                     <  1:1899171/150‑1 (MQ=255)
                     cACTGTTTGTCAGCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACg                                                                                                                    >  2:1874453/1‑150 (MQ=255)
                      aCTGTTTGTCAGCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGt                                                                                                                   <  1:777766/150‑1 (MQ=255)
                          tttGTCAGCCAGGGCTTCGTAGGTGCTCGAATCCTCGGGTTTGACCCGCTCGGTGGGGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGa                                                                                                               >  2:999072/1‑150 (MQ=255)
                          tttGTCAGCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACAAGATCGGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGa                                                                                                               >  1:1453251/1‑150 (MQ=255)
                             gTCAGCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGGTGTTCTGGTCGATGACGTCGGACTc                                                                                                            >  1:1167357/1‑150 (MQ=255)
                             gTCAGCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCAGTACCGCGcctc                                                                                                                                       >  2:1459605/1‑121 (MQ=255)
                             gTCAGCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCAGTACCGCGcctc                                                                                                                                       <  1:1459605/123‑3 (MQ=255)
                                aGCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCCCAGGCCGGTCCAGTAATGGGACGAGGCGCGGTACTGTGTCGGGATGTTCTGGTCGATGACGTCGGACTCGAt                                                                                                         <  1:1874453/150‑1 (MQ=255)
                                 gCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATc                                                                                                        >  2:1924732/1‑150 (MQ=255)
                                   cAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCgg                                                                                                      >  1:546685/1‑150 (MQ=255)
                                        cTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGACTCGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCa                                                                                                 <  1:1645781/150‑1 (MQ=255)
                                         ttCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCGGCAg                                                                                                >  1:1665989/1‑150 (MQ=255)
                                             tAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGATCCGGTCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATg                                                                                            <  1:1057880/150‑1 (MQ=255)
                                                  gCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCAt                                                                                       <  2:1375310/150‑1 (MQ=255)
                                                   cTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATc                                                                                      >  2:1617809/1‑150 (MQ=255)
                                                    tCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCt                                                                                     >  2:577516/1‑150 (MQ=255)
                                                         tcctcGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGAt                                                                                                         >  2:1910229/1‑125 (MQ=255)
                                                         tcctcGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGAt                                                                                                         <  1:1910229/125‑1 (MQ=255)
                                                           ctcGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGc                                                                              <  1:999072/150‑1 (MQ=255)
                                                           ctcGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGc                                                                              <  2:419249/150‑1 (MQ=255)
                                                            tcGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCgg                                                                                                                >  2:433951/1‑115 (MQ=255)
                                                            tcGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCgg                                                                                                                <  1:433951/115‑1 (MQ=255)
                                                            tcGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGcc                                                                             >  1:1577278/1‑150 (MQ=255)
                                                            tcGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCAGTACCGCGcctcgt                                                                                                                                     <  1:1478245/94‑5 (MQ=255)
                                                            tcGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCAGTACCGCGcctcgt                                                                                                                                     >  2:1478245/1‑90 (MQ=255)
                                                             cGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCt                                                                            <  2:813427/150‑1 (MQ=255)
                                                                 tttGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGa                                                                                                                         >  2:1682141/1‑101 (MQ=255)
                                                                 tttGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGa                                                                                                                         <  1:1682141/101‑1 (MQ=255)
                                                                    taCCTGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAg                                                                     >  2:968259/2‑150 (MQ=255)
                                                                    gACCCGCTCGGTGGAGGAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCGGCTCGGCCGGCCACAg                                                                     >  2:974514/1‑150 (MQ=255)
                                                                      cccGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGt                                                                   >  1:1636841/1‑150 (MQ=255)
                                                                       ccGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGtt                                                                  >  1:702159/1‑150 (MQ=255)
                                                                       ccGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGtt                                                                  <  2:702159/150‑1 (MQ=255)
                                                                          cTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGcc                                                               <  1:350583/150‑1 (MQ=255)
                                                                           tCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCt                                                              <  2:546685/150‑1 (MQ=255)
                                                                              gTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCa                                                           <  1:968259/150‑1 (MQ=255)
                                                                                   gTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGAc                                                      <  2:632039/150‑1 (MQ=255)
                                                                                       atgatgGTGCGGGCGCGCAGGCTCAGGCCGATCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGACGGTc                                                  >  2:805176/1‑150 (MQ=255)
                                                                                                    gcgcgcAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGACGGTCAGCAGCAGGTCg                                      >  1:829193/1‑150 (MQ=255)
                                                                                                    gcgcgcAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGACGGTCAGCAGCAGGTCg                                      <  2:378047/150‑1 (MQ=255)
                                                                                                          aGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAAt                                                                                             >  1:390659/1‑89 (MQ=255)
                                                                                                          aGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAAt                                                                                             <  2:390659/89‑1 (MQ=255)
                                                                                                          aGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGACGGTCAGCAGCAGGTCGGCCACg                                >  1:167783/1‑150 (MQ=255)
                                                                                                              tCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGACGGTcagcagcag                                          >  2:899744/1‑136 (MQ=255)
                                                                                                              tCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGACGGTcagcagcag                                          <  1:899744/136‑1 (MQ=255)
                                                                                                              tCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGACGGTCAGCAGCAGGTCGGCCACGCCGt                            >  2:1624238/1‑150 (MQ=255)
                                                                                                                  gCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGc                                                                          >  2:213381/1‑100 (MQ=255)
                                                                                                                  gCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGc                                                                          <  1:213381/100‑1 (MQ=255)
                                                                                                                  gCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGACGGTCAGCAGCAGGTCGGCCACGCCGTTATc                        <  2:1636841/150‑1 (MQ=255)
                                                                                                                     ggTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGACGGTCAGCAGCAGGTCGGCCACGCCGTTATCGcc                     >  2:386268/1‑150 (MQ=255)
                                                                                                                       tCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGACGGTCAGCAGCAGGTCGGCCACGCCGTTATCGCCCt                   <  1:1924732/150‑1 (MQ=255)
                                                                                                                              aTGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGCCTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTTCATCGACGGTCAGCAGCAGGTCGGCCACGCCGTTATCGCCCTCGGCCtt            <  2:1597394/150‑1 (MQ=255)
                                                                                                                               tGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGACGGTCAGCAGCAGGTCGGCCACGCCGTTATCGCCCTCg                 <  2:1578877/144‑1 (MQ=255)
                                                                                                                               tGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGACGGTCAGCAGCAGGTCGGCCACGCCGTTATCGCCCTCg                 >  1:1578877/1‑144 (MQ=255)
                                                                                                                                 ggACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGACGGTCAGCAGCAGGTCGGCCACGCCGTTATCGCCCTCGGCCTTGAt         >  1:75639/1‑150 (MQ=255)
                                                                                                                                  gACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGACGGTCAGCAGCAGGTCGGCCACGCCGTTATCGCCCTCGGCCTTGATg        >  1:1724462/1‑150 (MQ=255)
                                                                                                                                     gAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGACGGTCAGCAGCAGGTCGGCCACGCCGTTATCGCCCTCGGCCTTGATgcgc     >  1:1154153/1‑150 (MQ=255)
                                                                                                                                     gAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGACGGTCAGCAGCAGGTCGGCCACGCCGTTATCGCCCTCGGCCTTGATgcgc     <  2:75639/150‑1 (MQ=255)
                                                                                                                                     ccggCGCAGTACCGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGACGGTCAGCAGCAGGTCGGCCACGCCGTTATCGCCCTCGGCCTTGATgcgc     >  2:1770607/3‑150 (MQ=255)
                                                                                                                                       ggCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGACGGTCAGCAGCAGGTCGGCCACGCCGTTATCGCCCTCGGCCTTGATGCGCTg   >  2:1688844/1‑150 (MQ=255)
                                                                                                                                        gcgcGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGACGGTCAGCAGCAGGTCGGCCACGCCGTTATCGCCCTCGGCCTTGATGCGCTGc  >  1:689117/1‑150 (MQ=255)
                                                                                                                                            |                                                                                                                                                 
CGCAGGCACAGACGGCCTTCCCACTGTTTGTCAGCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGACGGTCAGCAGCAGGTCGGCCACGCCGTTATCGCCCTCGGCCTTGATGCGCTGC  >  NC_002947/5928632‑5928916

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