Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 675,063 A→G A236A (GCT→GCC ybeQ ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913675,0630AG99.7% 965.8 / ‑5.0 303A236A (GCT→GCCybeQhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (169/132);  total (170/133)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01

Only 100 of 318 total aligned reads displayed.

CCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGCAATTTGAAACTG  >  NC_000913/674973‑675152
                                                                                          |                                                                                         
ccATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGc                                                                                         <  1:2045061/93‑1 (MQ=255)
   ttCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGggcggc                                                                                      >  1:3592464/1‑93 (MQ=255)
   ttCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGggcggc                                                                                      <  2:8687064/93‑1 (MQ=255)
    tCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCa                                                                                     >  1:951167/1‑93 (MQ=255)
     cAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAg                                                                                    >  2:5065805/1‑93 (MQ=255)
      aaaGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGc                                                                                   >  1:296446/1‑93 (MQ=255)
       aaGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGcc                                                                                  <  2:4808088/93‑1 (MQ=255)
           cAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCa                                                                              <  2:6210729/93‑1 (MQ=255)
           cAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCa                                                                              <  2:503278/93‑1 (MQ=255)
            aGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAg                                                                             >  1:1722120/1‑93 (MQ=255)
            aGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAg                                                                             <  1:6028876/93‑1 (MQ=255)
             gCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGa                                                                            >  2:9934570/1‑93 (MQ=255)
             gCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGa                                                                            <  1:534721/93‑1 (MQ=255)
             gCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGa                                                                            <  2:5792284/93‑1 (MQ=255)
             gCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGa                                                                            >  1:464451/1‑93 (MQ=255)
             gCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGa                                                                            <  1:381577/93‑1 (MQ=255)
             gCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATAAATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGa                                                                            >  2:4604960/1‑93 (MQ=255)
              ccTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGAt                                                                           <  2:3718485/93‑1 (MQ=255)
                tGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATAc                                                                         >  2:8810644/1‑93 (MQ=255)
                tGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATAc                                                                         >  1:3997802/1‑93 (MQ=255)
                tGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATAc                                                                         >  2:9142942/1‑93 (MQ=255)
                tGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATAc                                                                         >  1:5322293/1‑93 (MQ=255)
                 gATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATAcc                                                                        <  1:4340218/93‑1 (MQ=255)
                   tAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAg                                                                      >  1:6952970/1‑93 (MQ=255)
                   tAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAg                                                                      >  2:2133323/1‑93 (MQ=255)
                    aaTCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAga                                                                     >  1:1004700/1‑92 (MQ=255)
                    aaTCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGt                                                                     >  2:5387524/1‑93 (MQ=255)
                     aTCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTa                                                                    >  1:9182850/1‑93 (MQ=255)
                      tCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCTCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTAc                                                                   >  2:181806/1‑93 (MQ=255)
                      tCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCGGCCTTCAGATACCAGTAc                                                                   >  1:4238611/1‑93 (MQ=255)
                       cTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACa                                                                  >  2:1982153/1‑93 (MQ=255)
                        tttCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACAt                                                                 >  1:4785381/1‑93 (MQ=255)
                        tttCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACAt                                                                 <  2:1113991/93‑1 (MQ=255)
                        tttCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACAt                                                                 <  1:3418936/93‑1 (MQ=255)
                         ttCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGAGGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATc                                                                >  1:4448816/1‑93 (MQ=255)
                         ttCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATc                                                                <  1:4893105/93‑1 (MQ=255)
                           ctcAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCgc                                                              <  2:7142040/93‑1 (MQ=255)
                            tcAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGTCTTCAGATACCAGTACATCgct                                                             >  2:2369146/1‑93 (MQ=255)
                            tcAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCgct                                                             >  2:3053948/1‑93 (MQ=255)
                            tcAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCCTCAGATACCAGTACATCgct                                                             >  2:3481859/1‑93 (MQ=255)
                              aaCGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCgcttg                                                           >  1:9935677/1‑93 (MQ=255)
                               aCGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCgcttgc                                                          >  2:9596141/1‑93 (MQ=255)
                                  ccTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTta                                                       <  1:9280673/93‑1 (MQ=255)
                                    ttGTCCGTGTTTATACATATAACCGATGTTGACGTACGCGCCGACGCTAACCTGGGCGGGAGCCTTCCGATACCAGTACATCGGTTGCTtatc                                                     >  2:7949919/1‑92 (MQ=255)
                                    ttGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTtata                                                     >  2:8158602/1‑93 (MQ=255)
                                    ttGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTtata                                                     >  2:4482877/1‑93 (MQ=255)
                                    ttGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTtata                                                     >  2:5567187/1‑93 (MQ=255)
                                      gTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGt                                                   <  2:329279/93‑1 (MQ=255)
                                       tCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTc                                                  >  1:7513798/1‑93 (MQ=255)
                                         cGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCt                                                <  1:3012870/93‑1 (MQ=255)
                                         cGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCt                                                <  1:9348636/93‑1 (MQ=255)
                                           tgtTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTgg                                              <  2:7489403/93‑1 (MQ=255)
                                             tttATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTc                                            >  1:1298824/1‑93 (MQ=255)
                                              ttATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTcc                                           >  2:7266235/1‑93 (MQ=255)
                                              ttATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTcc                                           <  1:2339111/93‑1 (MQ=255)
                                                ataCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCAc                                         >  2:1029811/1‑93 (MQ=255)
                                                ataCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCAc                                         >  1:6524371/1‑93 (MQ=255)
                                                ataCATATAACAGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCAc                                         <  2:9399134/93‑1 (MQ=255)
                                                  aCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCAccc                                       >  1:5655968/1‑93 (MQ=255)
                                                    atatAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCt                                     <  1:1430556/93‑1 (MQ=255)
                                                    atatAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCt                                     <  1:5486172/93‑1 (MQ=255)
                                                    atatAACCGATGTTGACATAAGCGCCGACGCTCCCCTGGGCGGCAACCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCt                                     <  1:8979920/93‑1 (MQ=255)
                                                        aaCCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGc                                 >  1:6794056/1‑93 (MQ=255)
                                                          ccGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCa                               <  2:1581152/93‑1 (MQ=255)
                                                          ccGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCa                               <  2:2700291/93‑1 (MQ=255)
                                                           cGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAg                              >  2:6044185/1‑93 (MQ=255)
                                                            gATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGc                             >  2:6875369/1‑93 (MQ=255)
                                                             aTGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCg                            >  1:8699054/1‑93 (MQ=255)
                                                              tGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGt                           >  1:9081736/1‑93 (MQ=255)
                                                               gTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGtt                          >  2:7643202/1‑93 (MQ=255)
                                                                ttGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTta                         >  1:8017034/1‑93 (MQ=255)
                                                                 tGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTtat                        <  1:7909430/93‑1 (MQ=255)
                                                                  gACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTtata                       <  1:7725846/93‑1 (MQ=255)
                                                                   aCGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATaa                      <  1:6326444/93‑1 (MQ=255)
                                                                      tAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATcc                   >  1:2073111/1‑93 (MQ=255)
                                                                      tAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATcc                   >  1:8330557/1‑93 (MQ=255)
                                                                      tAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATcc                   >  2:2490359/1‑93 (MQ=255)
                                                                       aaGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATccc                  >  1:2635967/1‑93 (MQ=255)
                                                                       aaGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATccc                  <  1:5241587/93‑1 (MQ=255)
                                                                        aGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCa                 <  2:9232632/93‑1 (MQ=255)
                                                                        aGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCa                 <  2:8270751/93‑1 (MQ=255)
                                                                        aGCGACGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCa                 >  2:6490754/1‑93 (MQ=255)
                                                                         gcgcCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCAt                <  1:310320/93‑1 (MQ=255)
                                                                          cgcCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATg               >  1:22415/1‑93 (MQ=255)
                                                                           gcCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGc              <  2:3989460/93‑1 (MQ=255)
                                                                            ccGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGCa             >  2:222046/1‑93 (MQ=255)
                                                                            ccGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGCa             >  1:7180744/1‑93 (MQ=255)
                                                                              gACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGCAAt           >  1:7068467/1‑93 (MQ=255)
                                                                                cGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGCAAttt         <  2:1683920/93‑1 (MQ=255)
                                                                                  cTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGCAATTTGa       >  1:4072658/1‑93 (MQ=255)
                                                                                  cTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGCAATTTGa       >  1:2835593/1‑93 (MQ=255)
                                                                                  cTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGCAATTTGa       >  1:499704/1‑93 (MQ=255)
                                                                                    ccccTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGCAATTTGaaa     <  2:928168/93‑1 (MQ=255)
                                                                                    ccccTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGCAATTTGaaa     >  1:8526697/1‑93 (MQ=255)
                                                                                     cccTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGCAATTTGAAAc    <  1:7340351/93‑1 (MQ=255)
                                                                                       cTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGCAATTTGAAACTg  >  1:2451700/1‑93 (MQ=255)
                                                                                       cTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGCAATTTGAAACTg  <  1:5889373/93‑1 (MQ=255)
                                                                                       cTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGCAATTTGAAACTg  >  1:5803936/1‑93 (MQ=255)
                                                                                       cTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGCAATTTGAAACTg  >  1:1887019/1‑93 (MQ=255)
                                                                                       cTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGCAATTTGAAACTg  >  1:1031707/1‑93 (MQ=255)
                                                                                          |                                                                                         
CCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGCAATTTGAAACTG  >  NC_000913/674973‑675152

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: