Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 961,423 A→G D69G (GAC→GGC)  rpsA → 30S ribosomal subunit protein S1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913961,4230AG99.5% 724.3 / ‑4.5 209D69G (GAC→GGC) rpsA30S ribosomal subunit protein S1
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (121/87);  total (121/88)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.03e-01

Only 100 of 216 total aligned reads displayed.

ACGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGACGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCACGAAGCCTGGATCACGCTGG  >  NC_000913/961333‑961512
                                                                                          |                                                                                         
aCGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCg                                                                                         >  2:4044805/1‑93 (MQ=255)
aCGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCg                                                                                         >  2:6046779/1‑93 (MQ=255)
aCGCTGGGCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGGAGTTGGCg                                                                                         >  2:3265228/1‑93 (MQ=255)
 cGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGt                                                                                        >  2:3412798/1‑93 (MQ=255)
 cGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGt                                                                                        <  1:3033118/93‑1 (MQ=255)
  gCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGtt                                                                                       >  2:4547893/1‑93 (MQ=255)
  gCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGtt                                                                                       <  1:3631323/93‑1 (MQ=255)
  ccTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGtt                                                                                       <  1:1418845/92‑1 (MQ=255)
      gTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGctc                                                                                   <  2:188829/93‑1 (MQ=255)
       tCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGctct                                                                                  <  1:5027259/93‑1 (MQ=255)
         tGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTgg                                                                                >  1:4903421/1‑93 (MQ=255)
          gAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGa                                                                               >  1:1000279/1‑93 (MQ=255)
           aaaTCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGGGACGAAGTTGGCGTTGCTCTGGAc                                                                              >  1:5156491/1‑93 (MQ=255)
            aaTCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACg                                                                             >  1:4545103/1‑93 (MQ=255)
             aTCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGc                                                                            <  1:5804841/93‑1 (MQ=255)
             aTCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGc                                                                            <  1:3531624/93‑1 (MQ=255)
             aTCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGc                                                                            >  1:3687867/1‑93 (MQ=255)
                 gAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTa                                                                        >  1:3900319/1‑93 (MQ=255)
                       gCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGAc                                                                  >  2:2094438/1‑93 (MQ=255)
                       gCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGAc                                                                  >  1:720163/1‑93 (MQ=255)
                           tCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCt                                                              <  1:254519/93‑1 (MQ=255)
                            cccGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCtt                                                             <  1:5553895/93‑1 (MQ=255)
                             ccGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTc                                                            <  2:2687205/93‑1 (MQ=255)
                              cGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCg                                                           <  1:6189558/93‑1 (MQ=255)
                              cGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGAAGGCTTCg                                                           >  2:4610874/1‑93 (MQ=255)
                                gCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGt                                                         <  1:1387997/93‑1 (MQ=255)
                                 cTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGGCGGCTTCGGTg                                                        <  1:1665977/93‑1 (MQ=255)
                                 cTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTg                                                        >  2:5614832/1‑93 (MQ=255)
                                  tGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGa                                                       >  2:2470486/1‑93 (MQ=255)
                                   gAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGaa                                                      >  1:3351028/1‑93 (MQ=255)
                                   gAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGaa                                                      >  2:2657544/1‑93 (MQ=255)
                                   gAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGaa                                                      >  2:5438257/1‑93 (MQ=255)
                                    aGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGaaa                                                     <  2:6299573/93‑1 (MQ=255)
                                    aGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGaaa                                                     >  1:5969693/1‑93 (MQ=255)
                                     gCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAAc                                                    >  1:2884830/1‑93 (MQ=255)
                                      cAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAAct                                                   <  1:3135734/93‑1 (MQ=255)
                                       aGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAActc                                                  >  1:690770/1‑93 (MQ=255)
                                       aGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAActc                                                  <  2:453563/93‑1 (MQ=255)
                                        gTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAActct                                                 >  1:5255188/1‑93 (MQ=255)
                                         ttCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTctg                                                <  1:1634945/93‑1 (MQ=255)
                                         ttCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTctg                                                <  2:109148/93‑1 (MQ=255)
                                          tCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTctgc                                               >  2:3532602/1‑93 (MQ=255)
                                           cAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTctgct                                              <  2:4993266/93‑1 (MQ=255)
                                            aaaaaCGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTctgctg                                             >  1:3104782/1‑93 (MQ=255)
                                               aaCGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTcc                                          >  1:970653/1‑93 (MQ=255)
                                                aCGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTccc                                         <  1:2191378/93‑1 (MQ=255)
                                                aCGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTccc                                         >  2:3760468/1‑93 (MQ=255)
                                                  gCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGt                                       >  1:1804940/1‑93 (MQ=255)
                                                   cccAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTg                                      <  2:5781549/93‑1 (MQ=255)
                                                   cccAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTg                                      <  2:6319146/93‑1 (MQ=255)
                                                    ccAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTga                                     <  1:2417317/93‑1 (MQ=255)
                                                       gggCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGaa                                  >  2:6075045/1‑93 (MQ=255)
                                                       gggCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGaa                                  >  2:1469618/1‑93 (MQ=255)
                                                       gggCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGaa                                  >  2:2975880/1‑93 (MQ=255)
                                                        ggCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGaaa                                 >  2:4900244/1‑93 (MQ=255)
                                                        ggCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGaaa                                 <  1:4826274/93‑1 (MQ=255)
                                                        ggCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGaaa                                 >  2:2613295/1‑93 (MQ=255)
                                                         gCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAg                                >  2:3844125/1‑93 (MQ=255)
                                                          cGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGc                               <  1:987635/93‑1 (MQ=255)
                                                          cGAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGc                               <  2:5651352/93‑1 (MQ=255)
                                                           gAGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCt                              <  2:3053169/93‑1 (MQ=255)
                                                            aGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTa                             >  2:2795051/1‑93 (MQ=255)
                                                            aGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTa                             >  1:1091855/1‑93 (MQ=255)
                                                            aGCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTa                             <  1:6412913/93‑1 (MQ=255)
                                                             gCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTaa                            >  2:5236858/1‑93 (MQ=255)
                                                             gCTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTaa                            >  2:3812382/1‑93 (MQ=255)
                                                              cTGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTaaa                           <  2:6192328/93‑1 (MQ=255)
                                                               tGGAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAAc                          <  2:6134358/93‑1 (MQ=255)
                                                                ggAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACg                         <  2:754501/93‑1 (MQ=255)
                                                                ggAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACg                         >  1:2347282/1‑93 (MQ=255)
                                                                ggAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACg                         >  1:5920622/1‑93 (MQ=255)
                                                                ggAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACg                         >  2:3910535/1‑93 (MQ=255)
                                                                 gAAATCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGt                        <  1:2025128/93‑1 (MQ=255)
                                                                   aaTCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCa                      >  2:1325969/1‑93 (MQ=255)
                                                                    aTCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCAc                     >  1:1510697/1‑93 (MQ=255)
                                                                    aTCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCAc                     <  2:50837/93‑1 (MQ=255)
                                                                    aTCCAGGTAGGTGACGAAGTTGACGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCAc                     <  2:2204269/93‑1 (MQ=255)
                                                                    aTCCAGGGAGGTGACGAAGTTGGCGTTGCTCTGGCCGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCAc                     >  1:4148617/1‑93 (MQ=255)
                                                                     tCCAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCACg                    >  2:2577978/1‑93 (MQ=255)
                                                                      ccAGGTAGGTGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCACGa                   >  1:4586933/1‑93 (MQ=255)
                                                                        aggtaggtGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCACGAAg                 <  2:4434623/93‑1 (MQ=255)
                                                                        aggtaggtGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCACGAAg                 >  1:774549/1‑93 (MQ=255)
                                                                         ggtaggtGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCACGAAGc                >  2:1833276/1‑93 (MQ=255)
                                                                          gtaggtGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCACGAAGcc               <  2:4168534/93‑1 (MQ=255)
                                                                          gtaggtGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCACGAAGcc               <  1:5540564/93‑1 (MQ=255)
                                                                          gtaggtGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCACGAAGcc               <  1:6259542/93‑1 (MQ=255)
                                                                            aggtGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCACGAAGCCTg             <  1:1237971/93‑1 (MQ=255)
                                                                            aggtGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCACGAAGCCTg             >  1:1843658/1‑93 (MQ=255)
                                                                             ggtGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTTTCCCGTGAGAAAGCTAAACGTCACGAAGCCTgg            <  2:1539859/93‑1 (MQ=255)
                                                                             ggtGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCACGAAGCCTgg            >  2:2052265/1‑93 (MQ=255)
                                                                             ggtGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCACGAAGCCTgg            >  1:5018641/1‑93 (MQ=255)
                                                                               tGACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCACGAAGCCTGGAt          <  1:273019/93‑1 (MQ=255)
                                                                                gACGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCACGAAGCCTGGATc         >  2:3617092/1‑93 (MQ=255)
                                                                                 aCGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCACGAAGCCTGGATCa        <  1:1463609/93‑1 (MQ=255)
                                                                                  cGAAGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCACGAAGCCTGGATCAc       <  1:3806002/93‑1 (MQ=255)
                                                                                    aaGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCACGAAGCCTGGATCACGc     >  2:1289464/1‑93 (MQ=255)
                                                                                    aaGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCACGAAGCCTGGATCACGc     >  1:4647716/1‑93 (MQ=255)
                                                                                    aaGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCACGAAGCCTGGATCACGc     >  1:1127419/1‑93 (MQ=255)
                                                                                     aGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCACGAAGCCTGGATCACGCt    >  1:6448065/1‑93 (MQ=255)
                                                                                     aGTTGGCGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCACGAAGCCTGGATCACGCt    >  2:1912584/1‑93 (MQ=255)
                                                                                          |                                                                                         
ACGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGACGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCACGAAGCCTGGATCACGCTGG  >  NC_000913/961333‑961512

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