Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 216,145 (C)6→7 intergenic (+166/+34) cutF → / ← yaeF copper homeostasis protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913216,1391.C97.2% 657.7 / 10.5 212intergenic (+160/+40)cutF/yaeFcopper homeostasis protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (3/3);  new base C (114/92);  total (117/95)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

Only 100 of 238 total aligned reads displayed.

CGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGG  >  NC_000913/216047‑216226
                                                                                              |                                                                                        
     tttGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAccccc                                                                                    <  2:269003/93‑1 (MQ=255)
       tGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAccccccc                                                                                  >  2:1361161/1‑92 (MQ=255)
       tGTAGGCCGAATAAGGCCTTACCGCCGCATCCGGCAATGGTGCTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAccccccc                                                                                  <  2:3801770/93‑2 (MQ=18)
          aGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAccccccctca                                                                               <  1:1222102/93‑5 (MQ=255)
          aGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAccccccctca                                                                               <  2:6000270/93‑5 (MQ=255)
             ccGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATAc                                                                             >  1:5086193/1‑93 (MQ=255)
             ccGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATAc                                                                             >  1:6296433/1‑93 (MQ=255)
             ccGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATAc                                                                             <  2:4424075/93‑1 (MQ=255)
              cGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATAcc                                                                            >  2:5355549/1‑93 (MQ=255)
               ggATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCt                                                                           >  2:6065976/1‑93 (MQ=255)
               ggATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCt                                                                           >  1:1547260/1‑93 (MQ=255)
                gATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCtt                                                                          <  2:2300351/93‑1 (MQ=255)
                gATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCtt                                                                          >  2:3082884/1‑93 (MQ=255)
                gATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCtt                                                                          >  1:2354872/1‑93 (MQ=255)
                 aTAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGTGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTa                                                                         >  1:496494/1‑93 (MQ=255)
                 aTAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTa                                                                         <  2:1439135/93‑1 (MQ=255)
                  tAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTaa                                                                        >  2:5606782/1‑93 (MQ=255)
                  tAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTaa                                                                        <  2:5929805/93‑1 (MQ=255)
                   aaGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAg                                                                       >  2:3640025/1‑93 (MQ=255)
                   aaGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAg                                                                       <  2:4276129/93‑1 (MQ=255)
                   aaGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAg                                                                       >  1:712186/1‑93 (MQ=255)
                   aaGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAg                                                                       <  1:539884/93‑1 (MQ=255)
                   aaGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAg                                                                       <  2:80303/93‑1 (MQ=255)
                   aaGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAg                                                                       <  1:2696581/93‑1 (MQ=255)
                   aaGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAg                                                                       >  1:2802791/1‑93 (MQ=255)
                   aaGGCGTATACGCCTCATCCGGCAATGGTGCTCATCGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAg                                                                       <  1:1820564/93‑1 (MQ=21)
                     ggCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGc                                                                     <  1:601077/93‑1 (MQ=255)
                       cGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCgg                                                                   >  1:5624125/1‑93 (MQ=255)
                       cGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCgg                                                                   <  2:3416030/93‑1 (MQ=255)
                         tttACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGAt                                                                 >  2:2012774/1‑93 (MQ=255)
                         tttACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGAt                                                                 >  1:2761968/1‑93 (MQ=255)
                          ttACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATa                                                                >  1:1330642/1‑93 (MQ=255)
                           tACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATaa                                                               <  1:3253090/93‑1 (MQ=255)
                             cgccccATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAgg                                                             <  1:192080/93‑1 (MQ=255)
                               ccgcATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCa                                                           >  2:2974905/1‑93 (MQ=255)
                               ccgcATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCa                                                           >  1:4814489/1‑93 (MQ=255)
                               ccgcATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCATAAGGCGGATAACGCa                                                           >  2:6211235/1‑93 (MQ=25)
                                cgcATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCAt                                                          <  1:2634540/93‑1 (MQ=255)
                                cgcATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCAt                                                          <  2:1050518/93‑1 (MQ=255)
                                 gcATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCAtt                                                         >  2:2876698/1‑93 (MQ=255)
                                 gcATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCAtt                                                         >  2:3497324/1‑93 (MQ=255)
                                  cATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCAttt                                                        <  2:1654366/93‑1 (MQ=255)
                                   aTCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTa                                                       >  1:385458/1‑93 (MQ=255)
                                    tCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTAc                                                      >  2:2491374/1‑93 (MQ=255)
                                      cGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACtt                                                    <  2:5316336/93‑1 (MQ=255)
                                      cGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACtt                                                    >  2:1679337/1‑93 (MQ=255)
                                      cGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACtt                                                    <  1:1411397/93‑1 (MQ=255)
                                       ggCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACttt                                                   >  1:5765548/1‑93 (MQ=255)
                                       ggCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACttt                                                   <  2:687296/93‑1 (MQ=255)
                                        gCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATTATGCCTACAACCCCCCCTCATACCGTAAGGCTGATAAAGCATTTACTTTa                                                  >  2:4614544/1‑93 (MQ=25)
                                        gCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTa                                                  <  1:4518655/93‑1 (MQ=255)
                                        gCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTa                                                  <  2:3337821/93‑1 (MQ=255)
                                         cAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTAt                                                 >  2:3897223/1‑93 (MQ=255)
                                          aaTGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATc                                                >  1:1035960/1‑93 (MQ=255)
                                          aaTGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATc                                                >  2:4220612/1‑93 (MQ=255)
                                           aTGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCa                                               >  1:353202/1‑93 (MQ=255)
                                           aTGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCa                                               <  2:541659/93‑1 (MQ=255)
                                             ggTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCAcc                                             <  2:584095/93‑1 (MQ=255)
                                              gTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCg                                            <  1:5549249/93‑1 (MQ=255)
                                                 cTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGt                                         <  2:5008782/93‑1 (MQ=255)
                                                  tCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTc                                        >  1:5256833/1‑93 (MQ=255)
                                                   cAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCa                                       >  1:3946567/1‑93 (MQ=255)
                                                   cAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCa                                       >  1:2246852/1‑93 (MQ=255)
                                                    aaCGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCaa                                      >  1:4377901/1‑93 (MQ=255)
                                                    aaCGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCaa                                      >  2:5531483/1‑93 (MQ=255)
                                                     aCGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCaaa                                     >  1:572675/1‑93 (MQ=255)
                                                       gCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAAcc                                   <  2:2542493/93‑1 (MQ=255)
                                                         cTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACcaa                                 >  2:4765478/1‑93 (MQ=255)
                                                          tGATGCGACGCTGGCGCGTCTTAGCATGCCTACAACCCCCCCTCATGCCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGGCAAACcaac                                >  2:5073387/1‑93 (MQ=37)
                                                           gATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACcaaca                               <  1:2536842/93‑1 (MQ=255)
                                                           gATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACcaaca                               >  1:2026982/1‑93 (MQ=255)
                                                           gATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACcaaca                               >  1:4657611/1‑93 (MQ=255)
                                                            aTGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACcaacaa                              >  1:900125/1‑93 (MQ=255)
                                                            aTGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCAGGTCAAACcaacaa                              <  2:1223530/93‑1 (MQ=255)
                                                             tGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACaaa                             >  1:4308542/1‑93 (MQ=255)
                                                               cGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGc                           >  2:3141559/1‑93 (MQ=255)
                                                                  cGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGAc                        >  1:3012015/1‑93 (MQ=255)
                                                                  cGCTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGAc                        <  2:5714834/93‑1 (MQ=255)
                                                                    cTGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGAccc                      >  1:3597003/1‑93 (MQ=255)
                                                                     tGGCGCGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCg                     >  1:4061529/1‑93 (MQ=255)
                                                                       gcgcGTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGcc                   <  1:4219328/93‑1 (MQ=255)
                                                                          cgTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTg                <  1:502765/93‑1 (MQ=255)
                                                                           gTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGa               >  2:1875244/1‑93 (MQ=255)
                                                                           gTCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGa               >  1:4173199/1‑93 (MQ=255)
                                                                            tCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGat              >  1:107869/1‑93 (MQ=255)
                                                                            tCTTATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGat              >  1:5957114/1‑93 (MQ=255)
                                                                              ttATCATGCCTACAACCCTCCCTCATACCTAAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGatat            >  2:1729469/1‑93 (MQ=255)
                                                                               tATCATGCCTACACCCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGAtata           <  2:2743826/93‑1 (MQ=255)
                                                                               tATCATGCCTACACCCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGAtata           <  1:6833/93‑1 (MQ=255)
                                                                               tATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGAtata           >  1:3410498/1‑93 (MQ=255)
                                                                               tATCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGAtata           <  1:1721673/93‑1 (MQ=255)
                                                                                aTCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATaa          >  2:4757280/1‑93 (MQ=255)
                                                                                 tCATGCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATaaa         >  1:4698142/1‑93 (MQ=255)
                                                                                     gCCTACACCCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATccc     <  2:3547484/93‑1 (MQ=255)
                                                                                     gCCTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATccc     <  1:2138448/93‑1 (MQ=255)
                                                                                      ccTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCa    <  1:1434073/93‑1 (MQ=255)
                                                                                      ccTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCa    >  1:5796135/1‑93 (MQ=255)
                                                                                      ccTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCa    <  1:2857368/93‑1 (MQ=255)
                                                                                      ccTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCa    <  1:311824/93‑1 (MQ=255)
                                                                                      ccTACAACCCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCa    <  1:1304562/93‑1 (MQ=255)
                                                                                              |                                                                                        
CGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGG  >  NC_000913/216047‑216226

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