Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 943,269 G→A 5.1% G771S (GGT→AGT)  dmsA → dimethyl sulfoxide reductase, anaerobic, subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913943,2690GA 5.1% 252.6 / 2.2 78G771S (GGT→AGT) dmsAdimethyl sulfoxide reductase, anaerobic, subunit A
Reads supporting (aligned to +/- strand):  ref base G (38/35);  new base A (2/2);  total (41/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

GGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGTACTGACCACTCAACGTCCGTCTCCTCTCGCTAAGGGGAATCCGTCACATACAAACCTTGTTCAGGTTGAAAAGGTGTAAGGAGTAA  >  NC_000913/943124‑943410
                                                                                                                                                 |                                                                                                                                             
ggAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGt                                                                                                                                             >  1:339911/1‑148 (MQ=255)
  aaaTGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTg                                                                                                                                            <  5:187710/147‑1 (MQ=255)
  aaaTGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTg                                                                                                                                            >  6:187710/1‑147 (MQ=255)
  aaaTGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGc                                                                                                                                           >  5:29088/1‑148 (MQ=255)
       tgGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGt                                                                                                                                      <  2:339911/148‑1 (MQ=255)
             aaCCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGAcc                                                                                                                                <  1:85274/148‑1 (MQ=255)
                ccGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGTCTGGTATGACCCg                                                                                                                              >  2:319705/1‑147 (MQ=255)
                ccGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCgg                                                                                                                             <  6:139780/148‑1 (MQ=255)
                ccGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCg                                                                                                                              >  5:365668/1‑147 (MQ=255)
                ccGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTAATAGTGAGGCGAAAGTGACGCCACGAATGATGCCTGGTGTGGTCGCACTGGGTGAATGTGCCTGGTATGACCCg                                                                                                                              >  5:372001/1‑147 (MQ=255)
                             aaaaaCGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTa                                                                                                                                     <  2:105165/127‑1 (MQ=255)
                             aaaaaCGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTa                                                                                                                                     >  1:105165/1‑127 (MQ=255)
                              aaaaCGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTc                                                                                                               >  3:331962/1‑148 (MQ=255)
                                aaCGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGa                                                                                                             <  4:331962/148‑1 (MQ=255)
                                   gcgGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGATGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACGGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATaa                                                                                                          <  2:147142/148‑1 (MQ=255)
                                          ccACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCaaa                                                                                                                      <  6:272930/129‑1 (MQ=255)
                                          ccACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCaaa                                                                                                                      >  5:272930/1‑129 (MQ=255)
                                          ccACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGggtggt                                                                                                   >  7:455432/1‑148 (MQ=255)
                                          ccACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGggtggt                                                                                                   >  1:299520/1‑148 (MQ=255)
                                          ccACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGggtggt                                                                                                   >  1:104924/1‑148 (MQ=255)
                                           cacaACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGggtggt                                                                                                   <  8:455432/147‑1 (MQ=255)
                                            acaACGGCGATAAAGTCAGGATCTTTAACGATCGTTGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGGGGTCGCACTGGGTGAAGGGGCCTGGTATGACCCGGAGGCaa                                                                                                                       <  5:106477/126‑1 (MQ=255)
                                            acaACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTg                                                                                                 >  8:292770/1‑148 (MQ=255)
                                             caACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGTTGCCTGGTATGACCCGGATGCaaa                                                                                                                      >  6:106477/1‑126 (MQ=255)
                                             caACGGAGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTg                                                                                                 >  2:244160/1‑147 (MQ=255)
                                               aCGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCAAGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTAt                                                                                              <  8:102093/148‑1 (MQ=255)
                                                  gCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTaa                                                                                           >  2:218330/1‑148 (MQ=255)
                                                     aTAAAGTCAGGATGTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGggtggt                                                                                                   <  7:345718/137‑1 (MQ=255)
                                                     aTAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGggtggt                                                                                                   >  8:345718/1‑137 (MQ=255)
                                                      tAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGt                                                                                        <  2:279816/147‑1 (MQ=255)
                                                      tAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGt                                                                                        >  1:279816/1‑147 (MQ=255)
                                                        aaGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGTAc                                                                                      <  7:221462/147‑1 (MQ=255)
                                                          gTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGTACTGa                                                                                   <  2:200626/148‑1 (MQ=255)
                                                              ggATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAATGGTGGTTGTATTAACGTACTGACAAc                                                                               >  8:94284/1‑148 (MQ=255)
                                                                aTCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGTACTGACCACt                                                                              >  3:149821/1‑147 (MQ=255)
                                                                     tAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGTGCTGACCACTc                                                                             <  1:244160/143‑1 (MQ=255)
                                                                          aTCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGTACTGACCACTCAAcgtccgtc                                                                   >  8:47052/1‑148 (MQ=255)
                                                                           tCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGTACTGACCACTCAACGTCCGtct                                                                  >  1:54462/1‑148 (MQ=255)
                                                                                     gTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGTACTGACCACTCAACGTCCGTCTCCTCTCGCTa                                                        <  2:104924/148‑1 (MQ=255)
                                                                                      ttCATATTGAGGCGACACTGCCCCCACGTATGATGCCGGGGTGGGGCGCACTGGGTGAAGGTTCCTGGTATGACCCGGATGCAAAACGTGTCGATAACGGTGGCTGTATTAACGTACTGACCACTCAACGTCCGTCTCCTCTCGCTaa                                                       <  6:29088/148‑1 (MQ=255)
                                                                                      ttCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGTACTGACCACTCAACGTCCGTCTCCTCTCGCt                                                         >  8:302539/1‑146 (MQ=255)
                                                                                      ttCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGTACTGACCACTCAACGTCCGTCTCCTCTCGCTaa                                                       <  6:365668/148‑1 (MQ=255)
                                                                                      ttCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTGGTATTAACGTACTGACCACTCAACGTCCGTCTCCTCTCGCTaa                                                       <  6:372001/148‑1 (MQ=255)
                                                                                       tCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGt                                                                                        <  6:227553/114‑1 (MQ=255)
                                                                                       tCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGt                                                                                        >  5:227553/1‑114 (MQ=255)
                                                                                        cATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGTACTGACCACTCAACGTCCGTCTCCTCTCGCTAAgg                                                     >  3:166977/1‑148 (MQ=255)
                                                                                        cATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAAGTGCCTGGTATGACCCGGATGCAAAACgtgt                                                                                                                >  7:348558/1‑89 (MQ=255)
                                                                                        cATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAAGTGCCTGGTATGACCCGGATGCAAAACgtgt                                                                                                                >  7:342373/1‑89 (MQ=255)
                                                                                        cATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAAGTGCCTGGTATGACCCGGATGCAAAACgtgt                                                                                                                <  8:342373/89‑1 (MQ=255)
                                                                                        cATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAAGTGCCTGGTATGACCCGGATGCAAAACgtgt                                                                                                                <  8:348558/89‑1 (MQ=255)
                                                                                         atatTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGTCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGTACTGACCACTCAACGTCCGTCTCCTCTCGCTAAgg                                                     <  1:319705/147‑1 (MQ=255)
                                                                                             tGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGTACTGACCACTCAACGTCCGTCTCCTCTCGCTAAGGGGAAt                                                >  4:379609/1‑148 (MQ=255)
                                                                                             tGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGTACTGACCACTCAACGTCCGTCTCCTCTCGCTAAGGGGAAt                                                >  6:375806/1‑148 (MQ=255)
                                                                                               aGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGTACTGACCACTCAACGTCCGTCTCCTCTCGCTAAGGGGAATc                                               <  7:292770/147‑1 (MQ=255)
                                                                                                    aaaGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGTACTGACCACTCAACGTCCGTCTCCTCTCGCTAAGGGGAATCCGTca                                          >  7:97539/1‑147 (MQ=255)
                                                                                                     aaGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGTACTGACCACTCAACGTCCGTCTCCTCTCGCTAAGGGGAATCCGTcaca                                        >  4:60542/1‑148 (MQ=255)
                                                                                                             ccACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGTACTGACCACTCAACGTCCGTCTCCTCTCGCTAAGGGGAATCCGTCACATACAAAcc                                >  8:14443/1‑148 (MQ=255)
                                                                                                              cACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGTACTGACCACTCAACGTCCGTCTCCTCTCGCTAAGGGGAATCCGTCACATACAAACCt                               <  2:299520/148‑1 (MQ=255)
                                                                                                               aCGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGTACTGACCACTCAACGTCCGTCTCCTCTCGCTAAGGGGAATCCGTCACATACAAACCtt                              >  7:394482/1‑148 (MQ=255)
                                                                                                                 gAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGTACTGACCACTCAACGTCCGTCTCCTCTCGCTAAGGGGAATCCGTCACATACAAACCTTGt                            <  5:416576/148‑1 (MQ=255)
                                                                                                                   atgatgCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGTACTGACCACTCAACGTCCGTCTCCTCTCGCTAAGGGGAATCCGTCACATACAAACCTTGTTc                          <  4:149821/148‑1 (MQ=255)
                                                                                                                        gCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGTACTGACCACTCAACGTCCGTCTCCTCTCGCTAAGGGGAATCCGTCACATACAAACCTTGTTCAGGtt                     >  3:8435/1‑148 (MQ=255)
                                                                                                                        gCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGTACTGACCACTCAACGTCCGTCTCCTCTCGCTAAGGGGAATCCGTCACATACAAACCTTGTTCAGGtt                     >  3:2618/1‑148 (MQ=255)
                                                                                                                            ggTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGTACTGACCACTCAACGTCCGTCTCCTCTCGCTAAGGGGAATCCGTCACATACAAACCTTGTTCAGGTTGaa                  >  8:336205/1‑147 (MQ=255)
                                                                                                                                 ggTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGTACTGACCACTCAACGTCCGTCTCCTCTCGCTAAGGGGAATCCGTCACATACAAACCTTGTTCAGGTTGAAAAGgtg            >  3:231978/1‑148 (MQ=255)
                                                                                                                                  gTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGTACTGACCACTCAAcgtcc                                                                      >  3:371602/1‑89 (MQ=255)
                                                                                                                                  gTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGTACTGACCACTCAAcgtcc                                                                      <  4:371602/89‑1 (MQ=255)
                                                                                                                                   tCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGTACTGACCACTCAACGTCCGTCTCCTCTCGCTAAGGGGAATCCGTCACATACAAACCTTGTTCAGGTTGAAAAGgtgt           >  8:237255/1‑147 (MQ=255)
                                                                                                                                   tCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGTACTGACCACTCAACGTCCGTCTCCTCTCGCTAAGGGGAATCCGTCACATACAAACCTTGTTCAGGTTGAAAAGGTGTa          <  2:54462/148‑1 (MQ=255)
                                                                                                                                   tCGCACTGGGTGAAGGTGCCTGGTATGACCCAGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGTACTGACCACTCAACGTCCGTCTCCTCTCGCTAAGGGGAATCCGTCACATACAAACCTTGTTCAGGTTGaa                  <  4:145761/140‑1 (MQ=255)
                                                                                                                                   tCGCACTGGGTGAAGGTGCCTGGTATGACCCAGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGTACTGACCACTCAACGTCCGTCTCCTCTCGCTAAGGGGAATCCGTCACATACAAACCTTGTTCAGGTTGaa                  >  3:145761/1‑140 (MQ=255)
                                                                                                                                    cGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGTACTGACCACTCAACGTCCGTCTCCTCTCGCTAAGGGGAATCCGTCACATACAAACCTTGTTCAGGTTGAAAAGGTGTaa         <  5:371712/148‑1 (MQ=255)
                                                                                                                                     gCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGTACTGACCACTCAACGTCCGTCTCCTCTCGCTAAGGGGAATCCGTCACATACAAACCTTGTTCAGGTTGAAAAGGTGTaa         <  7:336205/147‑1 (MQ=255)
                                                                                                                                     gCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGTACTGACCACTCAACGTCCGTCTCCTCTCGCTAAGGGGAATCCGTCACATACAAACCTTGTTCAGGTTGAAAAGGTGTaa         <  7:14443/147‑1 (MQ=255)
                                                                                                                                      cACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGTACTGACCACTCAACGTCCGTCTCCTCTCGCTAAGGGGAATCCGTCACATACAAACCTTGTTCAGGTTGAAAAGGTGTAAgg       <  7:237255/148‑1 (MQ=255)
                                                                                                                                           ggTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGTACTGACCACTCAACGTCCGTCTCCTCTCGCTAAGGGGAATCCGTCACATACAAACCTTGTTCAGGTTGAAAAGGTGTAAGGAGTaa  <  4:166977/148‑1 (MQ=255)
                                                                                                                                             tGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGTACTGACCACTCAACGTCCGTCTCCTCTCGCTAAGGGGAATCCGTca                                          >  6:291647/1‑106 (MQ=255)
                                                                                                                                             tGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGTACTGACCACTCAACGTCCGTCTCCTCTCGCTAAGGGGAATCCGTca                                          <  5:291647/106‑1 (MQ=255)
                                                                                                                                                 ggTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGTACTGACCACTCAACGTCCGTCTCCTCTCGCTAAGGGGAATCCGTCACATACAAACCTTGTTCAgg                       >  1:363246/1‑121 (MQ=255)
                                                                                                                                                 ggTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGTACTGACCACTCAACGTCCGTCTCCTCTCGCTAAGGGGAATCCGTCACATACAAACCTTGTTCAgg                       <  2:363246/121‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
GGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCCTGGTATGACCCGGATGCAAAACGTGTCGATAAGGGTGGTTGTATTAACGTACTGACCACTCAACGTCCGTCTCCTCTCGCTAAGGGGAATCCGTCACATACAAACCTTGTTCAGGTTGAAAAGGTGTAAGGAGTAA  >  NC_000913/943124‑943410

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