Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 942,026 C→G 5.5% A356A (GCC→GCG 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,0260CG 5.5% 231.9 / 2.6 68A356A (GCC→GCGdmsAdimethyl sulfoxide reductase, anaerobic, subunit A
Reads supporting (aligned to +/- strand):  ref base C (30/34);  new base G (2/2);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCTAACGGTGAAATCGCAACCCGTGCTATCTCG  >  NC_000913/941880‑942164
                                                                                                                                                  |                                                                                                                                          
aaCCTTGTGGATCATGAATTCCAAGATAAATATTGCGTTCGCAACGTTGACAAAACCCTGCCATCCATTGCGCCGAAAAATGGCCACTATAAATCTTATAATCTGGGTTAAGGGCCAGATTGCGTGGCTAAAACGCCGGAATTGGCCt                                                                                                                                           <  4:358752/148‑1 (MQ=255)
   cTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCg                                                                                                                                         >  3:329963/1‑147 (MQ=255)
   cTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCg                                                                                                                                         >  7:48057/1‑147 (MQ=255)
       tggATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGATAAAACCCTGCCAGCCAGTGAGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGACGGAATGGGCCTCGCAAAt                                                                                                                                    >  7:52263/1‑148 (MQ=255)
          aTCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCAc                                                                                                                                 <  4:247131/148‑1 (MQ=255)
             aGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTgg                                                                                                                              <  7:43934/148‑1 (MQ=255)
                  ttCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTc                                                                                                                         >  5:199128/1‑148 (MQ=255)
                   tCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTcc                                                                                                                        >  5:424651/1‑148 (MQ=255)
                   tCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTcc                                                                                                                        >  4:354851/1‑148 (MQ=255)
                    ccTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCGTCGCAAATCACTGGTGTTCCg                                                                                                                       >  6:338733/1‑148 (MQ=255)
                        gATAAATATTGCGTTGGCTACTATGAGAAAACCCTGCCAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAg                                                                                                                   <  7:326740/148‑1 (MQ=255)
                              tatTGCGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACaaaa                                                                                                             >  6:370438/1‑148 (MQ=255)
                                   cGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAatcat                                                                                                         <  1:18359/147‑1 (MQ=255)
                                       ggCTACGATGAGAAAACCCTGCCAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCGTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCaaa                                                                                                     <  5:338733/147‑1 (MQ=255)
                                              aTGAGAAAACCCTGCCAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCg                                                                                             >  3:260449/1‑148 (MQ=255)
                                              aTGAGAAAACCCTGCCAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAAGGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCg                                                                                             >  6:40907/1‑148 (MQ=255)
                                                 agaAAACCCTGCCAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTg                                                                                           >  5:158816/1‑147 (MQ=255)
                                                   aaaaCCCTGCCAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTTaaa                                                                                        <  2:149189/148‑1 (MQ=255)
                                                   aaaaCCCTGCCAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGaaa                                                                                        >  7:145488/1‑148 (MQ=255)
                                                       cacTGCCAGCCAGTGCGCCGAAAAATGGCAACTATAAAGCTTATATTCTGTGTGAAGGGCCAGATGGCTTGGCTAAAACGCAGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATAAAATTGGCTCGTGAAATCgg                                                                                    <  8:52263/146‑1 (MQ=255)
                                                        ccTGCCAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGgt                                                                                   <  6:325090/148‑1 (MQ=255)
                                                        ccTGCCAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGgt                                                                                   <  6:199128/148‑1 (MQ=255)
                                                        ccTGCCAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGgt                                                                                   <  6:386614/148‑1 (MQ=255)
                                                        ccTGCCAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGgt                                                                                   <  6:424651/148‑1 (MQ=255)
                                                           gccagccaGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGgtag                                                                                 >  4:112257/1‑147 (MQ=255)
                                                           gccagccaGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGgtag                                                                                 >  7:150898/1‑147 (MQ=255)
                                                            ccagccagTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGgtagta                                                                               <  7:170883/148‑1 (MQ=255)
                                                             cagccagTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTAc                                                                              >  3:143806/1‑148 (MQ=255)
                                                                     gcgcCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTAGTGAAATCGGTAGTACCAAACCgg                                                                      <  4:329963/148‑1 (MQ=255)
                                                                      cgcCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCgg                                                                      <  8:48057/147‑1 (MQ=255)
                                                                        ccGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGt                                                                   <  4:260449/148‑1 (MQ=255)
                                                                         cGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGtt                                                                  <  3:112257/148‑1 (MQ=255)
                                                                          gAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGttt                                                                 <  1:4982/148‑1 (MQ=255)
                                                                          gAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGtt                                                                  >  5:210875/1‑147 (MQ=255)
                                                                                 ggCCACTATAAAGCTTATATTCTGGGTGAAGGCCCAGATGGCGTGGCTAAAAAGCCGGAATGGGCCTCGCAAATAACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGc                                                           >  5:331516/1‑147 (MQ=255)
                                                                                  gCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGcc                                                          <  3:354851/147‑1 (MQ=255)
                                                                                    cACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAgg                                                       <  8:145488/148‑1 (MQ=255)
                                                                                        ataAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATg                                                   >  2:170679/1‑148 (MQ=255)
                                                                                         taAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATg                                                   <  5:370438/147‑1 (MQ=255)
                                                                                           aaGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATgggg                                                >  8:271606/1‑148 (MQ=255)
                                                                                             gCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGaaa                                                                                        <  3:106979/106‑1 (MQ=255)
                                                                                             gCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGaaa                                                                                        >  4:106979/1‑106 (MQ=255)
                                                                                              cTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGCCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGccc                                             >  2:296585/1‑148 (MQ=255)
                                                                                               ttATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCg                                            >  2:259731/1‑148 (MQ=255)
                                                                                                         ggTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTc                                     <  8:150898/145‑1 (MQ=255)
                                                                                                         ggTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCAc                                   <  1:259731/147‑1 (MQ=255)
                                                                                                          gTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGc                                 >  1:132704/1‑148 (MQ=255)
                                                                                                                  ccAGATGGCGTTGCTAATAAGCCGGATTGGGCCTCGCTAATCTCTTGGTTTCGGGCAGACAAAATAATCAAATTGTGTCGTGTAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTc                                     <  5:153891/136‑1 (MQ=255)
                                                                                                                  ccAGATGGCGTGGCTAAAACGCCGGAAGGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAACTTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCTAACGGTg                         <  2:421830/148‑1 (MQ=255)
                                                                                                                      aTGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCgg                                                                                    <  3:133641/85‑1 (MQ=255)
                                                                                                                      aTGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCgg                                                                                    >  4:133641/1‑85 (MQ=255)
                                                                                                                      aTGGCGTGGCTAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGgt                                                                                   <  1:257917/85‑1 (MQ=255)
                                                                                                                      aTGGCGTGGCTAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGgt                                                                                   >  2:257917/1‑85 (MQ=255)
                                                                                                                       tGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAAc                                                                         <  5:196239/95‑1 (MQ=255)
                                                                                                                       tGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCTAACGGTGAAATc                    <  1:170679/148‑1 (MQ=255)
                                                                                                                       tGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATAAAATTGGCTCGTGAAATCGGTAGTACCAAAc                                                                         >  6:196239/1‑95 (MQ=255)
                                                                                                                         gCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGTTGTTCCGTCAGACAAAATCATCAAATTCGCTCGTGAAAGCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCTAACGGTGAAATCg                   <  6:331516/147‑1 (MQ=255)
                                                                                                                             ggCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCTAACGGTGAAATCGCAAcc              <  6:158816/148‑1 (MQ=255)
                                                                                                                               cTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCTAACGGTGAAATCGCAACCCg            >  4:138526/1‑148 (MQ=255)
                                                                                                                                tAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAg                                                            >  5:270747/1‑99 (MQ=255)
                                                                                                                                tAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAg                                                            <  6:270747/99‑1 (MQ=255)
                                                                                                                                 aaaaCGCCGGAATGGGCGTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCTAACGGTGAAATCGCaa                <  5:385652/142‑1 (MQ=255)
                                                                                                                                 aaaaCGCCGGAATGGGCGTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCTAACGGTGAAATCGCaa                >  6:385652/1‑142 (MQ=255)
                                                                                                                                 aaaaCGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCTAACGGTGAAATCGCAACCCGTg          >  2:257977/1‑148 (MQ=255)
                                                                                                                                     cGCCGGAATGGGCCTCGCAAATCACTGGTGGTCCGGCAGACAAAATCATCAAATTGGCTCGTGCAATCGGTAGTACCAACCCGCCGTTGATCAGCCAGGGAGGGGGCCCGCCGCGTCACGCTGACGGTGCAAGCGCACCCCGGGCtgc      >  8:173723/1‑146 (MQ=255)
                                                                                                                                        cGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATg                                                   >  4:102907/1‑100 (MQ=255)
                                                                                                                                        cGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATg                                                   <  3:102907/100‑1 (MQ=255)
                                                                                                                                         ggAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCTAACGGTGAAATCGCAACCCGTGCTATCTCg  <  6:210875/148‑1 (MQ=255)
                                                                                                                                                  cTCCCAAATCACTGGTTTTCCGGCAGACAAAATCATCAAATTGGCGCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCTAACGGt                          <  8:298341/115‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                          
AACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCTAACGGTGAAATCGCAACCCGTGCTATCTCG  >  NC_000913/941880‑942164

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