Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,472,434 G→A 7.2% 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,4340GA 7.2% 440.1 / 15.4 138R34* (CGA→TGA) pyrLpyrBI operon leader peptide
Reads supporting (aligned to +/- strand):  ref base G (70/58);  new base A (5/5);  total (75/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.65e-01

Only 100 of 139 total aligned reads displayed.

TTGCGGGTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGC  >  NC_000913/4472288‑4472582
                                                                                                                                                   |                                                                                                                                                    
ttGCGGGTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGAt                                                                                                                                                    <  3:254510/149‑1 (MQ=255)
    gggTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGa                                                                                                                                                 <  2:813737/148‑1 (MQ=255)
    gggTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGa                                                                                                                                                 <  2:785746/148‑1 (MQ=255)
    gggTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGa                                                                                                                                                 <  2:28821/148‑1 (MQ=255)
    gggTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAg                                                                                                                                                <  4:681761/149‑1 (MQ=255)
    gggTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAg                                                                                                                                                <  3:179399/149‑1 (MQ=255)
     ggTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAg                                                                                                                                                >  1:130528/1‑148 (MQ=255)
       tttGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAggg                                                                                                                                              >  1:487615/1‑148 (MQ=255)
        ttGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCAATCGAGGGGc                                                                                                                                            <  3:445760/149‑1 (MQ=255)
          gCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCAATCGAGGGGCt                                                                                                                                           <  1:745882/148‑1 (MQ=255)
           cTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTgg                                                                                                                                         <  3:933843/149‑1 (MQ=255)
            tttCAGTTTGGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AGAAAAAAGCCCCTCGATTGAGGGGCTgga                                                                                                                                        >  4:822471/1‑148 (MQ=255)
            tttCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTggg                                                                                                                                        >  3:787674/1‑149 (MQ=255)
            tttCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTGATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCAAAATTTAGGGCCTggg                                                                                                                                        >  3:608991/1‑149 (MQ=255)
             ttCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTggg                                                                                                                                        >  1:967175/1‑148 (MQ=255)
             ttCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTggg                                                                                                                                        >  1:932994/1‑148 (MQ=255)
             ttCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCAATCGAGGGGCtttt                                                                                                                                        <  1:657290/148‑4 (MQ=255)
                aGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATg                                                                                                                                    >  4:535922/1‑149 (MQ=255)
                 tttccgccgcTGTCGCAAGAACAAGGTTAAGGTCATCGCGACAAGGGTGGTTTATGGAATTGATAGGTTTCTGATATAGCGGATTAGCCATCTTTAATCTCCTGCCGCTTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATg                                                                                                                                    <  1:330885/144‑1 (MQ=255)
                 gTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATgg                                                                                                                                   >  3:578487/1‑149 (MQ=255)
                  tttCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGTGGCTGGGAATggg                                                                                                                                  >  4:744940/1‑149 (MQ=255)
                  tttCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCAGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATgg                                                                                                                                   >  1:523874/1‑148 (MQ=255)
                   ttCGCCTCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGTTCGTTTATGGAAATGATATGTTTCTGAAATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAACAGCCCCTCGATTGAGGGGCTGGGAATggg                                                                                                                                  <  2:802862/148‑1 (MQ=255)
                   ttCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATggg                                                                                                                                  <  1:258725/148‑1 (MQ=255)
                     cgccgcTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGa                                                                                                                               <  4:714208/149‑1 (MQ=255)
                     cgccgcTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATCAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGa                                                                                                                               <  3:695372/149‑1 (MQ=255)
                         gcTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCaa                                                                                                                           <  3:219601/149‑1 (MQ=255)
                           tGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACg                                                                                                                         <  3:539466/149‑1 (MQ=255)
                             tCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACg                                                                                                                         <  2:487578/148‑1 (MQ=255)
                               gCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGga                                                                                                                      <  2:597149/148‑1 (MQ=255)
                                 cAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGgaag                                                                                                                    >  1:15956/1‑148 (MQ=255)
                                  aGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGgaagaa                                                                                                                  <  3:710154/149‑1 (MQ=255)
                                     aCCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGaaaaa                                                                                                               >  3:476745/1‑149 (MQ=255)
                                        aGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACg                                                                                                             >  1:509920/1‑148 (MQ=255)
                                               ggTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCa                                                                                                          >  4:440446/1‑144 (MQ=255)
                                                   aTCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGt                                                                                                 >  3:332069/1‑149 (MQ=255)
                                                    tCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATgg                                                                                                                                   <  4:216564/114‑1 (MQ=255)
                                                     cgcgACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTc                                                                                                >  1:365781/1‑148 (MQ=255)
                                                       cgACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCttt                                                                                             <  4:544160/149‑1 (MQ=255)
                                                         aCTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCAATCGAGGGGCtttttt                                                                                                                                      >  4:812462/1‑101 (MQ=255)
                                                         aCTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCAATCGAGGGGCtttttt                                                                                                                                      <  3:812462/106‑6 (MQ=255)
                                                          cTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCttttt                                                                                           <  1:644418/148‑1 (MQ=255)
                                                          cTAAGGTCGTTTATGGAAATGATATGTTACTGATATAGCGCATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCttttt                                                                                           <  2:130528/148‑1 (MQ=255)
                                                           tAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCTATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTGGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTc                                                                                         >  4:587378/1‑149 (MQ=255)
                                                           tAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTc                                                                                         <  4:332069/149‑1 (MQ=255)
                                                            aaGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGTGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCa                                                                                        <  3:744940/149‑1 (MQ=255)
                                                            aaGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTc                                                                                         >  1:395159/1‑148 (MQ=255)
                                                                tCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGAcg                                                                                    >  4:641449/1‑149 (MQ=255)
                                                                  gTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGAcgcg                                                                                  >  3:532137/1‑149 (MQ=255)
                                                                   tttATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGAcgcg                                                                                  >  2:284797/1‑148 (MQ=255)
                                                                      aTGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTaa                                                                              <  4:827663/149‑1 (MQ=255)
                                                                      aTGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTaa                                                                              <  3:521512/149‑1 (MQ=255)
                                                                       tGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTAATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAg                                                                             >  3:624252/1‑149 (MQ=255)
                                                                        ggAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAg                                                                             >  1:969121/1‑148 (MQ=255)
                                                                        ggAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGa                                                                            <  4:73038/149‑1 (MQ=255)
                                                                             tGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACaaaa                                                                       <  4:458851/149‑1 (MQ=255)
                                                                                    tttCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGa                                                                 <  2:406437/148‑1 (MQ=255)
                                                                                      tCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAac                                                               <  2:866042/148‑1 (MQ=255)
                                                                                           tataGCGGATTAGCCATCTTTTATCTCCTGACGCCGTGGC‑AGAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACt                                                          >  1:867827/1‑148 (MQ=255)
                                                                                            ataGCGGATTAGCAATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTg                                                         >  1:712769/1‑148 (MQ=255)
                                                                                              aGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAAc                                                      <  4:679825/149‑1 (MQ=255)
                                                                                                cGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAAcc                                                     >  2:744709/1‑148 (MQ=255)
                                                                                                       gCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACAtcc                                             >  4:274546/1‑149 (MQ=255)
                                                                                                        ccATCTTTTATCTCCTGACGCCTGGTC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACAtcct                                            >  4:885025/1‑149 (MQ=255)
                                                                                                        ccATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACAtcc                                             >  2:165083/1‑148 (MQ=255)
                                                                                                        ccATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACAtcc                                             >  2:241683/1‑148 (MQ=255)
                                                                                                           tCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACAtcctc                                           <  1:700032/147‑1 (MQ=255)
                                                                                                           tCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACAtcctc                                           >  2:700032/1‑147 (MQ=255)
                                                                                                           tCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTccc                                         <  3:641449/149‑1 (MQ=255)
                                                                                                             ttttATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCgg                                       >  3:659861/1‑149 (MQ=255)
                                                                                                             ttttATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCgg                                       >  3:577798/1‑149 (MQ=255)
                                                                                                             ttttATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCgg                                       >  3:624262/1‑149 (MQ=255)
                                                                                                              tttATCTCCTGACGCCTGGGC‑‑AAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACAtcctcc                                          <  1:541007/144‑1 (MQ=255)
                                                                                                              tttATCTCCTGACGCCTGGGC‑‑AAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACAtcctcc                                          >  2:541007/1‑144 (MQ=255)
                                                                                                                   ctcCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGgaa                                                                                                                     >  4:343301/1‑65 (MQ=255)
                                                                                                                      cTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGt                               <  2:395159/148‑1 (MQ=255)
                                                                                                                       tGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTcc                             <  4:162002/149‑1 (MQ=255)
                                                                                                                       tGACGCCTGGGC‑AAAAAAAAGCCCCGCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTACCGGAAAATTGTc                              >  2:232431/1‑148 (MQ=255)
                                                                                                                        gACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGc                                                                                                       <  2:745510/74‑1 (MQ=255)
                                                                                                                        gACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAg                                                                             <  3:353064/100‑1 (MQ=255)
                                                                                                                          cGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGc                          <  4:659861/149‑1 (MQ=255)
                                                                                                                           gCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCa                         <  4:62085/149‑1 (MQ=255)
                                                                                                                            ccTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCa                         >  2:87524/1‑148 (MQ=255)
                                                                                                                              tGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAg                                                                             >  4:596358/1‑94 (MQ=255)
                                                                                                                                ggC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATtat                     <  1:873071/148‑1 (MQ=255)
                                                                                                                                    aaaaaaaaGCCCCTGGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATAc                   <  4:42605/147‑1 (MQ=255)
                                                                                                                                    aaaaaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATAc                   >  3:42605/1‑147 (MQ=255)
                                                                                                                                    aaaaaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTc                 >  3:922722/1‑149 (MQ=255)
                                                                                                                                     attcccaGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTc                 >  2:395256/7‑148 (MQ=255)
                                                                                                                                       aaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCAt               >  2:548727/1‑148 (MQ=255)
                                                                                                                                          caGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGt            >  1:550365/2‑148 (MQ=255)
                                                                                                                                             ccccTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGa        >  4:851618/1‑149 (MQ=255)
                                                                                                                                               ccTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCttt                                                                                             <  1:115901/62‑1 (MQ=255)
                                                                                                                                                cTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAgc     >  3:448721/1‑149 (MQ=255)
                                                                                                                                                cTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAgc     >  3:853873/1‑149 (MQ=255)
                                                                                                                                                 tCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAgc     >  2:129861/1‑148 (MQ=255)
                                                                                                                                                 tCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAgc     >  1:638515/1‑148 (MQ=255)
                                                                                                                                                 tCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAgc     >  1:272834/1‑148 (MQ=255)
                                                                                                                                                  cGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAgcg    >  2:49793/1‑148 (MQ=255)
                                                                                                                                                   gATTGAGGGGCTGGGAATGGGTGATCAACGAGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAgcggc  >  3:115045/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TTGCGGGTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGC  >  NC_000913/4472288‑4472582

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑