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.7% 1010.4 / ‑4.8 309I280T (ATC→ACC) yicMputative transport protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (126/182);  total (126/183)
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.72e-01

Only 100 of 315 total aligned reads displayed.

GTTAAGCCATCCACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGATGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGCAGTGAAGGC  >  NC_000913/3838603‑3838782
                                                                                          |                                                                                         
gTTAAGCCATCCACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTg                                                                                         <  1:4451500/93‑1 (MQ=255)
 ttAAGCCATCCACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTgg                                                                                        >  1:529262/1‑93 (MQ=255)
 ttAAGCCATCCACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTgg                                                                                        >  1:2469769/1‑93 (MQ=255)
    aGCCATCCACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGa                                                                                     <  2:5354330/93‑1 (MQ=255)
    aGCCATCCACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGa                                                                                     <  2:712558/93‑1 (MQ=255)
     gCCATCCACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGat                                                                                    <  1:4538870/93‑1 (MQ=255)
      ccATCCACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGatc                                                                                   <  1:5919236/93‑1 (MQ=255)
      ccATCCACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGatc                                                                                   <  2:4262161/93‑1 (MQ=255)
       cATCCACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGatca                                                                                  <  2:3193690/93‑1 (MQ=255)
       cATCCACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGatca                                                                                  <  1:1358838/93‑1 (MQ=255)
          ccACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATcc                                                                               >  2:4858565/1‑93 (MQ=255)
          ccACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATcc                                                                               >  1:3067168/1‑93 (MQ=255)
           cACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATccc                                                                              <  1:6144430/93‑1 (MQ=255)
           cACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATccc                                                                              <  1:4292640/93‑1 (MQ=255)
           cACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATccc                                                                              <  2:6193015/93‑1 (MQ=255)
            aCGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCt                                                                             <  2:4653048/93‑1 (MQ=255)
             cGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTg                                                                            >  1:4172705/1‑93 (MQ=255)
                cGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCa                                                                         <  1:5762241/93‑1 (MQ=255)
                  aaTCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATc                                                                       >  1:5806556/1‑93 (MQ=255)
                  aaTCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATc                                                                       >  2:3425335/1‑93 (MQ=255)
                  aaTCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATc                                                                       <  2:2942064/93‑1 (MQ=255)
                  aaTCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATc                                                                       >  2:3004420/1‑93 (MQ=255)
                  aaTCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATc                                                                       >  2:5679722/1‑93 (MQ=255)
                   aTCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATca                                                                      >  2:2356792/1‑93 (MQ=255)
                    tCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATcac                                                                     <  1:468736/93‑1 (MQ=255)
                     cccGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATcaca                                                                    <  1:5215277/93‑1 (MQ=255)
                        gCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACAccc                                                                 <  1:4560813/93‑1 (MQ=255)
                         ccAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCg                                                                >  1:4946704/1‑93 (MQ=255)
                          cAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCgg                                                               >  1:5604275/1‑93 (MQ=255)
                           aGGTTCATCTACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGa                                                              >  1:6543142/1‑93 (MQ=255)
                           aGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGa                                                              >  2:1745951/1‑93 (MQ=255)
                            ggTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGAc                                                             >  2:5164589/1‑93 (MQ=255)
                             gTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACg                                                            >  2:645514/1‑93 (MQ=255)
                              gtCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGt                                                           <  2:3975529/92‑1 (MQ=255)
                               tCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGtt                                                          >  1:878397/1‑93 (MQ=255)
                                cATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTg                                                         >  2:3656916/1‑93 (MQ=255)
                                 atatACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGt                                                        >  2:5889756/1‑93 (MQ=255)
                                  tataCACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGta                                                       <  1:1505835/93‑1 (MQ=255)
                                   ataCACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGtaa                                                      >  2:2759255/1‑93 (MQ=255)
                                    taCACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGtaat                                                     <  1:2759360/93‑1 (MQ=255)
                                    taCACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGtaat                                                     >  2:1881921/1‑93 (MQ=255)
                                      cacCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGtaataa                                                   <  2:1508375/93‑1 (MQ=255)
                                        ccGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGc                                                 >  2:3501174/1‑93 (MQ=255)
                                        ccGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGc                                                 <  2:6117821/93‑1 (MQ=255)
                                          gggCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCgg                                               <  2:922234/93‑1 (MQ=255)
                                             cGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGaaa                                            <  2:1169080/93‑1 (MQ=255)
                                              gAATATACGTGTAGAAAGCAACCTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAg                                           <  1:2323092/93‑1 (MQ=255)
                                              gAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAg                                           >  2:1401030/1‑93 (MQ=255)
                                                 tataCGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTAt                                        <  1:2667617/93‑1 (MQ=255)
                                                  ataCGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTAtt                                       >  1:1239206/1‑93 (MQ=255)
                                                  ataCGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTAtt                                       >  1:5958672/1‑93 (MQ=255)
                                                  ataCGCGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTAtt                                       >  1:3545263/1‑93 (MQ=255)
                                                   taCGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTATTAAGCGGAAAGTAttt                                      <  2:1989651/93‑1 (MQ=255)
                                                    aCGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTAtttt                                     <  2:815494/93‑1 (MQ=255)
                                                      gTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGt                                   <  1:1020389/93‑1 (MQ=255)
                                                      gTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGt                                   <  2:5284186/93‑1 (MQ=255)
                                                       tGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAACGTATTTTGtt                                  >  1:96657/1‑93 (MQ=255)
                                                       ggaagaaAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACCCCCGGACGTTGTAATAAGCGGAAAGTATTTTGtt                                  <  1:6339866/92‑1 (MQ=255)
                                                        gaagaaAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGttt                                 <  1:3604877/93‑1 (MQ=255)
                                                         aagaaAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTc                                <  1:192812/93‑1 (MQ=255)
                                                         aagaaAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTc                                <  1:2880269/93‑1 (MQ=255)
                                                         aagaaAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTc                                >  1:3431578/1‑93 (MQ=255)
                                                          agaaAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCt                               <  1:4737940/93‑1 (MQ=255)
                                                           gaaAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTg                              <  2:2605051/93‑1 (MQ=255)
                                                           gaaAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTg                              <  1:1174799/93‑1 (MQ=255)
                                                            aaaGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGa                             <  1:1381088/93‑1 (MQ=255)
                                                             aaGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGAt                            <  1:1060963/93‑1 (MQ=255)
                                                             aaGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGAt                            >  1:5929639/1‑93 (MQ=255)
                                                             aaGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGAt                            <  2:2082041/93‑1 (MQ=255)
                                                              aGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATg                           <  1:1127186/93‑1 (MQ=255)
                                                                cAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCg                         <  1:318294/93‑1 (MQ=255)
                                                                cAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCg                         >  1:5390756/1‑93 (MQ=255)
                                                                 aaaCTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACTCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGa                        <  2:3570982/93‑1 (MQ=255)
                                                                 aaaCTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGa                        <  1:2819963/93‑1 (MQ=255)
                                                                 aaaCTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGa                        >  2:1568062/1‑93 (MQ=255)
                                                                  aaCTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAg                       >  2:5020301/1‑93 (MQ=255)
                                                                   aCTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAgg                      >  2:6398263/1‑93 (MQ=255)
                                                                    cTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAggg                     <  2:3318954/93‑1 (MQ=255)
                                                                     tGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGt                    <  2:4502728/93‑1 (MQ=255)
                                                                       cccGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTc                  <  2:6616314/93‑1 (MQ=255)
                                                                       cccGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTc                  <  2:5795761/93‑1 (MQ=255)
                                                                        ccGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCg                 >  2:5453984/1‑93 (MQ=255)
                                                                          ggCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGcc               <  2:2670205/93‑1 (MQ=255)
                                                                           gCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCt              >  2:4754701/1‑93 (MQ=255)
                                                                           gCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCt              >  2:6063903/1‑93 (MQ=255)
                                                                            cGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTg             <  2:595050/93‑1 (MQ=255)
                                                                             gAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTgg            <  1:6063330/93‑1 (MQ=255)
                                                                             gAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTgg            >  1:2177762/1‑93 (MQ=255)
                                                                             gAAAGACATGAAGGTGGCGATCATCCCTGCCAGCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTgg            <  2:2833811/93‑1 (MQ=255)
                                                                               aaGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGCa          >  2:3838127/1‑93 (MQ=255)
                                                                                aGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGCAg         >  1:1303592/1‑93 (MQ=255)
                                                                                aGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGCAg         >  1:3849619/1‑93 (MQ=255)
                                                                                  aCATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGCATTTTGTTTCTGATGCGAGGGTTCGCCTGGCAGta       >  2:4213336/1‑92 (MQ=255)
                                                                                   cATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGCAGTGa      >  2:1937726/1‑93 (MQ=255)
                                                                                    aTGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGCAGTGaa     <  2:3739600/93‑1 (MQ=255)
                                                                                    aTGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGCAGTGaa     >  1:1631051/1‑93 (MQ=255)
                                                                                    aTGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGCAGTGaa     <  1:4886043/93‑1 (MQ=255)
                                                                                    aTGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGCAGTGaa     <  2:5118071/93‑1 (MQ=255)
                                                                                     tGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGCAGTGAAg    >  1:408547/1‑93 (MQ=255)
                                                                                       aaGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGCAGTGAAGGc  >  2:1283821/1‑93 (MQ=255)
                                                                                          |                                                                                         
GTTAAGCCATCCACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGATGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGCAGTGAAGGC  >  NC_000913/3838603‑3838782

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