Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,077,285 A→C 10.6% intergenic (+32/+76) PP_4470 → / ← mgtE Arc domain‑containing transcriptional regulator/magnesium transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,077,2850AC10.6% 277.9 / 23.2 110intergenic (+32/+76)PP_4470/mgtEArc domain‑containing transcriptional regulator/magnesium transporter
Reads supporting (aligned to +/- strand):  ref base A (50/48);  new base C (5/7);  total (55/55)
Fisher's exact test for biased strand distribution p-value = 7.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 110 total aligned reads displayed.

TGTCCCTGCATGAGCGCGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAACCACCGCTGTCGGTCACAGCAGTGATCA  >  NC_002947/5077138‑5077431
                                                                                                                                                   |                                                                                                                                                  
tGTCCCTGCATGAGCGCGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGc                                                                                                                                                  <  1:1986199/150‑1 (MQ=255)
          tGAGCGCGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACAAAAAAGCCAGCTTACGCTGGCtttttc                                                                                                                                        <  2:2055284/150‑2 (MQ=255)
            aGCGCGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTtgt                                                                                                                                      <  2:940571/150‑1 (MQ=255)
             gcgcGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTtgtg                                                                                                                                     >  1:2187171/1‑150 (MQ=255)
              cgcgAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTgg                                                                                                                                    <  2:1056253/150‑1 (MQ=255)
              cgcgAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTgg                                                                                                                                    <  2:1642538/150‑1 (MQ=255)
              cgcgAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTgg                                                                                                                                    >  1:1300643/1‑150 (MQ=255)
                cgAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGc                                                                                                                                  <  2:258245/150‑1 (MQ=255)
                         gcAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTgg                                                                                                                        <  1:223323/150‑1 (MQ=255)
                          cAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTgg                                                                                                                        >  1:31849/1‑150 (MQ=255)
                          cAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTgg                                                                                                                        >  1:871816/1‑150 (MQ=255)
                          cAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTgg                                                                                                                        <  1:433426/150‑1 (MQ=255)
                             cGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGgcag                                                                                                                     >  2:1269158/1‑150 (MQ=255)
                              gTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGgcagc                                                                                                                    <  2:871816/150‑1 (MQ=255)
                              gTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGgcagc                                                                                                                    >  2:1946155/1‑150 (MQ=255)
                               tttCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAAAAGCGATAAAAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGgcagca                                                                                                                   >  1:1570704/1‑150 (MQ=255)
                                ttCCTCCAACTGCCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCaa                                                                                                                  <  2:492041/150‑1 (MQ=255)
                                    gCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAACAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGcc                                                                                                              >  2:883351/1‑150 (MQ=255)
                                    gCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGcc                                                                                                              >  2:1384419/1‑150 (MQ=255)
                                    gCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTAGCTTTTTTGTGGGCGTCTGACTGCCAGCAAGccc                                                                                                              >  2:492080/1‑150 (MQ=255)
                                      cAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCCGCGTACGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCAt                                                                                                            <  1:450583/150‑1 (MQ=255)
                                          tGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTcc                                                                                                        >  1:54084/1‑150 (MQ=255)
                                            tCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCt                                                                                                      >  1:1270681/1‑150 (MQ=255)
                                              ccACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCTTACGCTGGCttttt                                                                                                                                         <  1:705157/113‑1 (MQ=255)
                                              ccACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCTTACGCTGGCttttt                                                                                                                                         >  2:705157/1‑113 (MQ=255)
                                              ccACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGATGGATTTTCTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGc                                                                                                    >  1:2269951/1‑150 (MQ=255)
                                               cACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGc                                                                                                               <  2:1938363/138‑1 (MQ=255)
                                               cACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGc                                                                                                               >  1:1938363/1‑138 (MQ=255)
                                               cACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCg                                                                                                   <  1:883351/150‑1 (MQ=255)
                                                aCCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCgg                                                                                                  <  2:2187171/150‑1 (MQ=255)
                                                aCCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCgg                                                                                                  >  1:204038/1‑150 (MQ=255)
                                                 ccGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCggg                                                                                                 >  2:646809/1‑150 (MQ=255)
                                                 ccGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGACATTTCCCTGCggg                                                                                                 >  1:1060625/1‑150 (MQ=255)
                                                  cGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGa                                                                                                <  1:673206/150‑1 (MQ=255)
                                                    ccAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGATCTCGACGGCGTCCCCCCCTTAGAGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGAtt                                                                                              <  1:861509/150‑1 (MQ=255)
                                                    ccAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCt                                                                                                                                                 >  2:2296964/1‑99 (MQ=255)
                                                    ccAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCt                                                                                                                                                 <  1:2296964/99‑1 (MQ=255)
                                                        cagAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACt                                                                                                                          >  2:1531774/1‑118 (MQ=255)
                                                        cagAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACt                                                                                                                          <  1:1531774/118‑1 (MQ=255)
                                                        cacaaCGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTc                                                                                          >  1:2333222/4‑150 (MQ=255)
                                                              gCCTTGGTCGCCCTCATCGCTCACGCTGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAg                                                                                                                 >  1:1531664/1‑121 (MQ=255)
                                                              gCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAg                                                                                                                 <  2:1531664/121‑1 (MQ=255)
                                                              gCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAg                                                                                    <  1:468393/150‑1 (MQ=255)
                                                              gCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAg                                                                                    <  2:54084/150‑1 (MQ=255)
                                                                    gTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCACGGCCCGTTCCCTGCGGGATTGTTCCCACAGggaggg                                                                              >  1:622186/1‑149 (MQ=255)
                                                                      cGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGAtt                                                                            <  1:221901/150‑1 (MQ=255)
                                                                      cGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGAtt                                                                            <  1:1726962/150‑1 (MQ=255)
                                                                       gCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGAtta                                                                           <  1:2097753/150‑2 (MQ=255)
                                                                        cccTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTgg                                                                          <  2:1648372/150‑1 (MQ=255)
                                                                         ccTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGCTTGGc                                                                         >  2:1000654/1‑150 (MQ=255)
                                                                         ccTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGc                                                                         <  1:646809/150‑1 (MQ=255)
                                                                             aTCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTac                                                                     >  2:1922447/1‑150 (MQ=255)
                                                                                gCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGc                                                                                                                                  >  2:1266456/1‑86 (MQ=255)
                                                                                gCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGc                                                                                                                                  <  1:1266456/86‑1 (MQ=255)
                                                                                 cTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTacaaca                                                                 >  1:702964/1‑150 (MQ=255)
                                                                                  tCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCtttt                                                                                                                                          <  2:2226127/76‑1 (MQ=255)
                                                                                  tCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCtttt                                                                                                                                          >  1:2226127/1‑76 (MQ=255)
                                                                                  tCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAg                                                                >  2:487307/1‑150 (MQ=255)
                                                                                            gAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTg                                                      >  2:1292531/1‑150 (MQ=255)
                                                                                            gAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTg                                                      <  1:1269158/150‑1 (MQ=255)
                                                                                               aTGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGcc                                                   <  1:1384419/150‑1 (MQ=255)
                                                                                                 ggCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGAt                                                                                               <  1:938269/104‑1 (MQ=255)
                                                                                                 ggCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGAt                                                                                               >  2:938269/1‑104 (MQ=255)
                                                                                                  gCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGt                                                       >  2:1459804/1‑143 (MQ=255)
                                                                                                  gCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGt                                                       <  1:1459804/143‑1 (MQ=255)
                                                                                                      cAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCaa                                            <  1:1421230/150‑1 (MQ=255)
                                                                                                           gCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAg                                       >  2:580750/1‑150 (MQ=255)
                                                                                                             ctcctGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTtgtg                                                                                                                                     >  1:2059348/1‑54 (MQ=255)
                                                                                                             ctcctGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTtgtg                                                                                                                                     <  2:2059348/54‑1 (MQ=255)
                                                                                                              tcctGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATg                                    <  2:2269951/150‑1 (MQ=255)
                                                                                                                ctGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATgaa                                  <  2:204038/150‑1 (MQ=255)
                                                                                                                 tGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTATGACAGGGAGCAAGGGTATTTGCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCATCGCAAGGAAGATgaag                                 >  1:1450510/1‑150 (MQ=255)
                                                                                                                       cTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAacca                           <  2:2333222/150‑1 (MQ=255)
                                                                                                                         gCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCATCCTAAGGAAGATGAAGAAaccacc                         <  2:873392/150‑1 (MQ=255)
                                                                                                                         gCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAaccacc                         <  2:702964/150‑1 (MQ=255)
                                                                                                                          cAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAACCACCg                        >  1:306113/1‑150 (MQ=255)
                                                                                                                            gCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGAt                                     >  2:110822/1‑135 (MQ=255)
                                                                                                                            gCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGAt                                     <  1:110822/135‑1 (MQ=255)
                                                                                                                            gCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAaccaccact                      >  2:899490/1‑147 (MQ=255)
                                                                                                                               aaCACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAACCACCGCTGTc                   <  1:1145710/150‑1 (MQ=255)
                                                                                                                                 cacaAAAAAGCCAGCTTACGCTGGCTTTTTCGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCaa                                            >  1:35664/1‑123 (MQ=39)
                                                                                                                                 cacaAAAAAGCCAGCTTACGCTGGCTTTTTCGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCaa                                            <  2:35664/123‑1 (MQ=39)
                                                                                                                                 cacGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCTTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAg                                           <  2:948449/124‑1 (MQ=255)
                                                                                                                                 cacGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAg                                           >  1:948449/1‑124 (MQ=255)
                                                                                                                                  acaaaaaaGCCAGCTTACGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGAtt                                                                            >  2:1406741/4‑90 (MQ=38)
                                                                                                                                  acaaaaaaGCCAGCTTACGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGAtt                                                                            <  1:1406741/87‑1 (MQ=38)
                                                                                                                                   cGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAACCACCGCTGTCGGTc               >  2:495691/1‑150 (MQ=255)
                                                                                                                                      aaaaGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGgga                                                                                 <  2:1884216/81‑1 (MQ=255)
                                                                                                                                      aaaaGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGgga                                                                                 >  1:1884216/1‑81 (MQ=255)
                                                                                                                                       aaaGCCAGCTTACGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAg                                                                                    <  2:431844/77‑1 (MQ=39)
                                                                                                                                       aaaGCCAGCTTACGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAg                                                                                    >  1:431844/1‑77 (MQ=39)
                                                                                                                                        aaGCCAGCTTACGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATgaag                                 <  2:1922809/127‑1 (MQ=255)
                                                                                                                                        aaGCCAGCTTACGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATgaag                                 >  1:1922809/1‑127 (MQ=255)
                                                                                                                                          gCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTgg                                                     <  2:2152581/105‑1 (MQ=255)
                                                                                                                                          gCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTgg                                                     >  1:2152581/1‑105 (MQ=255)
                                                                                                                                           ccAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGAt                                     >  2:351734/1‑120 (MQ=255)
                                                                                                                                           ccAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGAt                                     <  1:351734/120‑1 (MQ=255)
                                                                                                                                             aGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAACCACCGCTGTCGGTCACAGCAGTGa     >  2:982849/1‑150 (MQ=255)
                                                                                                                                               cGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAACCACCGCTGTCGGTCACAGCAGTGatc   <  2:1326063/150‑1 (MQ=255)
                                                                                                                                                gTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGGGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAACCACCGCTGTCGGTCACAGCAGTGatca  >  1:1118723/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TGTCCCTGCATGAGCGCGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAACCACCGCTGTCGGTCACAGCAGTGATCA  >  NC_002947/5077138‑5077431

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