Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,187,091 T→A 17.1% intergenic (+21/+28) pepT → / ← roxA peptidase T/50S ribosomal protein L16 arginine hydroxylase; 2‑oxoglutarate oxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,187,0910TA17.1% 405.0 / 75.7 171intergenic (+21/+28)pepT/roxApeptidase T/50S ribosomal protein L16 arginine hydroxylase; 2‑oxoglutarate oxygenase
Reads supporting (aligned to +/- strand):  ref base T (76/65);  new base A (12/17);  total (89/82)
Fisher's exact test for biased strand distribution p-value = 2.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.59e-01

Only 100 of 178 total aligned reads displayed.

TTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAG  >  NC_000913/1186993‑1187190
                                                                                                  |                                                                                                   
tttGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCAtgtg                                                                                                    <  1:2744580/100‑1 (MQ=255)
 ttGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGa                                                                                                   >  1:113602/1‑100 (MQ=255)
    tgACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGc                                                                                                <  1:2923292/100‑1 (MQ=255)
    tgACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGc                                                                                                >  1:2093532/1‑100 (MQ=255)
    tgACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGc                                                                                                >  1:2073419/1‑100 (MQ=255)
     gACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGcc                                                                                               >  1:1900633/1‑100 (MQ=255)
      aCTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCg                                                                                              <  2:2653525/100‑1 (MQ=255)
       ctctGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGc                                                                                             >  1:631267/1‑100 (MQ=255)
       ctctGGAAGGGATGGAAAAAGCGGTGCAGGTGATCGGCCGGATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGcc                                                                                               <  1:1113751/98‑1 (MQ=255)
         ctGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCAt                                                                                           >  1:1561452/1‑100 (MQ=255)
          tGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGCGATGCCGCATc                                                                                          >  2:3122560/1‑100 (MQ=255)
          tGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCATc                                                                                          <  1:500685/100‑1 (MQ=255)
          tGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCATc                                                                                          <  1:3029733/100‑1 (MQ=255)
           ggAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATcc                                                                                         >  1:677505/1‑100 (MQ=255)
           ggAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATcc                                                                                         >  1:2259679/1‑100 (MQ=255)
           ggAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATcc                                                                                         <  1:645418/100‑1 (MQ=255)
           ggAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATc                                                                                          <  2:220162/99‑1 (MQ=255)
           ggAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCATcc                                                                                         <  1:2336962/100‑1 (MQ=255)
           ggAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCATcc                                                                                         <  1:1858408/100‑1 (MQ=255)
            gAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCATccc                                                                                        <  1:14486/100‑2 (MQ=39)
            gAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCATccc                                                                                        <  1:2698258/100‑2 (MQ=39)
            gAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCATccc                                                                                        <  1:2875839/100‑2 (MQ=39)
             aaGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCATcc                                                                                         <  2:3553861/98‑1 (MQ=39)
                 tATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCTAAAAGGTATGCGGCATGTGATGCCGCATCCGGCtt                                                                                    >  2:155758/1‑99 (MQ=255)
                  aTGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTaa                                                                                  >  2:1452360/1‑100 (MQ=255)
                  aTGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTaa                                                                                  <  1:3140650/100‑1 (MQ=255)
                  aTGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCATcccttttcg                                                                                  <  1:914063/100‑8 (MQ=37)
                    ggAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAAt                                                                                >  2:1893318/1‑100 (MQ=255)
                    ggAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGAAAGAAAGAGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAAt                                                                                <  2:569548/100‑1 (MQ=255)
                     gAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATc                                                                               >  2:2795153/1‑100 (MQ=255)
                       aaaaGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCa                                                                             <  2:1117161/100‑1 (MQ=255)
                        aaaGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAAGCCaa                                                                            <  2:313473/100‑1 (MQ=255)
                         aaGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCaaa                                                                           <  2:2893735/100‑1 (MQ=255)
                           gCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACt                                                                         >  2:659207/1‑100 (MQ=255)
                           gCGGTGCAGGTGATCGTCCGTATTGCAGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATGCGGCGTAAATCCAAACt                                                                         >  1:2221440/1‑100 (MQ=255)
                             ggTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTa                                                                       >  1:1927665/1‑100 (MQ=255)
                               tGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTAcc                                                                     >  2:3587416/1‑100 (MQ=255)
                               tGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTAcc                                                                     >  2:1287787/1‑100 (MQ=255)
                               tGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGAATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTAcc                                                                     >  2:2332159/1‑100 (MQ=255)
                                gCAGGTGATCGTCAGTATTGCGGAGTTAACGGCGCAAAAAAAAAAAAAAAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACt                                                                         <  2:965713/95‑1 (MQ=255)
                                  aGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCtt                                                                  >  2:2483383/1‑100 (MQ=255)
                                   ggTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTc                                                                 <  1:484291/100‑1 (MQ=255)
                                      gATCGTCCGTATTGCCGAGTTAACGGCCCAACGGCAGTAAGCGAACACGGCTGCGGCATGTGATCCCGCATCCGACCTCACTCCACACTTAccc                                                                    >  2:3643236/1‑94 (MQ=255)
                                           tCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGCAGAAcc                                                         >  2:420224/1‑100 (MQ=255)
                                           tCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAAcc                                                         <  2:2423302/100‑1 (MQ=255)
                                           tCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAAcc                                                         <  2:2991564/100‑1 (MQ=255)
                                            ccGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAActa                                                        <  1:2467029/100‑3 (MQ=255)
                                            ccGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCa                                                        <  2:13767/100‑1 (MQ=255)
                                              gTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAAt                                                      >  1:3543218/1‑100 (MQ=255)
                                                  tGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCCAACTTACCCTTCGAAGAACCAATAccc                                                  >  1:2830172/1‑100 (MQ=255)
                                                   gCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCg                                                 >  2:1752501/1‑100 (MQ=255)
                                                   gCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCg                                                 >  1:758324/1‑100 (MQ=255)
                                                   gCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCg                                                 >  1:1699365/1‑100 (MQ=255)
                                                    ccGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGc                                                <  1:898120/100‑1 (MQ=255)
                                                    ccGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGc                                                >  2:3355436/1‑100 (MQ=255)
                                                    ccGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGc                                                <  1:281977/100‑1 (MQ=255)
                                                      gAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTa                                              >  1:3301874/1‑100 (MQ=255)
                                                       aGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTAt                                             >  1:1329819/1‑100 (MQ=255)
                                                         ttAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTg                                           <  2:3369573/100‑1 (MQ=255)
                                                           aaCGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGAc                                         >  1:2017762/1‑100 (MQ=255)
                                                           aaCGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGAc                                         >  2:1526704/1‑100 (MQ=255)
                                                            aCGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGAcc                                        <  1:1081413/100‑1 (MQ=255)
                                                             cGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCa                                       <  1:2117498/100‑1 (MQ=255)
                                                             cGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACGTACCCTTCGAAGAACCAATACCCGCTATTGACCa                                       >  2:3162348/1‑100 (MQ=255)
                                                              ggCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAg                                      >  1:2845954/1‑100 (MQ=255)
                                                                cgcAACGGAATGAAGCGAAAAGGGATGCTGCATGTCACGCCGCATCCGGCTTAAATCCAAACTTACTCTTCGAACAACCAATACCCGCTAGTGAcc                                        <  2:3355907/96‑1 (MQ=255)
                                                                 gcAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAgcgc                                   >  2:507665/1‑100 (MQ=255)
                                                                     cGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCcgcg                               <  1:855675/100‑1 (MQ=255)
                                                                      ggAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGa                              <  2:2635494/100‑1 (MQ=255)
                                                                      ggAAGTAAGCGAAAAAGGATGCGGCATGTGATGCCGCATCCGACTTAAATCCAAACTTACACTTCAAAGAACCAATACCCGCTATTGACCACcaccgcaa                              >  2:1799232/1‑98 (MQ=255)
                                                                       gAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAg                             >  1:2503199/1‑100 (MQ=255)
                                                                       gAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCCAAGAACCAATACCCGCTATTGACCAGCGCCGCGAg                             <  1:343014/100‑1 (MQ=255)
                                                                       gAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAAGCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAg                             >  2:1753128/1‑100 (MQ=255)
                                                                           tAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATc                         >  1:173256/1‑100 (MQ=255)
                                                                           tAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATc                         >  2:3234263/1‑100 (MQ=255)
                                                                             aGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATcgc                       <  1:3619089/100‑1 (MQ=255)
                                                                                gAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAg                    >  1:213382/1‑100 (MQ=255)
                                                                                gAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAg                    >  1:1806278/1‑100 (MQ=255)
                                                                                 aaaaGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAgg                   <  1:718320/100‑1 (MQ=255)
                                                                                 aaaaGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACGTACCCTGCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAgg                   <  2:208687/100‑1 (MQ=255)
                                                                                   gccggATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGaa                 >  1:2243189/4‑100 (MQ=255)
                                                                                    ccggATGCGGCATCACATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAAt                >  1:1853568/3‑100 (MQ=38)
                                                                                     cggATGCGGCATCACATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATg               >  2:555396/2‑100 (MQ=39)
                                                                                     cggATCCGGCATCACATCCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGCCCCAATCCCCCCTATTGACCAGCGCCGCGAGCATCGCGAGGaa                 >  2:2343644/2‑98 (MQ=255)
                                                                                      ggATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGa              >  2:2116052/1‑100 (MQ=255)
                                                                                      ggATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGa              >  2:1575137/1‑100 (MQ=255)
                                                                                      ggATGCGGCATCACATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAAtaa              >  2:1227048/1‑98 (MQ=255)
                                                                                       gATGCGGCATGTGATGCAGCATCTGGCTTACATCCAAGCTTACGCTTCGTAGAACCAATAGCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGa              >  2:2381536/1‑99 (MQ=255)
                                                                                        aTGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACg            <  1:3355461/100‑1 (MQ=255)
                                                                                        aTGCGGCATGTGATGCCGCACCCGGCTTAAATCCAAGCTTACCCTTCGAAGAACCAATACCCGCGATTGACCAGCGCCGCGAGCATCGCGAGGAATGACg            <  1:2553432/100‑1 (MQ=255)
                                                                                         tGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACgg           >  1:20389/1‑100 (MQ=255)
                                                                                         tGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACgg           >  2:274890/1‑100 (MQ=255)
                                                                                            ggCATGTGATTCCGCATCCGGCATAAATCCAAACTTACCCTTCGAAGCAACAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGTATc        >  2:1048571/1‑100 (MQ=255)
                                                                                            ggCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATc        >  2:3504238/1‑100 (MQ=255)
                                                                                            ggCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATc        <  1:3110019/100‑1 (MQ=255)
                                                                                            ggCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATc        >  2:824813/1‑100 (MQ=255)
                                                                                            ggCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATc        <  1:2713357/100‑1 (MQ=255)
                                                                                             gCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATc        >  2:3317465/1‑99 (MQ=255)
                                                                                               aTGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTc     >  1:2148050/1‑100 (MQ=255)
                                                                                                tgtgATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTcc    >  1:875516/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
TTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAG  >  NC_000913/1186993‑1187190

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