Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,187,092 G→C 17.5% intergenic (+22/+27) 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,0920GC17.5% 438.3 / 84.2 171intergenic (+22/+27)pepT/roxApeptidase T/50S ribosomal protein L16 arginine hydroxylase; 2‑oxoglutarate oxygenase
Reads supporting (aligned to +/- strand):  ref base G (77/64);  new base C (12/18);  total (89/82)
Fisher's exact test for biased strand distribution p-value = 1.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.93e-01

Only 100 of 180 total aligned reads displayed.

TTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGC  >  NC_000913/1186993‑1187191
                                                                                                   |                                                                                                   
tttGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCAtgtg                                                                                                     <  2:1184728/100‑1 (MQ=255)
 ttGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGa                                                                                                    >  1:113602/1‑100 (MQ=255)
    tgACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGc                                                                                                 >  1:2093532/1‑100 (MQ=255)
    tgACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGc                                                                                                 >  1:2073419/1‑100 (MQ=255)
    tgACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGc                                                                                                 <  1:2923292/100‑1 (MQ=255)
    tgACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGc                                                                                                 >  2:3019372/1‑100 (MQ=255)
    tgAATCTGGAAGGTATGGAAAAAGCGGTGCAGCTGATCGTCCGTATTGCCAACTTAACCGCGCAACGGAAGTAAGCGCAAAGGGATGCGGCATGTGATGc                                                                                                 <  1:907574/100‑1 (MQ=255)
     gACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGcc                                                                                                >  1:1900633/1‑100 (MQ=255)
       ctctGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGc                                                                                              >  1:631267/1‑100 (MQ=255)
       ctctGGAAGGGATGGAAAAAGCGGTGCAGGTGATCGGCCGGATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGcc                                                                                                <  1:1113751/98‑1 (MQ=255)
         ctGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCAt                                                                                            >  1:896388/1‑100 (MQ=255)
         ctGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCAt                                                                                            >  1:1561452/1‑100 (MQ=255)
          tGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCATc                                                                                           <  1:3029733/100‑1 (MQ=255)
          tGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCATc                                                                                           <  1:2716535/100‑1 (MQ=255)
          tGGAAGGCATGAAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGGAAGCGAAAAGGGATGCGGCATGTGATGCCGCATc                                                                                           >  2:1909795/1‑100 (MQ=255)
           ggAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATcc                                                                                          >  1:721597/1‑100 (MQ=255)
           ggAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATcc                                                                                          >  1:2259679/1‑100 (MQ=255)
           ggAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATc                                                                                           <  2:220162/99‑1 (MQ=255)
           ggAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCATcc                                                                                          <  1:1858408/100‑1 (MQ=255)
           ggAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCATcc                                                                                          <  1:491919/100‑1 (MQ=255)
            gAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCATccc                                                                                         <  1:2698258/100‑2 (MQ=39)
             aaGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCATcc                                                                                          <  2:3553861/98‑1 (MQ=39)
                 tATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTa                                                                                    <  2:1561444/100‑1 (MQ=255)
                  aTGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCATcccttttcg                                                                                   <  2:13951/100‑8 (MQ=37)
                    ggAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAAt                                                                                 >  2:12310/1‑100 (MQ=255)
                    ggAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCATcccttttcgct                                                                                 <  1:3394696/100‑10 (MQ=37)
                    ggAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGAAAGAAAGAGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAAt                                                                                 <  2:569548/100‑1 (MQ=255)
                     gAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATc                                                                                >  1:219528/1‑100 (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)
                          aGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAAc                                                                           >  2:1528073/1‑100 (MQ=255)
                          aGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAAc                                                                           <  1:673564/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:1864549/100‑1 (MQ=255)
                               tGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTAcc                                                                      >  2:3587416/1‑100 (MQ=255)
                                gCAGGTGATCGTCAGTATTGCGGAGTTAACGGCGCAAAAAAAAAAAAAAAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACt                                                                          <  2:965713/95‑1 (MQ=255)
                                  aGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCtt                                                                   >  2:2483383/1‑100 (MQ=255)
                                      gATCGTCCGTATTGCCGAGTTAACGGCCCAACGGCAGTAAGCGAACACGGCTGCGGCATGTGATCCCGCATCCGACCTCACTCCACACTTAccc                                                                     >  2:3643236/1‑94 (MQ=255)
                                       aTCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCgaag                                                              <  1:2483399/100‑1 (MQ=255)
                                        tCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCgaaga                                                             <  1:542316/100‑1 (MQ=255)
                                           tCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGCAGAAcc                                                          >  2:420224/1‑100 (MQ=255)
                                           tCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAAcc                                                          <  2:2423302/100‑1 (MQ=255)
                                           tCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAAcc                                                          <  2:2991564/100‑1 (MQ=255)
                                               tATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAATGGATTCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCTTTCGATGAACCAAt                                                       >  2:1700/1‑99 (MQ=255)
                                                  tGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCCAACTTACCCTTCGAAGAACCAATAccc                                                   >  1:2830172/1‑100 (MQ=255)
                                                  tGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATAcc                                                    <  2:2363589/99‑1 (MQ=255)
                                                   gCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCg                                                  >  1:1699365/1‑100 (MQ=255)
                                                    ccGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGc                                                 <  1:2896502/100‑1 (MQ=255)
                                                    ccGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGc                                                 <  1:1441355/100‑1 (MQ=255)
                                                    ccGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGc                                                 <  1:281977/100‑1 (MQ=255)
                                                       aGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTAt                                              >  1:1329819/1‑100 (MQ=255)
                                                         ttAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTg                                            <  2:3369573/100‑1 (MQ=255)
                                                          tAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGa                                           >  2:2332832/1‑100 (MQ=255)
                                                           aaCGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGAc                                          >  1:2017762/1‑100 (MQ=255)
                                                           aaCGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGAc                                          >  1:3264451/1‑100 (MQ=255)
                                                           aaCGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGAc                                          >  2:1526704/1‑100 (MQ=255)
                                                            aCGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGAcc                                         <  1:1081413/100‑1 (MQ=255)
                                                             cGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGAcc                                         >  2:208976/1‑99 (MQ=255)
                                                             cGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCa                                        <  1:2117498/100‑1 (MQ=255)
                                                             cGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACGTACCCTTCGAAGAACCAATACCCGCTATTGACCa                                        >  2:3162348/1‑100 (MQ=255)
                                                                cgcAACGGAATGAAGCGAAAAGGGATGCTGCATGTCACGCCGCATCCGGCTTAAATCCAAACTTACTCTTCGAACAACCAATACCCGCTAGTGAcc                                         <  2:3355907/96‑1 (MQ=255)
                                                                 gcAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAgcgc                                    >  2:507665/1‑100 (MQ=255)
                                                                  cAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGcgcc                                   >  1:761376/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                              <  2:3328299/100‑1 (MQ=255)
                                                                       gAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAg                              >  1:2503199/1‑100 (MQ=255)
                                                                       gAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAg                              <  1:2815244/100‑1 (MQ=255)
                                                                       gAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCCAAGAACCAATACCCGCTATTGACCAGCGCCGCGAg                              <  1:343014/100‑1 (MQ=255)
                                                                       gAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAAGCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAg                              >  2:1753128/1‑100 (MQ=255)
                                                                        aaataaaaaaaaaGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAg                              <  2:2638247/90‑1 (MQ=255)
                                                                           tAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATc                          >  1:173256/1‑100 (MQ=255)
                                                                           tAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATc                          >  2:3234263/1‑100 (MQ=255)
                                                                                gAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAg                     >  1:3317111/1‑100 (MQ=255)
                                                                                gAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAg                     >  1:213382/1‑100 (MQ=255)
                                                                                gAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAg                     >  1:1806278/1‑100 (MQ=255)
                                                                                 aaaaGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAgg                    <  1:2428391/100‑1 (MQ=255)
                                                                                   gccggATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGaa                  >  1:2243189/4‑100 (MQ=255)
                                                                                   gccggATGCGGCATCACATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGaa                  >  1:70909/4‑100 (MQ=38)
                                                                                    ccggATGCGGCATCACATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAAt                 >  1:1853568/3‑100 (MQ=38)
                                                                                    ccggACGCGGCATCACATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACAAGCGCCGCGAGCATCGCGAGGAAt                 >  2:2695512/3‑100 (MQ=255)
                                                                                     cggATGCGGCATCACATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATg                >  2:555396/2‑100 (MQ=39)
                                                                                      ggATGCGGCATCACATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGa               >  1:2069844/1‑100 (MQ=255)
                                                                                      ggATGCGGCATCACATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGa               >  1:3599357/1‑100 (MQ=255)
                                                                                      ggATGCGGCATCACATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGa               >  1:286499/1‑100 (MQ=255)
                                                                                        aTGCGGCATGTGATGCCGCACCCGGCTTAAATCCAAGCTTACCCTTCGAAGAACCAATACCCGCGATTGACCAGCGCCGCGAGCATCGCGAGGAATGACg             <  1:2553432/100‑1 (MQ=255)
                                                                                         tGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACg             <  1:43291/99‑1 (MQ=255)
                                                                                         tGCGGCATGTGATACCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACgg            >  2:2253696/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)
                                                                                             gCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATc         >  2:3317465/1‑99 (MQ=255)
                                                                                                tgtgATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTcc     >  2:1010946/1‑100 (MQ=255)
                                                                                                tgtgATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTcc     <  2:1329811/100‑1 (MQ=255)
                                                                                                   gATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAgc  >  1:1526020/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGC  >  NC_000913/1186993‑1187191

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