Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,838,693 A→G I280T (ATC→ACC)  yicM ← putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,838,6930AG99.6% 1530.0 / ‑4.0 477I280T (ATC→ACC) yicMputative transport protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (209/265);  total (210/267)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

Only 100 of 492 total aligned reads displayed.

GTTAAGCCATCCACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGATGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGCAGTGAAGGC  >  NC_000913/3838603‑3838782
                                                                                          |                                                                                         
 ttAAGCCATCCACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTgg                                                                                        >  1:2426084/1‑93 (MQ=255)
    aGCCATCCACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGa                                                                                     >  1:4237517/1‑93 (MQ=255)
    aGCCATCCACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGa                                                                                     >  2:8468363/1‑93 (MQ=255)
      ccATCCACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGACGata                                                                                   >  1:7139800/1‑92 (MQ=255)
       cATCCACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGatca                                                                                  >  1:10018663/1‑93 (MQ=255)
        aTCCACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGatcat                                                                                 >  2:3850246/1‑93 (MQ=255)
         tCCACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGatcatc                                                                                >  1:1865801/1‑93 (MQ=255)
          ccACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATcc                                                                               <  1:2605038/93‑1 (MQ=255)
          ccACGCCGAAGCCCGCCAGGGCCAGATACACCGGGCGAATATACGAGAAGAAAGCAAACAGCCCGGCGAAAGACATGAAGGTGGCGAACATcc                                                                               >  1:8942717/1‑93 (MQ=255)
           cACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATccc                                                                              <  2:9741645/93‑1 (MQ=255)
            aCGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCACt                                                                             >  2:5424698/1‑93 (MQ=255)
              gCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGc                                                                           >  1:7001168/1‑93 (MQ=255)
               ccGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGcc                                                                          >  1:7841263/1‑93 (MQ=255)
                   aTCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATca                                                                      >  2:7781463/1‑93 (MQ=255)
                   aTCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATca                                                                      >  2:7956769/1‑93 (MQ=255)
                    tCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATcac                                                                     >  1:4336849/1‑93 (MQ=255)
                    tCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATACCTGCCATcac                                                                     >  1:8282389/1‑93 (MQ=255)
                       cGCCAGGTTCATATACCCCGGGCGAATATACGGGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACAcc                                                                  <  2:3299077/93‑1 (MQ=255)
                       cGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACAcc                                                                  <  2:5883583/93‑1 (MQ=255)
                        gacAGGTTCATATACACCGGGCGAAGATACTTGAAGAAAGCAAACTGCCCGGCAAAAGACATGAAGGTGGCGATCATCCCTGCCATCACAccc                                                                 <  2:5678555/91‑1 (MQ=255)
                        gCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACAccc                                                                 >  1:2102621/1‑93 (MQ=255)
                         ccAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCg                                                                <  1:6625909/93‑1 (MQ=255)
                          cAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCgg                                                               >  1:7265044/1‑93 (MQ=255)
                           aGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGa                                                              <  1:4976190/93‑1 (MQ=255)
                           aGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGa                                                              >  1:4896745/1‑93 (MQ=255)
                           aGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGa                                                              <  2:5148409/93‑1 (MQ=255)
                            ggTTCATATACACCGGGCGAATATCCGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACCCCCGGAc                                                             <  2:1387390/93‑1 (MQ=255)
                            ggTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGAc                                                             <  2:9565509/93‑1 (MQ=255)
                            ggTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGAc                                                             >  2:8560451/1‑93 (MQ=255)
                             gTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACg                                                            >  2:626107/1‑93 (MQ=255)
                               tCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGtt                                                          >  1:327810/1‑93 (MQ=255)
                               tCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGtt                                                          >  2:3977098/1‑93 (MQ=255)
                                 atatACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGt                                                        >  2:7267860/1‑93 (MQ=255)
                                  tataCACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGta                                                       >  2:4072651/1‑93 (MQ=255)
                                  tataCACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGta                                                       >  1:4721984/1‑93 (MQ=255)
                                   ataCCCCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGtaa                                                      <  2:5356436/93‑1 (MQ=255)
                                   ataCACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGtaa                                                      >  1:319016/1‑93 (MQ=255)
                                    taCACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGtaat                                                     <  2:3662691/93‑1 (MQ=255)
                                    taCACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGtaat                                                     >  1:6808396/1‑93 (MQ=255)
                                     acacCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGtaata                                                    >  1:5231323/1‑93 (MQ=255)
                                      cacCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGtaataa                                                   <  1:5921766/93‑1 (MQ=255)
                                      cacCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGtaataa                                                   <  1:3739209/93‑1 (MQ=255)
                                      cacCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGtaataa                                                   >  1:3555757/1‑93 (MQ=255)
                                      cacCGGGCGAAAATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGtaataa                                                   <  1:8473501/93‑1 (MQ=255)
                                        ccGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGc                                                 <  2:4723467/93‑1 (MQ=255)
                                        ccGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGc                                                 >  1:1674265/1‑93 (MQ=255)
                                        ccGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGc                                                 <  2:6854720/93‑1 (MQ=255)
                                         cGGGCGAATATACGTGAAGAAAGCAAGCTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCg                                                <  1:1350958/93‑1 (MQ=255)
                                         cGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCg                                                >  1:2939978/1‑93 (MQ=255)
                                         cGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCg                                                <  1:3363189/93‑1 (MQ=255)
                                         cGGGCGAATAGACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCg                                                <  2:3550876/93‑1 (MQ=255)
                                          gggCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCgg                                               >  1:7680014/1‑93 (MQ=255)
                                            gCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGaa                                             >  1:8233513/1‑93 (MQ=255)
                                            gCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGaa                                             <  2:1501143/93‑1 (MQ=255)
                                             cGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGaaa                                            <  1:4562112/93‑1 (MQ=255)
                                             cGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACACGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGaaa                                            <  2:10006795/93‑1 (MQ=255)
                                               aaTATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGt                                          >  2:3124199/1‑93 (MQ=255)
                                                atatACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTa                                         <  2:603038/93‑1 (MQ=255)
                                                atatACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTa                                         <  2:7140015/93‑1 (MQ=255)
                                                 tataCGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTAt                                        >  2:6553946/1‑93 (MQ=255)
                                                 tataCGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTAt                                        >  1:6284075/1‑93 (MQ=255)
                                                  ataCGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTAtt                                       <  2:4938969/93‑1 (MQ=255)
                                                   taCGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTAttt                                      >  2:1315322/1‑93 (MQ=255)
                                                    aCGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTAtttt                                     <  2:7615149/93‑1 (MQ=255)
                                                    aCGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTAtttt                                     <  1:9570981/93‑1 (MQ=255)
                                                     cGTGAAGAAAGCAACCTGCCGGGCGAAGGCCA‑GGAGGTGGCGATCATCCCTGCCGTCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGt                                   <  2:9267290/93‑1 (MQ=255)
                                                     cGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTg                                    <  1:1174006/93‑1 (MQ=255)
                                                      gTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGt                                   >  2:2138953/1‑93 (MQ=255)
                                                      gTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGt                                   <  1:548569/93‑1 (MQ=255)
                                                      gTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGt                                   <  1:9865102/93‑1 (MQ=255)
                                                       tGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGtt                                  <  2:189551/93‑1 (MQ=255)
                                                       tGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGtt                                  >  1:9275090/1‑93 (MQ=255)
                                                         aagaaAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTc                                <  1:3900080/93‑1 (MQ=255)
                                                            aaaGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGa                             <  2:8111250/93‑1 (MQ=255)
                                                            aaaGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGa                             <  1:8853084/93‑1 (MQ=255)
                                                            aaaGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGa                             <  1:8617192/93‑1 (MQ=255)
                                                             aaGCAAACTGCCCGGCGCAGGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTATCTGAt                            <  2:8983811/93‑1 (MQ=255)
                                                             aaGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGAt                            <  1:5788598/93‑1 (MQ=255)
                                                               gCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGc                          >  1:6384685/1‑93 (MQ=255)
                                                                 aaaCTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGa                        <  2:5569011/93‑1 (MQ=255)
                                                                 aaaCTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGa                        >  2:5625851/1‑93 (MQ=255)
                                                                  aaCTGCCGGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAg                       <  2:2518664/93‑1 (MQ=255)
                                                                  aaCTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAg                       <  1:8117195/93‑1 (MQ=255)
                                                                   aCTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAgg                      <  1:7749552/93‑1 (MQ=255)
                                                                   aCTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAgg                      <  2:4293476/93‑1 (MQ=255)
                                                                      gaccGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGtt                   <  2:1689748/91‑1 (MQ=255)
                                                                      gCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGtt                   <  1:1454329/93‑1 (MQ=255)
                                                                       cccGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTc                  >  1:6137186/1‑93 (MQ=255)
                                                                        ccGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCg                 >  2:2728865/1‑93 (MQ=255)
                                                                         cGGGGAAAGACATTTAGGTGGCGATCCTCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGCGTTTTGTTTCTTATGCGAGGGCTCGc                >  2:8301006/1‑93 (MQ=255)
                                                                         cGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGc                <  2:225148/93‑1 (MQ=255)
                                                                          ggCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGcc               >  2:6459370/1‑93 (MQ=255)
                                                                          ggCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGcc               <  1:9685324/93‑1 (MQ=255)
                                                                            cGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTg             <  1:1005516/93‑1 (MQ=255)
                                                                             gAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTgg            >  2:936936/1‑93 (MQ=255)
                                                                              aaaGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGc           >  2:8761555/1‑93 (MQ=255)
                                                                                aGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGCAg         >  1:107763/1‑93 (MQ=255)
                                                                                 gACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGCAGt        <  1:561722/93‑1 (MQ=255)
                                                                                   cATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGCAGTGa      <  1:9420487/93‑1 (MQ=255)
                                                                                       aaGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGCAGTGAAGGc  <  1:7505112/93‑1 (MQ=255)
                                                                                          |                                                                                         
GTTAAGCCATCCACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGATGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGCAGTGAAGGC  >  NC_000913/3838603‑3838782

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