Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,822,230 G→T 6.7% intergenic (+30/+28) 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,2300GT 6.7% 275.4 / 6.1 90intergenic (+30/+28)PP_mr22/fdxAPrrB_RsmZ ncRNA/ferredoxin 1
Reads supporting (aligned to +/- strand):  ref base G (45/39);  new base T (3/3);  total (48/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 = 6.66e-01

ATGCACAGGCTTTCAGGATGAAGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTCCGCCAGCTCGGCGTTCAGCTCGATGA  >  NC_002947/1822084‑1822379
                                                                                                                                                   |                                                                                                                                                     
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                                                                                                        <  1:2191057/150‑1 (MQ=255)
                                            gATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCaa                                                                                                        <  2:868771/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:2134822/1‑150 (MQ=255)
                                                              tttGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCAt                                                                                      >  1:1464662/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                                                                    <  2:2025261/146‑1 (MQ=255)
                                                                                    gAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCg                                                                    >  1:2025261/1‑146 (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                                                           >  2:400631/1‑150 (MQ=255)
                                                                                          aaaTGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGt                                                           >  1:170405/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                                                                                           >  1:2119770/4‑74 (MQ=25)
                                                                                                                                      gcgaaaaaGGCCCTTCGGGCCTTTTTTGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCtt                                                                                           <  2:2119770/71‑1 (MQ=25)
                                                                                                                                      aCAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTCCGCCAGCTCGGCg               <  2:1425898/150‑1 (MQ=255)
                                                                                                                                        aaaaaaGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGt                                              <  2:1707461/117‑1 (MQ=255)
                                                                                                                                        aaaaaaGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGt                                              >  1:1707461/1‑117 (MQ=255)
                                                                                                                                        aaaaaaGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTCCGCCAGCTCGGCGtt             >  2:487854/1‑150 (MQ=255)
                                                                                                                                        aaaaaaGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTCCGCCAGCTCGGCGtt             >  1:233667/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)
                                                                                                                                                   |                                                                                                                                                     
ATGCACAGGCTTTCAGGATGAAGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTCCGCCAGCTCGGCGTTCAGCTCGATGA  >  NC_002947/1822084‑1822379

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