Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,822,232 A→C 6.6% intergenic (+32/+26) PP_mr22 → / ← fdxA PrrB_RsmZ ncRNA/ferredoxin 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,822,2320AC 6.6% 244.5 / 8.0 91intergenic (+32/+26)PP_mr22/fdxAPrrB_RsmZ ncRNA/ferredoxin 1
Reads supporting (aligned to +/- strand):  ref base A (46/38);  new base C (3/3);  total (49/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.65e-01

ATGCACAGGCTTTCAGGATGAAGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTCCGCCAGCTCGGCGTTCAGCTCGATGAAG  >  NC_002947/1822084‑1822381
                                                                                                                                                     |                                                                                                                                                     
aTGCACAGGCTTTCAGGATGAAGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATTTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGAAAAACAAAAAAGGCCCGAAg                                                                                                                                                      >  1:2136881/1‑150 (MQ=255)
aTGCACAGGCTTTCAGGATGAAGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAg                                                                                                                                                      >  1:1984213/1‑150 (MQ=255)
    acaGGCTTTCAGGATGAAGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGcc                                                                                                                                                  <  1:1303118/150‑1 (MQ=255)
     caGGCTTTCAGGATGAAGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCt                                                                                                                                                 >  1:472820/1‑150 (MQ=255)
         cTTTCAGGATGAAGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCttttt                                                                                                                                             <  1:1091525/150‑1 (MQ=255)
                gctGAAGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTc                                                                                                                                            >  2:273193/3‑144 (MQ=255)
                gATGAAGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTc                                                                                                                                            <  1:273193/144‑1 (MQ=255)
                  tGAAGGCCAGTGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTc                                                                                                                                            >  1:1086810/1‑142 (MQ=255)
                  tGAAGGCCAGAGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTc                                                                                                                                            <  2:1086810/142‑1 (MQ=255)
                     aGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGGCTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGgcg                                                                                                                                 <  2:1126400/150‑1 (MQ=255)
                     aGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGgcg                                                                                                                                 >  1:1401925/1‑150 (MQ=255)
                      ggCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGgcg                                                                                                                                 >  1:202826/1‑150 (MQ=255)
                       gCCAGGGACACCGAAAGGAGCGACCTCATCAGGATGATGTTGGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTCCCCGTGGGAAAACAAAAGAGGCCCGAAGGTCCTTCTTCGCGTTTGGGCGCt                                                                                                                               >  2:386523/1‑150 (MQ=255)
                         cAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGCACAAGCAGGGACTATGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTAt                                                                                                                             <  2:1592796/150‑1 (MQ=255)
                            ggACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAg                                                                                                                          >  2:1171624/1‑150 (MQ=255)
                            ggACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAg                                                                                                                          <  2:570750/150‑1 (MQ=255)
                              aCATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAgcg                                                                                                                        >  1:976557/1‑150 (MQ=255)
                               cATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAgcgc                                                                                                                       <  1:2317046/150‑1 (MQ=255)
                                 tCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTc                                                                                                                     <  2:2107556/150‑1 (MQ=255)
                                 tCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTc                                                                                                                     >  1:675870/1‑150 (MQ=255)
                                            gATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCTTCGGGCCttttt                                                                                                                                             >  1:895667/1‑115 (MQ=255)
                                            gATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCTTCGGGCCttttt                                                                                                                                             <  2:895667/115‑1 (MQ=255)
                                            gATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCaa                                                                                                          <  2:868771/150‑1 (MQ=255)
                                            gATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCaa                                                                                                          <  1:2191057/150‑1 (MQ=255)
                                                tcatcaGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCtt                                                                                                      <  1:1181608/150‑1 (MQ=255)
                                                  atcaGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGc                                                                                                    >  1:1623218/1‑150 (MQ=255)
                                                  atcaGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCGGAGCGCGAGCAGCGCTCCAGGTCGGCAATCTTGc                                                                                                    >  2:1550433/1‑150 (MQ=255)
                                                      ggATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTAAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAgg                                                                                                <  2:768062/150‑1 (MQ=255)
                                                       gatgatGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGgc                                                                                                                                  >  1:179564/1‑115 (MQ=255)
                                                       gatgatGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGgc                                                                                                                                  <  2:179564/115‑1 (MQ=255)
                                                       gatgatGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGc                                                                                               >  1:1444550/1‑150 (MQ=255)
                                                        atgatgTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCt                                                                                              >  1:157208/1‑150 (MQ=255)
                                                         tgatgTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCtt                                                                                             >  2:2019159/1‑150 (MQ=255)
                                                         tgatgTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCtt                                                                                             <  2:1919557/150‑1 (MQ=255)
                                                           atgTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAAAAGGGCAGGGCTGAAGTGCCTGTTTCGCGTCTGGGCGCTATCAGCACTCCAGTTCGGCAATCTTGCCAGGCTTGc                                                                                           >  1:1796439/1‑150 (MQ=255)
                                                              tttGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCTTCGGGCCttttt                                                                                                                                             <  1:1822730/97‑1 (MQ=255)
                                                              tttGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCTTCGGGCCttttt                                                                                                                                             >  2:1822730/1‑97 (MQ=255)
                                                              tttGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCAt                                                                                        >  1:1464662/1‑150 (MQ=255)
                                                              tttGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCAt                                                                                        >  1:2134822/1‑150 (MQ=255)
                                                                  ggACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCa                                                                                    <  1:1771583/150‑1 (MQ=255)
                                                                  ggACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATATGTTCTACAGGTCGGCAATCTTGCCAGGCTTGCCATCCCa                                                                                    >  1:1289990/1‑150 (MQ=255)
                                                                    aCAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACNGCCCCTNTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCTCGGTCGGCAATCTTGCCAGGCTTGCCATCCCAtt                                                                                  >  1:1537654/1‑150 (MQ=255)
                                                                    aCAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCAtt                                                                                  <  2:675870/150‑1 (MQ=255)
                                                                    aCAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCAtt                                                                                  >  2:2030714/1‑150 (MQ=255)
                                                                     cAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCAttc                                                                                 <  1:2030714/150‑1 (MQ=255)
                                                                        gCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCAttcttc                                                                              >  2:1326254/1‑150 (MQ=255)
                                                                         cAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCGCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCa                                                                             <  1:386523/150‑1 (MQ=255)
                                                                           gggACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATATTGACAGGCTTGCCATCCACTTCTTAAGc                                                                           >  1:2232508/1‑150 (MQ=255)
                                                                               cTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTcc                                                                       <  2:157208/150‑1 (MQ=255)
                                                                                  cGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGc                                                                    >  2:341915/1‑150 (MQ=255)
                                                                                    gAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCg                                                                      >  1:2025261/1‑146 (MQ=255)
                                                                                    gAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCg                                                                      <  2:2025261/146‑1 (MQ=255)
                                                                                      aaaaaaaTGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAgg                                                                 <  2:472820/150‑1 (MQ=255)
                                                                                      aaaaaaaTGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAgg                                                                 <  1:2019159/150‑1 (MQ=255)
                                                                                       aaaaaaTGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCa                                                                                                                           >  1:189235/1‑91 (MQ=255)
                                                                                       aaaaaaTGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCa                                                                                                                           <  2:189235/91‑1 (MQ=255)
                                                                                        aaaaaTGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGc                                                               <  1:1326254/150‑1 (MQ=255)
                                                                                          aaaTGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGt                                                             >  1:170405/1‑150 (MQ=255)
                                                                                          aaaTGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGt                                                             >  2:400631/1‑150 (MQ=255)
                                                                                                 ggCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCCTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACg                                                      >  1:1579551/1‑150 (MQ=255)
                                                                                                      gTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCgg                                                 <  1:400631/150‑1 (MQ=255)
                                                                                                         aaaCCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGa                                              <  1:1550433/150‑1 (MQ=255)
                                                                                                          aaCCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGat                                             >  1:1262459/1‑150 (MQ=255)
                                                                                                           aCCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGata                                            <  2:170405/150‑1 (MQ=255)
                                                                                                            ccGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAg                                                                  <  2:420982/127‑1 (MQ=255)
                                                                                                            ccGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAg                                                                  >  1:420982/1‑127 (MQ=255)
                                                                                                               ccccTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTTGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAg                                        >  1:1904649/1‑150 (MQ=255)
                                                                                                               ccccTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCGGGCAGGGCGTCCTTACGCTCGGTGATATTAg                                        >  1:981185/1‑150 (MQ=255)
                                                                                                               ccccTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAg                                        >  1:1714828/1‑150 (MQ=255)
                                                                                                               ccccTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAg                                        >  2:352221/1‑150 (MQ=255)
                                                                                                                cccTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAgg                                       <  1:352221/150‑1 (MQ=255)
                                                                                                                  cTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGcc                                     <  2:191965/150‑1 (MQ=255)
                                                                                                                          gCCCGTAGGAAAACAAAAAAGGCCAGATGCTCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTcc                             <  1:1171624/150‑1 (MQ=255)
                                                                                                                            ccGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTCCGc                           <  2:1714828/150‑1 (MQ=255)
                                                                                                                            ccGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTCCGc                           <  2:1984213/150‑1 (MQ=255)
                                                                                                                            ccGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTCCGc                           <  2:1444550/150‑1 (MQ=255)
                                                                                                                             cGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTCCGcc                          >  2:1726849/1‑150 (MQ=255)
                                                                                                                                     aaCAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTCCGCCAGCTCGGc                  >  2:50418/1‑150 (MQ=255)
                                                                                                                                      gcgaaaaaGGCCCTTCGGGCCTTTTTTGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCtt                                                                                             <  2:2119770/71‑1 (MQ=25)
                                                                                                                                      gcgaaaaaGGCCCTTCGGGCCTTTTTTGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCtt                                                                                             >  1:2119770/4‑74 (MQ=25)
                                                                                                                                      aCAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTCCGCCAGCTCGGCg                 <  2:1425898/150‑1 (MQ=255)
                                                                                                                                        aaaaaaGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGt                                                >  1:1707461/1‑117 (MQ=255)
                                                                                                                                        aaaaaaGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGt                                                <  2:1707461/117‑1 (MQ=255)
                                                                                                                                        aaaaaaGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTCCGCCAGCTCGGCGtt               >  1:233667/1‑150 (MQ=255)
                                                                                                                                        aaaaaaGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTCCGCCAGCTCGGCGtt               >  2:487854/1‑150 (MQ=255)
                                                                                                                                                cccGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTCCGCCAGCTCGGCGTTCAGCTCGa       <  2:1518244/150‑1 (MQ=255)
                                                                                                                                                cccGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTCCGCCAGCTCGGCGTTCAGCTCGa       >  2:1898130/1‑150 (MQ=255)
                                                                                                                                                 ccGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTCAGCCAGATAGGAGATAAGCTCGAt      >  2:1297605/1‑150 (MQ=255)
                                                                                                                                                 ccGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTACTTACGCTCGGTGATATCAGGCCAGATTTCCg                            >  1:827172/1‑128 (MQ=255)
                                                                                                                                                 ccGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTACTTACGCTCGGTGATATCAGGCCAGATTTCCg                            <  2:827172/128‑1 (MQ=255)
                                                                                                                                                   gAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTCCGCCAGCTCGGCGTTCAGCTCGATGa    >  1:1223107/1‑150 (MQ=255)
                                                                                                                                                     aGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTCCGCCAGCTCGGCGTTCAGCTCGATGAAg  <  2:2136881/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATGCACAGGCTTTCAGGATGAAGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTCCGCCAGCTCGGCGTTCAGCTCGATGAAG  >  NC_002947/1822084‑1822381

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: