Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 941,229 G→T G91C (GGC→TGC)  dmsA → dimethyl sulfoxide reductase, anaerobic, subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913941,2290GT95.7% 199.0 / NA 69G91C (GGC→TGC) dmsAdimethyl sulfoxide reductase, anaerobic, subunit A
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base T (40/26);  total (41/28)
Fisher's exact test for biased strand distribution p-value = 5.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGCGCACGCTGTCGATAGCGCCATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCGGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGACCGCCTGAAATATCCGATGAAACGAGTCGGGGCGCGCGGTGAAGGCAAATTCGAGCGCATTA  >  NC_000913/941083‑941373
                                                                                                                                                  |                                                                                                                                                
ttGCGCACGCTGTCGATAGCGCCATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCTGc                                                                                                                                                >  2:1145460/1‑149 (MQ=255)
   cgcACGCTGTCGATAGCGCCATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCTGCGac                                                                                                                                             >  2:490055/1‑148 (MQ=255)
   cgcACGCGGTCGATAGCGCCATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGAAGGTGAAATCAAATCTGACGACACGGACAATACCTGCGAt                                                                                                                                             <  2:125547/149‑1 (MQ=255)
    gcACGCTGTCGATAGCGCCCTGCCAACAAACTCAGACGGAAAGGTTATCTGGAGCGCCTGTACAGTTAACTGTGGTAGTCCCTGCCCGCTACGTCTGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCTGCGATg                                                                                                                                            <  2:854679/149‑1 (MQ=255)
      aCGCTGTCGATAGCGCCATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCTGCGATGAc                                                                                                                                          <  2:876703/149‑1 (MQ=255)
           ggcGATAGCGCCATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCTGCGATGACAATTa                                                                                                                                     <  2:1115950/147‑1 (MQ=255)
                  gcgcCATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGGCAATACCTGCGGTGGCAATTACGACGGc                                                                                                                              >  1:361294/1‑149 (MQ=255)
                  gcgcCATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCTGCGATGACAATTACGACGGc                                                                                                                              >  2:945918/1‑149 (MQ=255)
                    gcCATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACGATACCTGCGATGACAATTAAGACGGCCt                                                                                                                            >  1:1321127/1‑149 (MQ=255)
                       aTTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCTGCGATGACAATTACGACGGCCTGCa                                                                                                                         >  1:1153964/1‑149 (MQ=255)
                           cacaaaaaGCAGACGAAAGTGTTATCTTGACCGCCTGTACAGGTAACTGTGGTAGTCGCTGCCCGTTATGTATGCTCGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCTGCGATGACAATTACGACGGCCTGCACCAg                                                                                                                     <  2:153058/145‑1 (MQ=255)
                           caacaaAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCTGCGATGACAATTACGACGGCCTGCACCAg                                                                                                                     >  2:980519/1‑149 (MQ=255)
                            aacaaAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACTGTGGTAGTCGCTGCCCGCTACGTCTGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCTGCGATGACAATTACGACGGCCTGCACCAgg                                                                                                                    >  2:1278132/1‑149 (MQ=255)
                                     aGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCTGCGATGACAATTACGACGGCCTGCACCAGGTTCGCgcct                                                                                                           <  2:1338093/149‑1 (MQ=255)
                                     aGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCTGCGATGACAATTACGACGGCCTGCACCAGGTTCGCgcct                                                                                                           <  1:460162/149‑1 (MQ=255)
                                      gACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACATGCGATGACAATTACGACGGCCTGCACCAGGTTCGCggctg                                                                                                          >  2:364156/1‑149 (MQ=255)
                                       aCGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAACTACCTCCTGACGTACAATACCTGCGATGACAATTACGACGGCCTGCACCAGGTTCGCgcctgc                                                                                                         <  1:1059440/149‑1 (MQ=255)
                                         gAAAAGGTTATCTGGAGCGCCTGTACAGTTAACTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATCCCTGCGATGACAATTACGACGGACTGCACCAGGTCCGCgcctgcct                                                                                                       >  1:1141527/1‑149 (MQ=255)
                                         gAAAAGGTTATCTGGAGCGCCTGTACAGTTAACTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCTGCGATGACAATTACGACGGCCTGCACCAGGTTCGCgcctgcct                                                                                                       >  2:851579/1‑149 (MQ=255)
                                         gAAAAGGTTATCTGGAGCGCCTGTACAGTTAACTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCTGCGATGACAATTACGACGGCCTGCACCAGGTTCGCgcctgcct                                                                                                       >  2:835417/1‑149 (MQ=255)
                                            aaGGTTATCTGGAGCGCCTGTACAGTTAACTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGAGAATACCGGCGAAGACAATTACGACGGCCTGCACCAGGTTCCCCGCTGCCcgcg                                                                                                    >  2:886462/1‑149 (MQ=255)
                                             aGGTTATCTGGAGCGCCTGTACAGTTAACTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCTGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGt                                                                                                   <  2:1345620/149‑1 (MQ=255)
                                                    cgggggcgcCCGGACAGTTAACGGTGGTAGTCGCTGCCCGCCACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCTGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGtt                                                                                            <  2:712274/144‑1 (MQ=255)
                                                     tGGTGCGCCTGTACCGTTAACTGTGGTAGTCGCTGCCCGCTACGTATCCAGTGCTTGTGCTGTGATATCAAATATTTCGAAACGGACAATACCTGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTc                                                                                           <  1:597477/149‑1 (MQ=255)
                                                     tGGAGCGCCTGTACAGTTAACTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCTGCGATGGCAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTc                                                                                           >  1:819989/1‑149 (MQ=255)
                                                             cTGTACAGTTAACTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCTGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTc                                                                                   >  1:964458/1‑149 (MQ=255)
                                                                  cAGTTAACTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCTGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCgt                                                                              <  2:1153964/149‑1 (MQ=255)
                                                                   aGTTAACTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCTGCGATGCCAATTACGACGGCCTGCACCGGGTTCGCGCCGGCCTGCGTGGGCGCTCCTTGCGTCGCCgtg                                                                             >  1:364595/1‑149 (MQ=255)
                                                                   aGTTAACTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCTGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCgtg                                                                             >  1:410604/1‑149 (MQ=255)
                                                                    gTTATCTGTGGTAGTCGCTGCCCGCTCCGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGTACAACACCTGCGATGACAATTACGACGGCCTGCACCAGGGTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCgtgt                                                                            <  1:366729/149‑1 (MQ=255)
                                                                    gTTAACTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCTGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCgtgt                                                                            >  2:964197/1‑149 (MQ=255)
                                                                      tAACTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCTGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCt                                                                          >  2:1240768/1‑149 (MQ=255)
                                                                       aaCTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCTGCGATGACAATTACGACGGCCTGCACCAGGTTCGCCCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGtata                                                                         >  1:315699/1‑149 (MQ=255)
                                                                          tgtgGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCTGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCGGCCTGGGCGTTCCCCGCGCCGCCGTGTCCAcca                                                                      >  1:536556/1‑147 (MQ=255)
                                                                              gTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGGGAAATCAAATATGTCGGAACGGACAATACCTGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCg                                                                  >  1:1303905/1‑149 (MQ=255)
                                                                                aGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCTGGGATGACAATTACGACGGCATGCACCAGGTTCGCGCCTGCATGAGTGGGCGTTCCATGCGCCGCCGTGTCTAAAGTCCGGa                                                                >  1:1292014/1‑149 (MQ=255)
                                                                                 gTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCTGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTCCCATGCGTCGCCGCGTCCACAAACCGGAc                                                               >  2:681193/1‑149 (MQ=255)
                                                                                   cGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCTGCGATGGCAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGGTCCATGCGTCGCCGTGTCTACAATCCGGCCCg                                                             >  2:251037/1‑149 (MQ=255)
                                                                                   cGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCTGCGATGACAATTACGACGGCCTGCACCAGGCTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGACCg                                                             >  2:590513/1‑149 (MQ=255)
                                                                                    gCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCTGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCCACAATCCGGAcgcc                                                            >  2:651826/1‑146 (MQ=255)
                                                                                      tGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATTTGTCGAAGGGGACGATACCTGCGATGATCCTTTCGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGACCGCCt                                                          <  1:261202/149‑1 (MQ=255)
                                                                                      tGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCTGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCCTGCGGCGCCGTGTCTACAGTGCGGAACGcgc                                                          >  2:29811/1‑147 (MQ=255)
                                                                                         ccGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCTGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGACCGCCTGaa                                                       >  2:600031/1‑149 (MQ=255)
                                                                                          cGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCTGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGACCGCCTGaaa                                                      >  1:1302029/1‑149 (MQ=255)
                                                                                          cGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCTGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTCTCTACAATCCGGACCACCTGaaa                                                      >  1:1222462/1‑149 (MQ=255)
                                                                                            cTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCTTCGGTGGCAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGGCCGCCCGAAata                                                    >  2:1022249/1‑149 (MQ=255)
                                                                                            cTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCTGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGGGTCTACAATCCGGACCGCATGAAata                                                    >  2:356805/1‑149 (MQ=255)
                                                                                            cTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCTGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGACTGCGAGGGCGTTCCATGCGTCGCCGTGTCGGCAATCCGGAACGGCTGaaaaa                                                    >  2:886902/1‑147 (MQ=255)
                                                                                              aCGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCTGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGACCGCCTGAAATATc                                                  <  1:337951/149‑1 (MQ=255)
                                                                                                 tATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCTGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGACCGCCTGAAATATCCGa                                               >  1:945100/1‑149 (MQ=255)
                                                                                                   tggcggtcgtGGACGGGGAAATCAAGTAGGCCGAAACGGACAATACCGGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCGGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGACCGCCTGAAATATCCGATg                                             <  2:819739/144‑1 (MQ=255)
                                                                                                      aCGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCTGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGACCGCCTGAAATATCCGATGaaa                                          <  1:600031/149‑1 (MQ=255)
                                                                                                             ggACGGTGAAATCAAATATGTCGAAACGGACAATACCTGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGCCGCCGTGTCTACAATCCGGACCGCCTGAAATATCCGATGAAACGAGTCg                                   <  1:1278132/149‑1 (MQ=255)
                                                                                                                 ggggAAATCAAATATGTCGAAACGGACAATACCTGCGCTGACAATTACGCCGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAACCCGGCCCGCCTGAAATATCCGATGAAACGAGTCGGGgc                               <  2:41812/146‑1 (MQ=255)
                                                                                                                   tGAAATCAAATATGTCGAAACGGACAATACCTGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGGTCCATGCGTCGCCGTGTCTACAATCCGGACCGCCTGAAATATCCGAGGAAACGCGTCGGGgcgc                             >  1:450213/1‑149 (MQ=255)
                                                                                                                     aaaTCAAATATGTCGAAACGGACAATACCTGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGACCGCCTGAAATAGCCGATGAAACGAGTCGGGgcgcgc                           >  1:203501/1‑149 (MQ=255)
                                                                                                                        tCAAATATGTCGAAACGGACAATACCTGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGACCGCCTGAAATATCCGATGAAACGAGTCGGGGCGCGCGGt                        <  1:1240768/149‑1 (MQ=255)
                                                                                                                           aaTATGTCGAAACGGACAATACCTGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGACCGCCTGAAATATCCGATGAAACGAGTCGGGGCGCGCGGTGaa                     <  1:945918/149‑1 (MQ=255)
                                                                                                                           aaTATGCCGAAACGGACAATCCCTGCGATGACGATTACGACGGCCTGCACCAGGTCCGCGCCGGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGACCGCCTGAAATATCCGATGAAACGAGTCGGGGCGCGCGGTGaa                     <  2:315338/149‑1 (MQ=255)
                                                                                                                            agcggTCGAAACGGACAATACCGGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCTCCGTGTCTACAATCCGGACCGCCTGAAATATCCGATGAAACGAGTCGGGGCGCGCGGTGAAg                    <  1:356805/145‑1 (MQ=255)
                                                                                                                                 tCGAAACGGACAATACCTGCGATGACAATTACGATGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGACCGCCTGAAATATCCGATGAAACGAGTCGGGGCGCGCGGTGAAGGCaga               >  1:687196/1‑147 (MQ=255)
                                                                                                                                    aaagcgACAATACCTGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCAATGCGTCGCCGTGTCTACAATCCGGACCGCCTGAAATATCCGATGAAACGAGTCGGGGCGCGCGGTGAAGGCAAATTc            <  2:1187830/144‑1 (MQ=255)
                                                                                                                                      aCGGACAATACCTGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTCTCTACAATCCGGACCACCTGAAATATCCGATGAAACGAGTCGGGGCGCGCGGTGAAGGCAAATTCGa          <  2:1222462/149‑1 (MQ=255)
                                                                                                                                      aCGGACAATACCTGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCATGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGCCCGCCTGAAATATCCGATGAAACGAGTCGGGGGGCGCGCTGAAGGCAAATTCgg          >  2:185520/1‑148 (MQ=255)
                                                                                                                                      aCGGACAAGACCTGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGACCGCCTGAAATATCCGATGAACCGCGTCGGGGCGCGCGGTGAAGGCAAATTCGa          <  1:273407/149‑1 (MQ=255)
                                                                                                                                        ggACAATACCTGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGGGTCTACAATCCGGACCGCCTGGAATATCCGATGAAACGAGTCGGGGCGCGCGGTGAAGGCAAATTCGAgc        >  2:1199251/1‑149 (MQ=255)
                                                                                                                                           cAATACCTGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGACCGCCTGAAATATCCGATGAAACGAGCCGGGGCGCGCGGCGAAGGCAAATTCGAGCGCa     >  1:1236084/1‑149 (MQ=255)
                                                                                                                                            aaTACCTGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGACCGCCTGAAATATCCGATGAAACGATTCGGGGCGCGCGGTGAAGGCAAATTCGAGCGCAt    <  2:964458/149‑1 (MQ=255)
                                                                                                                                              tACCTGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGACCGCCTGAAATATCCGATGAAACGAGTCGGGGCGCGCGGTGAAGGCAAATTCGAGCGCATTa  <  1:851579/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
TTGCGCACGCTGTCGATAGCGCCATTCCAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCGGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGACCGCCTGAAATATCCGATGAAACGAGTCGGGGCGCGCGGTGAAGGCAAATTCGAGCGCATTA  >  NC_000913/941083‑941373

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