Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 942,080 C→T 5.7% I374I (ATC→ATT dmsA → dimethyl sulfoxide reductase, anaerobic, subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913942,0800CT 5.7% 179.6 / 2.4 71I374I (ATC→ATTdmsAdimethyl sulfoxide reductase, anaerobic, subunit A
Reads supporting (aligned to +/- strand):  ref base C (29/37);  new base T (2/2);  total (32/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

AAACCCTGCCAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCTAACGGTGAAATCGCAACCCGTGCTATCTCGATGCTGGCGATTCTGACCGGTAACGTTGGTATTAACGGAGGCAACAGCGGCGCGCGTGAAGGT  >  NC_000913/941932‑942227
                                                                                                                                                    |                                                                                                                                                   
aaaCCCTGCCAGCCAGTGCGCCGAAAAAAGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGAAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATc                                                                                                                                                     <  1:119618/149‑1 (MQ=255)
  ccccTGCCAGCCACTGCGCCGAAATAGGGCCACTATAAAGCTTATATTCTTGTGGAAGGGCCAGATGGCGTGGCTAAAAGGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCgg                                                                                                                                                   <  2:697031/148‑1 (MQ=255)
  aCCCTGCCAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCgg                                                                                                                                                   >  1:645903/1‑149 (MQ=255)
    ccTGCCAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTATGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTGCGGCAGCCAAAATCATCAAATTGGCTCGTGAAATCGgta                                                                                                                                                 >  1:903668/1‑149 (MQ=255)
             cagTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAAc                                                                                                                                        >  1:43043/1‑149 (MQ=255)
                     ctaaaaaTGGCCACTATAAAGCTTATATTCTGGGTGACGGGCCAGATGGCGGGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGttt                                                                                                                                <  2:87182/147‑1 (MQ=255)
                      gAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTa                                                                                                                               <  1:65433/149‑1 (MQ=255)
                        aaaaTGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCCCTGGTGTTCCGGCAGACAAAGTCATGAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATc                                                                                                                             >  2:1520353/1‑149 (MQ=255)
                         aaaTGGCCACTATAAAGCTTATATTCTGGGTGAAGGTACAGATGTCGTGGCTAAAACTCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCa                                                                                                                            <  1:1055104/149‑1 (MQ=255)
                         aaaTGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCa                                                                                                                            <  2:302106/149‑1 (MQ=255)
                          aaTGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAg                                                                                                                           <  1:1434193/149‑1 (MQ=255)
                           aTGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGc                                                                                                                          >  1:873946/1‑149 (MQ=255)
                            tGGCCACTAAAAAGATAATATTCTGCGTGAAGGGCCAGATGGCGTGGCTAAAACTCAGGAATGGGCCTCGCAAATCAATGGTGTTCCCGCAGACAACATCATCAAAATGTCTCGTGAAATAGGCAGTACCAAACCGGGGTTTAACAccc                                                                                                                         >  2:1225074/1‑149 (MQ=255)
                                cACTATAACGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCGCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGCCTCGTGAAATCGGTAGTACCATACCGGCGTTTATCAGCCAggg                                                                                                                     <  2:1537267/149‑1 (MQ=255)
                                  cTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGAt                                                                                                                   >  1:515949/1‑149 (MQ=255)
                                  cTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGGTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGAt                                                                                                                   <  2:1390222/149‑1 (MQ=255)
                                     tacaGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATggg                                                                                                                <  2:1409506/146‑1 (MQ=255)
                                        aGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGcc                                                                                                             <  2:262807/149‑1 (MQ=255)
                                         gCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGccc                                                                                                            <  2:636122/149‑1 (MQ=255)
                                         gCTTATATTCGGGGGGCAGGGCCAGCTGGCGTGGCTAAAACGCCGGAATGGGGCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGccc                                                                                                            <  1:562283/149‑1 (MQ=255)
                                          cTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAACACGCCGGAATGGGCCTCGCAAAACACTGGTGTTCCGGCTGACAAAATCATCAAATTGGCGCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCg                                                                                                           <  2:1285623/149‑1 (MQ=255)
                                            tataTTCTGGGTGAAGGGCCAGATGGCGGGGCTAAAAAGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAACTTGGCTCGTGAAAGCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCa                                                                                                         <  2:1478472/149‑1 (MQ=255)
                                                ttCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCCGCGt                                                                                                     >  2:389422/1‑149 (MQ=255)
                                                   tGGGTCAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGTCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCAc                                                                                                  <  1:928138/149‑1 (MQ=255)
                                                    gggTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATTGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACg                                                                                                 >  2:378659/1‑149 (MQ=255)
                                                    gggTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATTGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACg                                                                                                 >  2:384222/1‑149 (MQ=255)
                                                            ggCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCGCGCAAAGCACTGGGGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCTGCCAGGGATGGGGCCCGCAGCGTCACGCTAACGGt                                                                                         <  1:933179/149‑1 (MQ=255)
                                                               cAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTAACACCCAGGGATGGGGCCCGCAGCGTCACGCTAACGCTGaa                                                                                      >  2:844445/1‑149 (MQ=255)
                                                                  aTGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTGCGGCAGCCAAAATCATCAAATTGGCTCGTGAAATCGGTAGTAGCAAGCCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCCCGCTAACGGTGAAATc                                                                                   >  1:208929/1‑149 (MQ=255)
                                                                      cGTGGCTAAAACGCCGGAATGGTCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTCGCTCGGGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCTAACGGTGAAATCGCaa                                                                               <  1:1270331/149‑1 (MQ=255)
                                                                                 cGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCTAACGGTGAAATCGCAACCCGTGCTAtc                                                                    <  1:688232/149‑1 (MQ=255)
                                                                                  gCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCTAACGGTGAAATCGCAACCCGTGCTAtct                                                                   <  2:1426293/149‑1 (MQ=255)
                                                                                  gCCGGAATCGGCCTCGCAAATCACTGGTGTTCCGGCAGAAATAATCATCAACTTGTCTCGTGAAATCGGTAGTACCAAACGGGCGTTAATCAGCCAGGGATGGGGCCCGCAGCGTCACGCTAACGGTGAAATCGCAACCCGTGCTAtct                                                                   <  1:416991/149‑1 (MQ=255)
                                                                                   ccGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGGACCAAACCGGCGGGTATCAGCCAGGGAAGGGGCCCGCAGCGTCACGCTAACGGTGAAATCGGAACCCGTGCTAtctc                                                                  >  2:843525/1‑149 (MQ=255)
                                                                                      gAATGGGCCTCGCAAATCACTGTTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCTAACGGTGAAATCGCAACCCGTGCTATCTCGAt                                                               <  1:1032088/149‑1 (MQ=255)
                                                                                      gAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCTAACGGTGAAATCCCAACCCGTGCTATCTCGAt                                                               >  1:312234/1‑149 (MQ=255)
                                                                                       aaTGGGCCTCGCAAAGCACTGGTGTTCCGGCAGACAAATTCATCAAATGGGCTCGGGAAATCGGTAGTAGCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCTAACGGTGAAATCGCAACCCGTGCTATCTCGATg                                                              <  1:844445/149‑1 (MQ=255)
                                                                                            gCCTGGCATATCACTGGTGTTCCGGCAGACAAAATCATCAAATTCGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCTAACGGTGAAATCGCCACCCGTGCTATCTCGATTCTGGc                                                         <  2:1356473/149‑1 (MQ=255)
                                                                                              cTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTGACATACCGGCGTTTTTACGCCAGGTTTGGGTACATCAGCTTGACGCTAACTGTGAAATAGCAACCCGTGCTTTCTCGATgattgaga                                                       >  1:1014956/1‑142 (MQ=255)
                                                                                                cGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATTGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCTAACGGTGAAATCGCAACCCGTGCTATCTCGATGCTGGCGAtt                                                     <  1:378659/149‑1 (MQ=255)
                                                                                                cGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATTGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCTAACGGTGAAATCGCAACCCGTGCTATCTCGATGCTGGCGAtt                                                     <  1:384222/149‑1 (MQ=255)
                                                                                                 gCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGAGTTTATCAGCCAGGCATGGGGCCCGCAGCGTCACGCTAACGGTACTATCGCAACACGGGCTATCTCTATGCTGGCGAtta                                                    >  2:311354/1‑148 (MQ=255)
                                                                                                     aTCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCTAACGGTGAAATCGCAACCCGTGCTATCTCGATGCTGGCGATTCTGAc                                                >  1:18592/1‑149 (MQ=255)
                                                                                                      tCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCTAACGGTGAAATCGCAACCCGTGCTATCTCGATGCTGGCGATTCTGAcc                                               >  1:677790/1‑149 (MQ=255)
                                                                                                       cACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTGTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCTAACGGTGAAATCGCAACCCGTGCTATCTCGATGCTGGCGATTCTGACCg                                              <  2:1428111/149‑1 (MQ=255)
                                                                                                        aCTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCTAACGGTGAAATCGCAACCCGTGCTATCTCGATGCTGGCGATTCTGACCgg                                             <  2:472238/149‑1 (MQ=255)
                                                                                                           ggTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCTAACGGTGAAATCGCAACCCGTGCTATCTCGATGCTGGCGATTCTGACCGGTaa                                          <  1:446033/149‑1 (MQ=255)
                                                                                                           ggTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGATAACGGTGAAATCGCAACCCGTGCTATCTCGTTGCTGGCGATTCTGACCGGTaa                                          >  2:619558/1‑149 (MQ=255)
                                                                                                            gtgtTCCGGCAGACAAAATCATCAAATTGGCTCGGGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCTAACGGTGAAATCGCAACCCGTGCTATCTCGATGCTGGCGATTCTGACCGGTAAc                                         >  2:213410/1‑149 (MQ=255)
                                                                                                                 ccGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATCGGGCCAGCAGCGTAACGCTAACGGTGAAATCGCAACCCGTGATATCTCGATGCTGGCGATACTGACCGGTAACGTTgg                                    >  1:299195/1‑149 (MQ=255)
                                                                                                                  cgccAGACAAACTCATCAAATTGGCTCGTGAAATCGGTAGTAGCAAACCGGCGTTTTTCAGCCAGGGCTGGGGCCCGCAGCGTCACGATAACGGTGAAATCGCAACCCGTGCTATCTCGATGCTGGCGATTCTGACCGGTAACGTTGGt                                   <  2:1252447/146‑1 (MQ=255)
                                                                                                                  cGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCTAACGGTGAAATCGCAACCCGTGCTATCTCGATGCTGGCGATTCTGACCGGTAACGTTGGt                                   <  1:1134826/149‑1 (MQ=255)
                                                                                                                  cGGCAGACAAAATCATCAAATAGGCTCGTGAAATCGGGAGTACCAAACCGGCGTTTATCAGCCAGGGAAGGGGCCCGCAGCGTCACGCTAACTGTGAAATCGCAACCCGTGCTATCTCGATGCTGGCGATTCTGACCGGTAACGTTGGt                                   <  1:141099/149‑1 (MQ=255)
                                                                                                                   ggCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCTAACGGTGAAATCGCAACCCGTGCTATCTCGATGCTGGCGATTCTGCCCGGTAACGTTGGTa                                  >  1:672685/1‑149 (MQ=255)
                                                                                                                    gCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCTAACGGTGAAATCGCCACCCGTGCTATCTCGATGCTGGCGCTTCTGACCGGTATCGTTGGTAt                                 >  1:129107/1‑149 (MQ=255)
                                                                                                                    gCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCTAACGGTGAAATCGCAACCCGTGCTATCTCGATGCTGGCGATTCTGACCGGTAACGTTGGTAt                                 <  2:515949/149‑1 (MQ=255)
                                                                                                                    gCAGACAAAAACATCAAATTGGCGCGTGAAATCGGTAGTACCAAACCGGCGGGTATCCGCCAGGGATGGGCCCCGCCGCGTCACGCTCGGGGTGGAATCGCAACCCGTCCTAGCTCGCTGCTGGCGATTGTGCCCGGTACCGGTGGTAt                                 >  2:382261/1‑149 (MQ=255)
                                                                                                                       gACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCTAACGGTGAAATCGCAACCCGTGCTATCTCGATGCTGGCGATTCTGACCGGTAACGTTGGTATTaa                              <  2:18592/149‑1 (MQ=255)
                                                                                                                        aCAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCCCGCTAACGGTGAAATCGCAACCCGTGCTATCTCGATGCTGGCGATTCTGACCGGTAACGTTGGGATTAAc                             >  1:562631/1‑149 (MQ=255)
                                                                                                                        aCAAAATCACCAACGTGGCTCGTGAACTCGGTAGGACCAAACCGGCGTTTATCAGCCAGGGACGGGGCCCGCAGCGGCCCGCTAACCGGGAAATCGCCACCCGTGCTATCTCGATGCTGGCGATTCTGACCGGTAACGTTGGTATTAAc                             <  2:903668/149‑1 (MQ=255)
                                                                                                                             atcatcAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGAGTGGGACCGCAGCTTCAAGCTAACGGTGAAATCGCAACCCGTGCTATCTCGATGCTGGCGATTCTGACCGGTAACGTTGGTATTAACGGAgg                        <  2:873946/149‑1 (MQ=255)
                                                                                                                             atcatcAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGGTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCTAACGGTGAAATCGCAACCCGTGCTATCTCGATGCTGGCGATTCTGACCGGTAACGTTGGTATTAACGGAgg                        <  1:754294/149‑1 (MQ=255)
                                                                                                                               catcaAATTGTCTCGTGAAATCGGTAGTACCAATCCGGCGTTTATAAGACAGGGATCGGGCCCGCAGCGTCACGCTAACGGTAAAAACGCAACCCGTTATATCTCGCTGCTGGCGATTCTGACCGGTAACGTTGTTATTAACgagggac                      >  2:1021379/1‑143 (MQ=255)
                                                                                                                                atcaAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCTAACGGTGAAATCGCAACCCGTGCTATCTCGATGCTGGCGATTATGACCGGTAACGGTGTTATTAACGGAGGCaa                     >  2:908063/1‑149 (MQ=255)
                                                                                                                                atcaAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCAAGGGATGGGGCCCGCAGCGTCACGCTAACGGTGAAATCGCCACCCGTGCTATCTCGATGCTGGCGCTTCTGACCGGTAACGATGGTATTACCGGAGTCaa                     >  2:1067127/1‑149 (MQ=255)
                                                                                                                                 tcaAATTGGATCGTGAAATCGTTAGTACCCAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCAAGCTAACGGTGATATCGCAACACGTGCTATCTAGATGCTGGCGCTTCTGAACGGTAACGGTGGTATTAACGGAGGCatc                    >  1:1414463/1‑147 (MQ=255)
                                                                                                                                    aaTTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTGTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCTAACGGTGAAATCGCAACCCGTGCTATCTCGATGCTGGCGATTCTGACCGGTAACGTTGGTATTAACCGAGGCAACAgc                 >  1:91595/1‑149 (MQ=255)
                                                                                                                                     aTTGGCTCGCGAAAGCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCTAACGGTGAAATCGCAACCCGTGCTATCTCGATGCTGGCGATTCTGACCGGTAACGTTGGTATTAACGGAGGCAACAgcg                <  2:405972/149‑1 (MQ=255)
                                                                                                                                              ggAAATCGGTAGTACCAAACCGGCGTTTATCAGGCAGGGATGGGGCCCGCAGCGTCAAGCTAACGGTGAAATCGCAACCCGTGGTATCTCGATGCTGGCGATTCTGACCGGTAACGTTGGTATTAACGGAGGCAACAGCGGCGCGCGTg       <  2:33201/148‑1 (MQ=255)
                                                                                                                                               gAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCTAACGGTGAAATCGCCACCCGTGCTATCTCGATGCTGGCGATTCTGCCCGGTAACGTGGGTATTAACGGAGGCAACCGCGGCGCGCgaga      >  2:1423251/1‑149 (MQ=255)
                                                                                                                                                   tCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCTAACGGTGAAATCGCAACCCGTGCTATCTCGATGCTGGCGATTCTCACCGGTAACGTTGGTATTAACGGAGGCAACAGCGGCGCGCGTGAAGGt  >  1:1302557/1‑149 (MQ=255)
                                                                                                                                                   tCGGCAGTACCACACCGGCGTTTCTCAGCCAGGGATGGGGCCCGCATCGTCAGGCTAACGGTGAAATCGCAACCCGTGCTATCTCGATGCTGGCGATTCTGACCGGTAACGTTGGTATTAACGGAGGCAACAGCGGCGCGCGTGAAGGt  <  1:1225074/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AAACCCTGCCAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCTAACGGTGAAATCGCAACCCGTGCTATCTCGATGCTGGCGATTCTGACCGGTAACGTTGGTATTAACGGAGGCAACAGCGGCGCGCGTGAAGGT  >  NC_000913/941932‑942227

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