Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,123,943 T→A 9.6% intergenic (‑138/+83) 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,9430TA 9.6% 188.0 / 7.7 63intergenic (‑138/+83)PP_mr13/gcvT‑IcspA regulatory region/aminomethyltransferase
Reads supporting (aligned to +/- strand):  ref base T (25/31);  new base A (2/4);  total (27/36)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCAAGGCCACCTTCT  >  NC_002947/1123794‑1124087
                                                                                                                                                     |                                                                                                                                                
tACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAAGACAGTGCAGTGATTCTACGCCAAGGCCCGAGCTAAAAACAAGCCCTTGGTCGTATTGAAAACGGGCGGGTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCgcct                                                                                                                                                  >  1:1097241/1‑150 (MQ=255)
 aCTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCgcctg                                                                                                                                                 <  1:697431/150‑1 (MQ=255)
  cTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCgcctgg                                                                                                                                                >  1:175835/1‑150 (MQ=255)
     gCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCgcctggcct                                                                                                                                             >  2:1189434/1‑150 (MQ=255)
            tCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCAGGCCAGGCGCCa                                                                                                                                      <  1:283178/150‑1 (MQ=255)
            tCACGAACTGTTGAGACCCATGCAAACAAAGGGCAGTGACTATACGCCAACTGTCGAGCAAAAAACAAGTCCTTGGTCGTATGGAAAACTTGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCAGGCCAGGCTCCa                                                                                                                                      <  2:635588/150‑1 (MQ=255)
             cACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCAt                                                                                                                                     <  1:862365/150‑1 (MQ=255)
             cACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCAt                                                                                                                                     <  2:1358432/150‑1 (MQ=255)
             cACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCAt                                                                                                                                     <  1:664102/150‑1 (MQ=255)
                gAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCAGGCCAGGCGCCAt                                                                                                                                     <  2:35117/147‑1 (MQ=255)
                gAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCAGGCCAGGCGCCAt                                                                                                                                     >  1:35117/1‑147 (MQ=255)
                gAACTGTTGCGACCCATGAAAACACAGTGAAGTGACTCTACGCAACCTCTCGAGCAAAAAACAGGCCTTTGGTCGTATGAAAAACGGGGGGTTTGCCGTAACCCGCCCGTTTTGTGGGAAACGGCAGGGCCCCGGGCCGGGCGCCAtttt                                                                                                                                  <  2:175835/150‑1 (MQ=255)
                     gTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCaaaa                                                                                                                             <  1:846046/150‑1 (MQ=255)
                     cttGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCaaaa                                                                                                                             <  2:167765/149‑1 (MQ=255)
                            cccATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTg                                                                                                                      <  2:1097241/150‑1 (MQ=255)
                            cccATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTg                                                                                                                      <  1:831190/150‑1 (MQ=255)
                             ccATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCaaaa                                                                                                                             >  2:905306/1‑142 (MQ=255)
                             ccATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCaaaa                                                                                                                             <  1:905306/142‑1 (MQ=255)
                             ccATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTgg                                                                                                                     >  1:359056/1‑150 (MQ=255)
                                 gCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAg                                                                                                                 >  1:1377759/1‑150 (MQ=255)
                                   aaaCGCCGTGCAGTGACTCTACGCCCACTTTCTGACGAATAACAAGCCCTTGGTCATATGGAAGACTGGCGGTTTGCCGATAACTGACCGATTTGTCGTAAACGGCATGGCGCTTGGCCTGGTGCCATGTTCAAAAGTTATTGGGCAGGt                                                                                                               >  2:171230/1‑150 (MQ=255)
                                   aaaCACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGt                                                                                                               >  2:171216/1‑150 (MQ=255)
                                   aaaCACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGt                                                                                                               >  1:84074/1‑150 (MQ=255)
                                    aaCACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAAAAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGtt                                                                                                              <  1:198794/150‑1 (MQ=255)
                                          gtgcagtgACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGaaaa                                                                                                        >  2:987254/1‑150 (MQ=255)
                                            gcagtgACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACg                                                                                                      >  2:1377966/1‑150 (MQ=255)
                                             cagtgACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCAGGCCAGGCGCCAt                                                                                                                                     >  1:583769/1‑118 (MQ=255)
                                             cagtgACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCAGGCCAGGCGCCAt                                                                                                                                     <  2:583769/118‑1 (MQ=255)
                                               gtgACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAat                                                                                                   >  1:1041171/1‑150 (MQ=255)
                                                 gACTCTACGCCAACTCTCGAGCGACCATAGAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAatat                                                                                                 >  1:37233/1‑150 (MQ=255)
                                                    tctACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTc                                                                                              >  2:855146/1‑150 (MQ=255)
                                                     ctACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCt                                                                                             <  2:1377759/150‑1 (MQ=255)
                                                      tACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTa                                                                                            <  1:855146/150‑1 (MQ=255)
                                                      tACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTa                                                                                            <  2:359056/150‑1 (MQ=255)
                                                       aCGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCCGGTTCGAAAACGAATATGTCTAt                                                                                           <  1:404570/150‑1 (MQ=255)
                                                            aaCTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTAtgctgc                                                                                      >  2:720263/1‑150 (MQ=255)
                                                               tctcGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCAATAACCGACCGGTTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCa                                                                                   >  1:858958/1‑150 (MQ=255)
                                                                        aaaaaCAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTca                                                                          <  1:1189434/150‑1 (MQ=255)
                                                                         aaaaCAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATCCACAGCCGACGGCTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTcac                                                                         <  2:37233/150‑1 (MQ=255)
                                                                         aaaaCAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAAc                                                                                                       >  1:471941/1‑120 (MQ=255)
                                                                         aaaaCAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAAc                                                                                                       <  2:471941/120‑1 (MQ=255)
                                                                          aaaCAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTcaca                                                                        >  1:780665/1‑150 (MQ=255)
                                                                             cAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAg                                                                     >  2:197200/1‑150 (MQ=255)
                                                                                 cccTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGt                                                                 <  1:720263/150‑1 (MQ=255)
                                                                                  ccTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGtt                                                                <  1:197200/150‑1 (MQ=255)
                                                                                    ttGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTca                                                                          <  2:1280353/138‑1 (MQ=255)
                                                                                    ttGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTca                                                                          >  1:1280353/1‑138 (MQ=255)
                                                                                              ggAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAat                                                                                                   <  1:304487/103‑1 (MQ=255)
                                                                                              ggAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAat                                                                                                   >  2:304487/1‑103 (MQ=255)
                                                                                                       gCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACgtgg                                           >  1:980703/1‑150 (MQ=255)
                                                                                                        cGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACgtggt                                          <  1:987254/150‑1 (MQ=255)
                                                                                                                cGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGaa                                  >  1:1293145/1‑150 (MQ=255)
                                                                                                                 gATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGaaa                                 <  2:84074/150‑1 (MQ=255)
                                                                                                                          ccgaTTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAGCGCTAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTg                        <  1:449724/150‑1 (MQ=255)
                                                                                                                            ggtttGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCt                      <  2:858958/148‑1 (MQ=255)
                                                                                                                             aTTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGcc                                                                                     >  1:167423/1‑86 (MQ=255)
                                                                                                                             aTTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGcc                                                                                     <  2:167423/86‑1 (MQ=255)
                                                                                                                              tttGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGa                    <  1:1377966/150‑1 (MQ=255)
                                                                                                                                 gTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACtt                 >  1:273784/1‑150 (MQ=255)
                                                                                                                                  tCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAACCGAATATGTCGATACTCACCAGCCCTTTGACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTGTGCTGACttt                >  1:678286/1‑150 (MQ=255)
                                                                                                                                  tCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTc                                                                                              <  1:816625/72‑1 (MQ=255)
                                                                                                                                  tCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTc                                                                                              >  2:816625/1‑72 (MQ=255)
                                                                                                                                     gAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCaa             <  2:1041171/150‑1 (MQ=255)
                                                                                                                                     gAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCaa             <  2:273784/150‑1 (MQ=255)
                                                                                                                                       aaCGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTTGTGTCACGAAAGGCATTTTGCTGACTTTCAAgg           <  1:1266935/150‑1 (MQ=255)
                                                                                                                                            cATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCAAGGCCAcc      >  2:603931/1‑150 (MQ=255)
                                                                                                                                                gcgcCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCAAGGCCACcttct  <  2:780665/150‑1 (MQ=255)
                                                                                                                                                  gcctggcctgGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCaatgc          >  1:687654/1‑137 (MQ=255)
                                                                                                                                                  gcctggcctgGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCAATGc          <  2:687654/140‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
TACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCAAGGCCACCTTCT  >  NC_002947/1123794‑1124087

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