Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,123,948 T→A 6.9% intergenic (‑143/+78) PP_mr13 ← / ← gcvT‑I cspA regulatory region/aminomethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,123,9480TA 6.9% 260.2 / 8.1 87intergenic (‑143/+78)PP_mr13/gcvT‑IcspA regulatory region/aminomethyltransferase
Reads supporting (aligned to +/- strand):  ref base T (40/40);  new base A (3/3);  total (43/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

CAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCAAGGCCACCTTCTTTCCA  >  NC_002947/1123802‑1124092
                                                                                                                                                  |                                                                                                                                                
cAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCgcctggcctgac                                                                                                                                               >  2:621481/1‑148 (MQ=255)
    tCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCAGGCCAGGCGCCa                                                                                                                                           <  1:891282/150‑1 (MQ=255)
           cTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGACCGAAAAACAAGCGCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGACGCCATGTTCaa                                                                                                                                    >  2:1326783/1‑150 (MQ=255)
             gTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCaaaa                                                                                                                                  <  2:2111135/150‑1 (MQ=255)
             gTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCaaaa                                                                                                                                  <  2:2216214/150‑1 (MQ=255)
                   aCCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTAtt                                                                                                                            <  1:773125/150‑1 (MQ=255)
                    cccATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGATCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCGTCCCTGGCGCCATGTTCAAAAGTTATTg                                                                                                                           <  1:65731/150‑1 (MQ=255)
                     ccATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTgg                                                                                                                          >  2:431969/1‑150 (MQ=255)
                     ccATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTgg                                                                                                                          <  1:1571004/150‑1 (MQ=255)
                       aTGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGCCTGGCAGTCAGTTGGTCGAATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGc                                                                                                                        >  2:82617/1‑150 (MQ=255)
                       aTGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCAt                                                                                                                                          <  1:1768917/132‑1 (MQ=255)
                       aTGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCAt                                                                                                                                          >  2:1768917/1‑132 (MQ=255)
                        tGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCa                                                                                                                       >  1:2350449/1‑150 (MQ=255)
                        tGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGCCGCGGGGCCTGGCGCCATGTTCAAAAGTTATTGGGCa                                                                                                                       <  1:2019062/150‑1 (MQ=255)
                          cAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAgg                                                                                                                     >  2:465816/1‑150 (MQ=255)
                           aaaCACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGt                                                                                                                    >  2:1850119/1‑150 (MQ=255)
                                        tgACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAata                                                                                                       <  2:2350449/150‑1 (MQ=255)
                                         gACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAatat                                                                                                      >  2:261636/1‑150 (MQ=255)
                                            tctACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGaa                                                                                                         >  1:1788974/1‑144 (MQ=255)
                                            tctACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGaa                                                                                                         <  2:1788974/144‑1 (MQ=255)
                                                   cAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAACACGAATATGTCTAtgctg                                                                                            >  2:823739/1‑150 (MQ=255)
                                                      ctctCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGccc                                                                                         >  2:1847617/1‑150 (MQ=255)
                                                       tctcGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCa                                                                                        >  2:209850/1‑150 (MQ=255)
                                                       tctcGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCa                                                                                        >  2:528156/1‑150 (MQ=255)
                                                         tcGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCaaa                                                                                      >  1:2261457/1‑150 (MQ=255)
                                                          cGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCaaaa                                                                                     <  2:2025109/150‑1 (MQ=255)
                                                          cGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCaaaa                                                                                     <  1:817344/150‑1 (MQ=255)
                                                          cGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCaaaa                                                                                     <  1:1847617/150‑1 (MQ=255)
                                                           gAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAatat                                                                                                      <  1:1074184/132‑1 (MQ=255)
                                                           gAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAatat                                                                                                      >  2:1074184/1‑132 (MQ=255)
                                                           gAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAAc                                                                                    >  1:953560/1‑150 (MQ=255)
                                                             gCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCt                                                                                  >  2:1328752/1‑150 (MQ=255)
                                                              cGACAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTAGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGATAACCCCTCTGGCCATTTTGCCCAAAACCtt                                                                                 <  2:1403118/150‑1 (MQ=255)
                                                                aaaaaCAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTca                                                                               >  1:1831482/1‑150 (MQ=255)
                                                                aaaaaCAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTca                                                                               <  1:261636/150‑1 (MQ=255)
                                                                aaaaaCAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTca                                                                               >  1:1831860/1‑150 (MQ=255)
                                                                aaaaaCAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTca                                                                               <  1:209850/150‑1 (MQ=255)
                                                                    aCAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACaaa                                                                           >  2:370928/1‑150 (MQ=255)
                                                                    aCAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACaaa                                                                           <  1:528156/150‑1 (MQ=255)
                                                                     cAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAg                                                                          >  1:2165950/1‑150 (MQ=255)
                                                                      aaGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGc                                                                         <  1:1850119/150‑1 (MQ=255)
                                                                       aGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGcc                                                                        <  2:2129015/150‑1 (MQ=255)
                                                                          ccTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGtt                                                                     <  2:1831860/150‑1 (MQ=255)
                                                                          ccTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGtt                                                                     <  2:1831482/150‑1 (MQ=255)
                                                                           cTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGc                                                                         <  2:199560/145‑1 (MQ=255)
                                                                           cTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGc                                                                         >  1:199560/1‑145 (MQ=255)
                                                                           cTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGttt                                                                    >  2:2372629/1‑150 (MQ=255)
                                                                                tCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGACATGGCGCCAGGCCAGGCGCCATgcc                                                                                                                                       <  1:1372529/78‑3 (MQ=37)
                                                                                tCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGACATGGCGCCAGGCCAGGCGCCATgcc                                                                                                                                       >  2:1372529/1‑76 (MQ=38)
                                                                                       gAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACaaa                                                                           >  1:177488/1‑131 (MQ=255)
                                                                                       gAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACaaa                                                                           <  2:177488/131‑1 (MQ=255)
                                                                                       gAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGgtag                                                        >  2:2286601/1‑150 (MQ=255)
                                                                                        aaaaCTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGgtagt                                                       >  2:2303049/1‑150 (MQ=255)
                                                                                          aaCTGGCGGTGTGCCGATAACCGGCCGGTTTGTCGGAAACGGCATGGCGGCTGGCCTGGGGCCATGTTCAAAAGCTATTGGGCAGGTTCGAAAACGGATATGTCTATACTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGgtagtca                                                     >  2:374557/1‑148 (MQ=255)
                                                                                             tGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTa                                                                                                 >  2:1281998/1‑103 (MQ=255)
                                                                                             tGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTa                                                                                                 <  1:1281998/103‑1 (MQ=255)
                                                                                              ggCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGTCCTGGCGCCATGTTCAAAAGTTATTGGGCATGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACgtg                                                 <  2:1231962/150‑1 (MQ=255)
                                                                                                cGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACgtgga                                               <  1:2303049/150‑2 (MQ=255)
                                                                                                 ggTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGgta                                                         <  2:1202772/139‑1 (MQ=255)
                                                                                                 ggTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGgta                                                         >  1:1202772/1‑139 (MQ=255)
                                                                                                  gTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGgtgt                                             >  2:253700/1‑150 (MQ=255)
                                                                                                         gATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTAt                                                                                                <  2:1290878/92‑1 (MQ=255)
                                                                                                         gATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTAt                                                                                                >  1:1290878/1‑92 (MQ=255)
                                                                                                         gATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGaaa                                      <  2:2165950/150‑1 (MQ=255)
                                                                                                          aTAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAg                                     <  1:2286601/150‑1 (MQ=255)
                                                                                                            aaCCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGc                                   >  1:1999804/1‑150 (MQ=255)
                                                                                                                 accgaTTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCAtttt                              >  1:2136983/1‑150 (MQ=255)
                                                                                                                 accgaTTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCAtttt                              >  1:588563/1‑150 (MQ=255)
                                                                                                                    gaTTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCt                           <  1:2037701/150‑1 (MQ=255)
                                                                                                                    gaTTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCt                           <  2:2113270/150‑1 (MQ=255)
                                                                                                                    gaTTTGTCGGAAACGGAATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCt                           <  1:2063661/150‑1 (MQ=255)
                                                                                                                      tttGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGt                                                                                                    >  2:236348/1‑75 (MQ=255)
                                                                                                                      tttGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGt                                                                                                    <  1:236348/75‑1 (MQ=255)
                                                                                                                      tttGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGa                         <  1:1328752/150‑1 (MQ=255)
                                                                                                                       ttGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGAc                        <  1:1326783/150‑1 (MQ=255)
                                                                                                                           cGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTc                    <  1:253700/150‑1 (MQ=255)
                                                                                                                            ggAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAAc                                                    >  2:2216331/1‑117 (MQ=255)
                                                                                                                            ggAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAAc                                                    <  1:2216331/117‑1 (MQ=255)
                                                                                                                             gAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCACAAATTCTGCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGCAAGGCATTTTGCTGACTTTCaa                  >  2:791039/1‑150 (MQ=255)
                                                                                                                              aaaCGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCAAg                 <  1:370928/150‑1 (MQ=255)
                                                                                                                                aCGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCAAGGc               <  1:465816/150‑1 (MQ=255)
                                                                                                                                  ggCATGGCGCCTGGCCGGGCGCCATGTTCAAAAGTTATTGGGCATGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCAAGGCCa             <  1:431969/150‑1 (MQ=255)
                                                                                                                                   gCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCAAGGCCAc            >  1:2310189/1‑150 (MQ=255)
                                                                                                                                   acATGGCGCCAGGCCAGGCGCCATGCTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCAAGGCCAc            <  1:402417/149‑1 (MQ=255)
                                                                                                                                    cATGGCGCCAGGCCAGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCAAGGCCAcc           >  1:132559/1‑150 (MQ=255)
                                                                                                                                     aTGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCAAGGCCACct          >  2:696604/1‑150 (MQ=255)
                                                                                                                                     aTGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCAAGGCCACct          >  2:842215/1‑150 (MQ=255)
                                                                                                                                     aTGGCGCCAGGCCAGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCAAGGCCACct          >  1:1158749/1‑150 (MQ=255)
                                                                                                                                          gccaggccaggcGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGaaa                                      <  2:1168193/113‑1 (MQ=255)
                                                                                                                                          gccaggccagGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGaaa                                      >  1:1168193/1‑117 (MQ=255)
                                                                                                                                            caggccaggcGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCAAGGCCACCTTCTTTcc   <  2:1158749/148‑1 (MQ=255)
                                                                                                                                             tggcctggcGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCAAGGCCACCTTCTTTCCa  >  2:308797/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
CAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCAAGGCCACCTTCTTTCCA  >  NC_002947/1123802‑1124092

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