Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,187,090 G→C 15.0% intergenic (+20/+29) 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,0900GC15.0% 479.9 / 71.1 174intergenic (+20/+29)pepT/roxApeptidase T/50S ribosomal protein L16 arginine hydroxylase; 2‑oxoglutarate oxygenase
Reads supporting (aligned to +/- strand):  ref base G (78/69);  new base C (11/15);  total (89/85)
Fisher's exact test for biased strand distribution p-value = 3.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.79e-01

Only 100 of 186 total aligned reads displayed.

AGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCC  >  NC_000913/1186991‑1187188
                                                                                                   |                                                                                                  
aGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCAtg                                                                                                    <  2:3763664/100‑1 (MQ=255)
 gTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCAtgt                                                                                                   <  1:1335526/100‑1 (MQ=255)
 gTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCAtgt                                                                                                   <  2:213363/100‑1 (MQ=255)
 gTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCAtgt                                                                                                   >  1:3307182/1‑100 (MQ=255)
 gTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCAtgt                                                                                                   >  2:3406439/1‑100 (MQ=255)
  tttGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCAtgtg                                                                                                  <  2:1949438/100‑1 (MQ=255)
  tttGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCAtgtg                                                                                                  <  1:545554/100‑1 (MQ=255)
   ttGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGa                                                                                                 >  1:1380509/1‑100 (MQ=255)
      ttaCTCTGAAAAGAATAAAAAAAGCGGTGCAGGTGATCGTCCATATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGAt                                                                                                <  1:2435282/96‑1 (MQ=255)
      tgACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGc                                                                                              >  2:3701993/1‑100 (MQ=255)
       gACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGcc                                                                                             >  2:2364061/1‑100 (MQ=255)
       gACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGcc                                                                                             >  2:2922376/1‑100 (MQ=255)
        aCTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCg                                                                                            <  2:2222611/100‑1 (MQ=255)
         ctctGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGc                                                                                           >  2:1637918/1‑100 (MQ=255)
         ctctGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGc                                                                                           <  2:2895850/100‑1 (MQ=255)
           ctGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCAt                                                                                         >  2:3147444/1‑100 (MQ=255)
           ctGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCAt                                                                                         <  2:1848157/100‑1 (MQ=255)
            tGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCATc                                                                                        <  2:1825287/100‑1 (MQ=255)
            tGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCATc                                                                                        <  1:2224234/100‑1 (MQ=255)
            tGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCATc                                                                                        <  2:2109125/100‑1 (MQ=255)
            tGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCATc                                                                                        <  2:2292595/100‑1 (MQ=255)
            tGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCATc                                                                                        <  2:406028/100‑1 (MQ=255)
            tGGAAGGTATGGAAAAAGCGGTGCAGCCGATCGTCGTTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCATc                                                                                        <  1:1948891/100‑1 (MQ=255)
             ggaaggaaAGAAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCATcc                                                                                       <  1:2850546/100‑1 (MQ=255)
             ggAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGAAAGAGAAAAGGGATGCGGCATGTGATGCCGCATcc                                                                                       >  2:201716/1‑100 (MQ=255)
                 ggTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCt                                                                                   >  2:3836104/1‑100 (MQ=255)
                   tATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTa                                                                                 >  1:3710492/1‑100 (MQ=255)
                   tATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTa                                                                                 >  1:1151961/1‑100 (MQ=255)
                   tATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTa                                                                                 <  1:3809082/100‑1 (MQ=255)
                   tATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTa                                                                                 <  2:2699061/100‑1 (MQ=255)
                     tGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTaaa                                                                               <  2:3655617/100‑1 (MQ=255)
                     tGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCATcccttttcgc                                                                               <  1:696082/100‑9 (MQ=37)
                      ggAAAAATGGGGGCAGGTGATCGTCCGTATTGCCGATTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTaaa                                                                               <  1:3832206/99‑1 (MQ=255)
                      ggAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAAt                                                                              <  1:4331423/100‑1 (MQ=255)
                       gAAAAAGCGGTGCAGGTTATCGTCCATATTGCCGATTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATc                                                                             <  1:461455/100‑1 (MQ=255)
                       gAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATc                                                                             <  2:594458/100‑1 (MQ=255)
                       gAAAAAACGGTACAGGTGATCGTCCGTATTGCCAAGTTAACGCCGCAACGGGAGTATGCGAGAAGGGATGCGGCATGTGATGCCGCATCCGGCTTACATc                                                                             >  1:4111535/1‑100 (MQ=255)
                        aaaaaGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATcc                                                                            <  2:4177012/100‑1 (MQ=255)
                         aaaaGCGGTGCAGCTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGTATGCGGCATGTGATGCCGCATCCGGCTTAAATCCa                                                                           >  1:813199/1‑100 (MQ=255)
                          aaaGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCaa                                                                          >  1:3666835/1‑100 (MQ=255)
                          aaaGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCaa                                                                          <  1:4256397/100‑1 (MQ=255)
                            aGCGGTGCAGGTGATCGGCCGTATTTCCGAGTTAACGGCGCAACGGAAGTAAGCGAGAAGGGATGCGTCCTGGGATGCCGCATCCGGCTTAAATCCAAAc                                                                        >  1:2813013/1‑100 (MQ=255)
                                gTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTa                                                                     >  2:4001465/1‑99 (MQ=255)
                                gTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTAc                                                                    <  2:3789166/100‑1 (MQ=255)
                                  gCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTAccc                                                                  <  2:3332160/100‑1 (MQ=255)
                                  gCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTAccc                                                                  <  1:2482221/100‑1 (MQ=255)
                                     ggTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCtt                                                                >  2:678859/1‑99 (MQ=255)
                                     ggTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTc                                                               >  2:1770272/1‑100 (MQ=255)
                                       tGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCga                                                             <  1:1797470/100‑1 (MQ=255)
                                       tGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCga                                                             <  1:2254796/100‑1 (MQ=255)
                                        gATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCgaa                                                            >  2:4256241/1‑100 (MQ=255)
                                        gATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCgaa                                                            >  2:3213389/1‑100 (MQ=255)
                                         aTCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCgaag                                                           <  1:919858/100‑1 (MQ=255)
                                         aTCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCgaag                                                           >  1:3230873/1‑100 (MQ=255)
                                           cGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCgaagaa                                                         <  2:1764823/100‑1 (MQ=255)
                                            gTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAAc                                                        <  1:3197390/100‑1 (MQ=255)
                                             tCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAAcc                                                       <  1:3077414/100‑1 (MQ=255)
                                             tCCGTAGTGCCGAGTTAACGGCGCAACGGAATTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCgaatgacc                                                       <  1:2624909/100‑6 (MQ=255)
                                              ccGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCa                                                      <  1:1001060/100‑1 (MQ=255)
                                                gTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAAt                                                    >  1:1809729/1‑100 (MQ=255)
                                                  aTTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATAc                                                  <  1:4357167/100‑1 (MQ=255)
                                                   ttGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATAcc                                                 <  2:2616492/100‑1 (MQ=255)
                                                   tggCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATAcc                                                 <  2:2555208/98‑1 (MQ=255)
                                                     gCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATAcc                                                 <  1:3381820/98‑1 (MQ=255)
                                                     gCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCg                                               >  2:310247/1‑100 (MQ=255)
                                                      ccGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGc                                              >  2:2411920/1‑100 (MQ=255)
                                                       cGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCt                                             >  2:652779/1‑100 (MQ=255)
                                                        gAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTa                                            >  1:3314275/1‑100 (MQ=255)
                                                         aGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTAt                                           <  2:4345301/100‑1 (MQ=255)
                                                          gTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTAtt                                          >  2:1533977/1‑100 (MQ=255)
                                                           ttAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTg                                         <  1:2553179/100‑1 (MQ=255)
                                                              aCGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGAcc                                      <  1:1686890/100‑1 (MQ=255)
                                                                ggCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAg                                    <  2:1126986/100‑1 (MQ=255)
                                                                  cgcAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCa                                     >  1:2669561/1‑97 (MQ=255)
                                                                   gcAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAgcgc                                 <  1:127781/100‑1 (MQ=255)
                                                                        ggAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGa                            >  1:3610991/1‑100 (MQ=255)
                                                                           aGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCa                         >  1:2831505/1‑100 (MQ=255)
                                                                           aGTAAGAGAAAAGGGATGCGGCATGTGATTCCGCATCCGGCTTAAATCCAAACTTACCCTTCGCAGAACCAATACCCGCTATTGACCAGCGCCGCGAGc                          <  2:2497996/99‑1 (MQ=255)
                                                                               aGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATcgc                     >  1:3909223/1‑100 (MQ=255)
                                                                                gCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTAAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATcgcg                    >  2:3582667/1‑100 (MQ=255)
                                                                                  gAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTGGAAGAACCAATACCCGCTATTGACCAGCGACGCGAGCATCGCGa                   >  1:3975301/1‑99 (MQ=255)
                                                                                   aaaaGGGATGCGGCATGTGATGCCGCATCCGGCTTACATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAgg                 <  2:4206007/100‑1 (MQ=255)
                                                                                     gccggATGCGGCATCACATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGaa               >  2:3144412/4‑100 (MQ=38)
                                                                                     gccggATGCGGCATCACATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGaa               >  1:2865168/4‑100 (MQ=38)
                                                                                       gggATGCTGCATGTGATGCCGCATCCGTCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAAt              <  1:611894/99‑1 (MQ=255)
                                                                                       gggATGCTGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATATCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATg             <  1:2364051/100‑1 (MQ=255)
                                                                                       gggATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATg             >  1:65014/1‑100 (MQ=255)
                                                                                       cggATGCGGCATCACATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATg             >  1:4267117/2‑100 (MQ=39)
                                                                                        ggATGCGGCATCACATGCCGCATCCGGCTTAAATCCACACTTACCCTTCGAAGAACCAATACCCGCTATTGACCCGCGCCGCGCGCATCGCGacccatga            >  1:3440569/1‑93 (MQ=255)
                                                                                         gATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGAc           >  2:800166/1‑100 (MQ=255)
                                                                                         gATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGAc           >  1:2149238/1‑100 (MQ=255)
                                                                                         gATGCGGCATCACATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGAc           >  2:4406190/1‑100 (MQ=255)
                                                                                         gATGCGGCATCACATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGAc           >  2:1295398/1‑100 (MQ=255)
                                                                                          aTGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACg          <  2:4275029/100‑1 (MQ=255)
                                                                                            atctcATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGa        <  2:3493240/96‑1 (MQ=255)
                                                                                              ggCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATc      >  2:2956316/1‑100 (MQ=255)
                                                                                                cATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGaagat    >  1:3493217/1‑96 (MQ=255)
                                                                                                cATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGAAGGATCtt    <  1:1704418/100‑1 (MQ=255)
                                                                                                  tgtgATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTcc  <  1:2094673/100‑1 (MQ=255)
                                                                                                  tgtgATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTcc  >  2:3768950/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
AGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCC  >  NC_000913/1186991‑1187188

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: