Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,472,434 G→A 17.6% R34* (CGA→TGA)  pyrL ← pyrBI operon leader peptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,472,4340GA17.6% 448.3 / 78.0 194R34* (CGA→TGA) pyrLpyrBI operon leader peptide
Reads supporting (aligned to +/- strand):  ref base G (88/71);  new base A (25/9);  total (114/80)
Fisher's exact test for biased strand distribution p-value = 5.68e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01

Only 100 of 196 total aligned reads displayed.

GGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCA  >  NC_000913/4472335‑4472532
                                                                                                    |                                                                                                  
 gTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGa                                                                                                   >  2:1459621/1‑100 (MQ=255)
   cATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGAtt                                                                                                 >  2:1262227/1‑100 (MQ=255)
   cATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGAtt                                                                                                 <  2:3969045/100‑1 (MQ=255)
       gcgACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAgg                                                                                             >  1:3141287/1‑100 (MQ=255)
        cgACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAggg                                                                                            <  2:4027053/100‑1 (MQ=255)
        cgACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAggg                                                                                            >  1:466473/1‑100 (MQ=255)
        cgACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAggg                                                                                            >  1:101532/1‑100 (MQ=255)
        cgACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAggg                                                                                            <  1:1744089/100‑1 (MQ=255)
        cgACTAAGGTCGTTTATGGAAATGACATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAgg                                                                                             <  2:3080032/99‑1 (MQ=255)
        cgACTAAGGTAGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCGGGGCAAAAAAAAGCCCCTCGATTGAgg                                                                                             <  1:1882054/99‑1 (MQ=255)
          aCTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGc                                                                                          >  2:1601404/1‑100 (MQ=255)
           cTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCAATCGAGGGGCt                                                                                         <  2:1077658/100‑1 (MQ=255)
             aaGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTgg                                                                                       >  1:3445841/1‑100 (MQ=255)
             aaGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCAATCGAGGGGCttt                                                                                       <  1:1560276/100‑3 (MQ=255)
             aaGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCAATCGAGGGGCttt                                                                                       <  1:1337260/100‑3 (MQ=255)
              aGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTggg                                                                                      <  2:214640/100‑1 (MQ=255)
              aGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTggg                                                                                      <  2:493336/100‑1 (MQ=255)
                  cGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCAATCGAGGGGCtttttttt                                                                                  <  2:310203/100‑8 (MQ=39)
                   gTTTATGGAAGTGCTATGTTTCTGATATAGCTGATTATCCTTCTTTTATCTCCTGACGCCTTAGCAAAAAAAAGCCCCTCAATCGAGGGGCttttttttg                                                                                 <  2:822896/100‑9 (MQ=25)
                    tttATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAACAAGCCCCTCGATTGAGCGGCTGGGAAt                                                                                   >  1:4211900/1‑97 (MQ=255)
                      tacggAAATCATATGTTTCTGATATAGCGAATTAGCCATCTTTTATCTAAAAAAAAAAAAAAAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTg                                                                              <  1:3683589/97‑1 (MQ=255)
                      tATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTg                                                                              >  2:3435535/1‑100 (MQ=255)
                      tATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTg                                                                              >  2:2923981/1‑100 (MQ=255)
                       aaaggaaTGATATGTTTATGATATAGCGGAATAGCCAGCTTGTATCTCCTGACGCAGGAGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTg                                                                              <  1:2300073/94‑1 (MQ=255)
                        tGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAGAGCCCCTCGATTGAGGGGCTGGGAATGGGTGAt                                                                            >  2:3913670/1‑100 (MQ=255)
                          gAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCa                                                                          <  1:3508991/100‑1 (MQ=255)
                          gAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCa                                                                          <  2:705329/100‑1 (MQ=255)
                             aTGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACg                                                                       <  2:3467118/100‑1 (MQ=255)
                                  atGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGgaan                                                                  <  1:3000097/100‑2 (MQ=255)
                                    gTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGgaag                                                                  <  1:1192168/98‑1 (MQ=255)
                                     tttCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGaaa                                                               >  2:3809873/1‑100 (MQ=255)
                                      ttCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGaaaa                                                              >  2:1139952/1‑100 (MQ=255)
                                      ttCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGaaaa                                                              <  2:3169678/100‑1 (MQ=255)
                                      ttCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGaaaa                                                              >  1:953743/1‑100 (MQ=255)
                                         tGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACg                                                           <  2:2016321/100‑1 (MQ=255)
                                         tGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACg                                                           >  1:1475115/1‑100 (MQ=255)
                                          gATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACgg                                                          >  2:3739165/1‑100 (MQ=255)
                                           atatAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGc                                                         >  2:4087053/1‑100 (MQ=255)
                                            tataGCGGATTAGCCATTTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCa                                                        <  1:296284/100‑1 (MQ=255)
                                            tataGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCa                                                        >  1:640601/1‑100 (MQ=255)
                                               aGCGGATGAGCCATTTTTTATTTCCTGAGGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGc                                                     <  1:2551679/100‑1 (MQ=255)
                                                 cGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCa                                                   >  2:1631061/1‑100 (MQ=255)
                                                  ggATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAg                                                  <  2:315784/100‑1 (MQ=255)
                                                    aTTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCg                                                <  2:3780352/100‑1 (MQ=255)
                                                        gCCATCTTTTATCTCCTGACGCCTGGGCAAAAAACAGCCCCTCGCTTGAGGGGCTGGGAATGGGTGACCAACGGGAAGAAAAACGGCAGGCCAGCGTCtt                                            >  1:994909/1‑100 (MQ=255)
                                                        gCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCtt                                            >  1:1071868/1‑100 (MQ=255)
                                                        gCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCtt                                            >  1:2459590/1‑100 (MQ=255)
                                                             cTTTTATCTCCTGACGCCTGAGCAAAAAAAAACCCATCGATTGAGGGGCTGGCAATGGGTAAAAAAAAAAAAAAAAAACGGCAGGCCAGCGTCTTTTTTc                                       <  1:4053632/100‑1 (MQ=255)
                                                              ttttATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCa                                      <  2:4395145/100‑1 (MQ=255)
                                                              ttttATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCa                                      <  2:2460145/100‑1 (MQ=255)
                                                               tttATCTCCTGACGCCTGGGCAAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCa                                      <  2:838725/100‑1 (MQ=255)
                                                                 tATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGAc                                   >  1:1199932/1‑100 (MQ=255)
                                                                     tcCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCgg                               >  2:637771/1‑100 (MQ=255)
                                                                      ccTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGt                              >  2:1809325/1‑100 (MQ=255)
                                                                      ccTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGt                              >  2:4219217/1‑100 (MQ=255)
                                                                      ccTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGt                              >  2:3641952/1‑100 (MQ=255)
                                                                      ccTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGt                              <  2:48556/100‑1 (MQ=255)
                                                                       cTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGAcgctggc                             <  2:349831/100‑5 (MQ=255)
                                                                       cTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTa                             <  1:1459615/100‑1 (MQ=255)
                                                                        tGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTaa                            >  2:3367018/1‑100 (MQ=255)
                                                                         gACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAg                           <  2:159022/100‑1 (MQ=255)
                                                                          aCGCCTAGGCAAAAAAAAGCCCCTCGATTGAGGCGCTGGTACTGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTac                            <  1:2734248/98‑2 (MQ=255)
                                                                           cGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGAc                         >  2:1309675/1‑100 (MQ=255)
                                                                            gCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACa                        >  2:2450274/1‑100 (MQ=255)
                                                                             ccTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACaa                       <  2:1751148/100‑1 (MQ=255)
                                                                              cTGGGCAAAAAAAAGCCCCTCGAGTGAGGGGCCAGGCAGACGTCATCCACGGGCAGAACAACGGCAGGTCAGCGGCATGTTTCAGACGCGGTAAGACa                        >  1:2060953/1‑98 (MQ=255)
                                                                                gggCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAAt                    >  2:936901/1‑100 (MQ=255)
                                                                                 ggCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATg                   >  2:2386773/1‑100 (MQ=255)
                                                                                 ggCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATg                   >  2:2603900/1‑100 (MQ=255)
                                                                                 ggCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATg                   >  1:872450/1‑100 (MQ=255)
                                                                                  gCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATg                   >  2:232546/1‑99 (MQ=255)
                                                                                   cccattcccaGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATg                   >  1:4276566/10‑99 (MQ=38)
                                                                                   cAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCTGACTTTGTAAGACAAAATGTc                 >  1:76401/1‑100 (MQ=255)
                                                                                   cAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTc                 >  2:3990063/1‑100 (MQ=255)
                                                                                     cattcccaGCCCCTCAATCGAGGGGCTTGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCg                >  2:375234/8‑100 (MQ=38)
                                                                                     cattcccaGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGCCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCg                >  1:564775/8‑100 (MQ=38)
                                                                                      attcccaGCCCCTCAATCGAGGGGCTGGGAATGGATGATCAACGAGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGt                  >  1:512759/7‑97 (MQ=38)
                                                                                      aaaaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGa               >  2:1573557/1‑100 (MQ=255)
                                                                                       ttcccaGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGaa              >  2:2406614/6‑100 (MQ=255)
                                                                                        cttccgatcttCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAac             <  1:720072/90‑1 (MQ=255)
                                                                                        aaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAac             >  2:683063/1‑100 (MQ=255)
                                                                                         cccaGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAaca            >  2:218399/4‑100 (MQ=255)
                                                                                         cccaGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAaca            >  1:2635748/4‑100 (MQ=255)
                                                                                         cccaGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAaca            >  1:3716250/4‑100 (MQ=255)
                                                                                          ccaGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCATACGTGGTAAGACAAAATGTCGAacac           >  1:2216465/3‑100 (MQ=255)
                                                                                          ccaGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAacac           >  2:332087/3‑100 (MQ=255)
                                                                                          ccaGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAacac           >  2:2707688/3‑100 (MQ=255)
                                                                                          ccaGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAacac           >  1:3988805/3‑100 (MQ=255)
                                                                                          ccaGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAaca            >  1:689224/3‑99 (MQ=255)
                                                                                          aaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAacac           >  2:961666/1‑100 (MQ=255)
                                                                                           caGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAacaca          >  2:415594/2‑100 (MQ=255)
                                                                                           caGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAacaca          >  2:151673/2‑100 (MQ=255)
                                                                                           aaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAacaca          >  1:1716005/1‑100 (MQ=255)
                                                                                            aGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAacacac         >  2:1373909/1‑100 (MQ=255)
                                                                                            aGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAacacac         >  2:2769818/1‑100 (MQ=255)
                                                                                              ccccTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTg       >  1:2616910/1‑100 (MQ=255)
                                                                                                ccTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTTTTACTTTTTAAGACAAAATGTCGAACACACTGaa     >  1:431352/1‑100 (MQ=255)
                                                                                                ccTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGaa     >  1:2791410/1‑100 (MQ=255)
                                                                                                ccTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGaa     >  2:1978808/1‑100 (MQ=255)
                                                                                                   cGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCa  <  2:1180220/100‑1 (MQ=255)
                                                                                                    |                                                                                                  
GGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCA  >  NC_000913/4472335‑4472532

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: