Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,703,251 A→T 16.0% intergenic (+17/+17) add → / ← ydgJ adenosine deaminase/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,703,2510AT16.0% 333.4 / 51.5 131intergenic (+17/+17)add/ydgJadenosine deaminase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (62/48);  new base T (6/15);  total (68/63)
Fisher's exact test for biased strand distribution p-value = 3.04e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.72e-01

Only 100 of 132 total aligned reads displayed.

CACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTG  >  NC_000913/1703155‑1703350
                                                                                                |                                                                                                   
cacaGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATg                                                                                                  <  2:383955/100‑1 (MQ=255)
   aGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCAc                                                                                               <  1:3334583/100‑1 (MQ=255)
    gATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCAc                                                                                               >  1:1972055/1‑99 (MQ=255)
       tAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCATTTGCAccatc                                                                                           <  1:29067/100‑1 (MQ=255)
        aaTGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCATTTGCAccatcc                                                                                          <  2:2991655/100‑1 (MQ=255)
        aaTGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCATTTGCAccatcc                                                                                          <  1:2269842/100‑1 (MQ=255)
        aaTGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCATTTGCAccatcc                                                                                          <  1:3459011/100‑1 (MQ=255)
        aaTGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCAccatcc                                                                                          >  1:3049368/1‑100 (MQ=255)
        aaTGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCAccatcc                                                                                          >  1:3499296/1‑100 (MQ=255)
        aaTGGTATACAAATGGCTTTCCTAAGAGCTGAAGAAAAACGAAAAATAAAAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCATTTGCAccatcc                                                                                          <  2:1272247/100‑1 (MQ=16)
         aTGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAGAATGGATGGTGCATTTGCAccatcca                                                                                         <  1:3397941/100‑1 (MQ=255)
         aTGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCATTTGCAccatcca                                                                                         <  2:272399/100‑1 (MQ=255)
          tGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCATTTGCAccatccat                                                                                        <  2:741040/100‑1 (MQ=255)
          tGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCATTTGCAccatccat                                                                                        <  1:424423/100‑1 (MQ=255)
          tGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCATTTGCAccatccat                                                                                        <  1:2727930/100‑1 (MQ=255)
          tGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCAccatccat                                                                                        <  2:698101/100‑1 (MQ=255)
           ggTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCAtt                                                                                       >  1:2016595/1‑100 (MQ=255)
           ggTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCAtt                                                                                       >  1:1441268/1‑100 (MQ=255)
           ggTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCAtt                                                                                       >  2:3074283/1‑100 (MQ=255)
           ggTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCAtt                                                                                       <  1:1363908/100‑1 (MQ=255)
            gTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCAttt                                                                                      >  2:2326864/1‑100 (MQ=255)
               tgaaaaTTGCTTTCCTCAGCGCTGATGAAAAACGCTCACTGCAAGAAAAAGTCGCCGCGAAGTAACAAAATGGATAGTGCATTTTCACCATCCAttttg                                                                                    <  2:1323352/96‑2 (MQ=25)
                ggAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCa                                                                                  >  1:985863/1‑100 (MQ=255)
                ggAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCa                                                                                  <  2:1951864/100‑1 (MQ=255)
                 gAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCAt                                                                                 <  2:816496/100‑1 (MQ=255)
                   cgcgtttttCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGc                                                                               <  1:2753684/94‑1 (MQ=255)
                   aaTGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGc                                                                               >  2:1995651/1‑100 (MQ=255)
                     cgtttttCCTCAGCGCTGAGGAAAGACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCaa                                                                             >  2:3490382/5‑100 (MQ=255)
                     cgtttttCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCaa                                                                             >  1:2645758/5‑100 (MQ=255)
                      gtttttCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAg                                                                            >  1:1286341/4‑100 (MQ=255)
                       gCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCGATTTTTCATGCAAgg                                                                           >  1:3889831/1‑100 (MQ=255)
                        cTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGc                                                                          >  1:260546/1‑100 (MQ=255)
                        cTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGc                                                                          >  2:124811/1‑100 (MQ=255)
                            ccTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACaa                                                                      <  2:1398396/100‑1 (MQ=255)
                             cTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACaaa                                                                     >  1:350911/1‑100 (MQ=255)
                                 gcgcTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTcg                                                                 >  2:292833/1‑100 (MQ=255)
                                   gcTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTcgcg                                                               <  2:68661/100‑1 (MQ=255)
                                    cTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTcgcgc                                                              <  2:350911/100‑1 (MQ=255)
                                     tGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTcgcgcg                                                             <  1:1236914/100‑1 (MQ=255)
                                     tGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTcgcgcg                                                             <  1:3240472/100‑1 (MQ=255)
                                      gAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGa                                                            >  1:3416859/1‑100 (MQ=255)
                                      gAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGAGGGTGCAAATGAACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGa                                                            <  2:3476328/100‑1 (MQ=255)
                                       aGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGAt                                                           >  1:1513871/1‑100 (MQ=255)
                                        ggAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATg                                                          >  2:2870771/1‑100 (MQ=255)
                                         gAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGt                                                         >  2:1118423/1‑100 (MQ=255)
                                          aaaaaCGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGtt                                                        <  1:2848511/100‑1 (MQ=255)
                                             aaCGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTgg                                                     >  1:650472/1‑100 (MQ=255)
                                              aCGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGTACCATCCATTTTTCATGCAAGGCGCAAAGTCGCTCGATGTTTg                                                      >  2:1707169/1‑98 (MQ=255)
                                               cgcgCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCg                                                   >  2:493815/1‑100 (MQ=255)
                                               cgcgCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGAGCAATGTTTGGCg                                                   >  1:2069260/1‑100 (MQ=255)
                                                   cACTGCGAGAAAAAGTCGCCGCGAATTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGAtt                                               <  1:599220/100‑1 (MQ=255)
                                                    aCTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTc                                              >  1:2739353/1‑100 (MQ=255)
                                                    aCTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTc                                              <  1:2154178/100‑1 (MQ=255)
                                                     cTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCg                                             >  2:3840709/1‑100 (MQ=255)
                                                       gCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGa                                            <  2:2068105/99‑1 (MQ=255)
                                                        cGAGAAAAAGTCGCCGCGAAGTAACAAAATGGCTGGTGCAAATGCCCCATCCATTTTTCATGCATGGCACCAAGTCGCGCGATGTTTGGCGGATTCGATg                                          >  1:124717/1‑100 (MQ=255)
                                                        cGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATg                                          >  1:3305987/1‑100 (MQ=255)
                                                        cGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATg                                          >  2:3834204/1‑100 (MQ=255)
                                                             aaaaGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAg                                     >  1:3784937/1‑100 (MQ=255)
                                                              aaaGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGc                                    <  1:2318981/100‑1 (MQ=255)
                                                              aaaGTCGCCGCGAAATAACGAAATAGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGc                                    <  2:2391740/100‑1 (MQ=255)
                                                                aGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCtc                                  <  1:313677/100‑1 (MQ=255)
                                                                 gTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCtca                                 >  1:157467/1‑100 (MQ=255)
                                                                 gTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCtca                                 >  1:701594/1‑100 (MQ=255)
                                                                    gccgcGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCtcaatc                              <  1:3527014/100‑1 (MQ=255)
                                                                    gccgcGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCtcaatc                              >  1:346450/1‑100 (MQ=255)
                                                                    gccgcGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCtcaatc                              <  1:2512622/100‑1 (MQ=255)
                                                                    gccgcGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCtcaatc                              <  1:242934/100‑1 (MQ=255)
                                                                     ccgcGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCtcaatca                             <  2:2734734/100‑1 (MQ=255)
                                                                      cgcgAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCtcaatcaa                            >  1:2645399/1‑100 (MQ=255)
                                                                          aaGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTcc                        <  2:3821593/100‑1 (MQ=255)
                                                                          aaGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTcc                        <  1:3878505/100‑1 (MQ=255)
                                                                           gcatgaaaaaTGGATGGTGCATTTGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCa                       >  1:3735400/7‑100 (MQ=37)
                                                                           aGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCa                       <  2:1551338/100‑1 (MQ=255)
                                                                            gTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTccgt                      >  1:557793/1‑98 (MQ=255)
                                                                            gTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCCATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCAt                      >  2:2775532/1‑100 (MQ=255)
                                                                            catgaaaaaTGGATGGTGCATTTGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCAt                      >  1:1291104/6‑100 (MQ=37)
                                                                             tAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCAtt                     >  2:847101/1‑100 (MQ=255)
                                                                              tgaaaaaTGGATGGTGCATTTGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTa                    >  2:3294752/4‑100 (MQ=38)
                                                                              aaCAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTa                    >  2:1419600/1‑100 (MQ=255)
                                                                              aaCAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTa                    >  2:3616082/1‑100 (MQ=255)
                                                                                cAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTagc                  >  2:2395367/1‑98 (MQ=255)
                                                                                cAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTAcc                  >  1:1578318/1‑100 (MQ=255)
                                                                                cAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTAcc                  >  1:1614222/1‑100 (MQ=255)
                                                                                 aaaaTGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCt                 <  2:3148863/100‑1 (MQ=255)
                                                                                   aaTGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTgg               <  1:2950969/100‑1 (MQ=255)
                                                                                    atggatggTGCATTTGCACCATCCATTTTTTTTTTTTTTTACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTg                >  2:1717600/1‑98 (MQ=21)
                                                                                    atggatggTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGa              <  2:2860147/100‑1 (MQ=255)
                                                                                      ggatggTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGAtt            >  1:25307/1‑100 (MQ=255)
                                                                                      ggatggTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGAtt            >  1:1433313/1‑100 (MQ=255)
                                                                                      ggatggTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGAtt            >  2:3335449/1‑100 (MQ=255)
                                                                                       gatggTGCATTTGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGAtt            >  2:3721032/1‑99 (MQ=255)
                                                                                       gatggTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTg           <  1:806246/100‑1 (MQ=255)
                                                                                          ggTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCt        <  1:2099063/100‑1 (MQ=255)
                                                                                          ggTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCt        <  2:3411488/100‑1 (MQ=255)
                                                                                            tGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTgg      >  1:652402/1‑100 (MQ=255)
                                                                                            tGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTgg      <  2:2695224/100‑1 (MQ=255)
                                                                                             gCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGc     <  1:1822659/100‑1 (MQ=255)
                                                                                              cAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCt    >  1:2817948/1‑100 (MQ=255)
                                                                                                aaTGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTg  >  1:366390/1‑100 (MQ=255)
                                                                                                |                                                                                                   
CACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTG  >  NC_000913/1703155‑1703350

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