Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,620,216 C→T 8.9% intergenic (‑56/+30) uraA ← / ← upp uracil permease/uracil phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,620,2160CT 8.9% 330.3 / 15.9 112intergenic (‑56/+30)uraA/uppuracil permease/uracil phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (53/49);  new base T (5/5);  total (58/54)
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.44e-01

Only 100 of 112 total aligned reads displayed.

ACCAAACATGGCGAACAAATGTTGCAAACTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCGACGTCCGGGTGCGCTTT  >  NC_000913/2620071‑2620362
                                                                                                                                                  |                                                                                                                                                  
aCCAAACATGGCGAACAAATGTTGCAAACTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTaa                                                                                                                                                >  1:410193/1‑150 (MQ=255)
        tGGCGAACAAATGTTGCAAACTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAACGCCGACTCTTAAAGTCGGCt                                                                                                                                        <  1:1560018/150‑1 (MQ=255)
        tGGCGAACAAATGTTGCAAACTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCt                                                                                                                                        <  1:1220920/150‑1 (MQ=255)
          gCGAACAAATGTTGCAAACTAAGCGGGATTGGCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCGGTATTATGTGTTATAGGCGCTTTACTCAAAAAAACGCCGACTCTTAAAGTCGGCttt                                                                                                                                      <  2:263897/150‑1 (MQ=255)
           cGAACAAATGTTGCAAACTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTa                                                                                                                                     >  2:1363602/1‑150 (MQ=255)
           cGAACAAATGTTGCAAACTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTa                                                                                                                                     >  1:1334238/1‑150 (MQ=255)
             aaCAATTGTTGCAATCTAAGCGGTATTGTCTGTAAATGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTTTTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTTTAAGAGTCGGCtttttt                                                                                                                                   <  2:214600/150‑4 (MQ=255)
              aCAAATTGTGGAAACCTACCGGGATTGTCTGTAAAAGTGGAGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTTTAAGAGTCGGCttttttt                                                                                                                                  <  1:917830/150‑5 (MQ=255)
              aCAAATGTTGCAAACTAAGCGGGATTGTCTGTAACAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTGTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAatt                                                                                                                                  <  1:304356/150‑1 (MQ=255)
                aaaTGTTGCAAACTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTTTAAGAGTCGGCttttttttg                                                                                                                                <  1:1074867/150‑7 (MQ=255)
                aaaTCTTGCAAACTTAGCGGGCTTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGGGTCATAGTATTATCCTCTGTATTTGGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTTTAAGAGTCGGCttttttttg                                                                                                                                <  1:1308603/150‑7 (MQ=255)
                   tGTTGCAACCTAACAGGGATTTTCTGTAAAAGTGTCGGTCTTTCACTCACCCCTATAGCACGGCGCGTCTTAGTATTATCCTCTGTATTTTTGGTTATAGGCGCTTAACTCAAAAAAACGCCGACTCTTAAAGTCGGCTTTAATTAtttt                                                                                                                             <  2:481986/150‑1 (MQ=255)
                    gTTGCAAACTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTAttttt                                                                                                                            <  2:414352/150‑1 (MQ=255)
                      tGCAAACTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTAt                                                                                                                          >  2:374538/1‑150 (MQ=255)
                       gCAAACTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTAtt                                                                                                                         >  2:873944/1‑150 (MQ=255)
                       gCAAACTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTAATATCCTCTGTATTATGTGTTATAGGCGCGTTACTCAAAAAAAAGCCGCATCTTAAAGTCGGCTTTAattattattatt                                                                                                                         >  1:1182042/1‑150 (MQ=255)
                        cAAACTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCt                                                                                                                       >  1:701183/1‑150 (MQ=255)
                         aaaCTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAATTCGGCTTTAATTATTTTTATTCt                                                                                                                       >  2:469795/1‑150 (MQ=255)
                         aaaCTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCt                                                                                                                       >  1:1234828/1‑150 (MQ=255)
                         aaaCTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGAATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCt                                                                                                                       >  1:123367/1‑150 (MQ=255)
                          aaCTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCtt                                                                                                                      <  2:1199546/150‑1 (MQ=255)
                           aCTAAGCGGGATTGGCTGTAAAAGTGGCGGTCTTTCACACACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTAAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCttt                                                                                                                     <  2:438406/150‑1 (MQ=255)
                             tAAGCGGGATTGTCTGTAAATGTGGCGGTCTTTCCCTCACGCCGATAGCACGGCGCGTCATATTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTAAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTAt                                                                                                                   <  2:1117430/150‑1 (MQ=255)
                              aaGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTAtt                                                                                                                  >  1:1186692/1‑150 (MQ=255)
                              aaGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTAATTATTTTTATTCTTTAttt                                                                                                                 >  1:1061016/1‑150 (MQ=255)
                               aGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTAttt                                                                                                                 <  2:1334238/150‑1 (MQ=255)
                                 cGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCg                                                                                                               >  1:869084/1‑150 (MQ=255)
                                  gggATTGTCTGTAAAAGTGGCTGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGt                                                                                                              <  1:1602705/150‑1 (MQ=255)
                                  gggATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTTTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGt                                                                                                              <  2:1354008/150‑1 (MQ=255)
                                   ggATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATATTCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTa                                                                                                             >  1:1689795/1‑150 (MQ=255)
                                   ggATTGTCTCTAAAAGTGGCGGTCTTTCACTAACCCCGATAGCACGGCGCGTCATAGTATTATAATCTGTATTAAGTGTTATAGGCGCTTTACTCAAAAAAAAACCGACTCTTAAAGAAGGCTTTAATTATTTTTATtatttattaataa                                                                                                             >  2:740457/1‑145 (MQ=255)
                                       tGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATCGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAACGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCaa                                                                                                         <  2:295607/150‑1 (MQ=255)
                                         tCTGTAAAAGTGGCGTTCTTTCACTCACCCCTATAGCACGGCGCGTAATATTTTTATCCTATGTTTTTTGTGTTATTGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAg                                                                                                       <  2:613281/150‑1 (MQ=255)
                                          cTGTAAAAGTGGCGGTCTTTCACTCCCCCCGATAGCACGGCGCGTCAAAGTATTATCCTCTGTATTATGTGTTAAAGGCGCTTTACTCAAAAAAAACCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGa                                                                                                      <  2:616001/150‑1 (MQ=255)
                                          cTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGa                                                                                                      >  2:113875/1‑150 (MQ=255)
                                          cTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGa                                                                                                      >  1:1085535/1‑150 (MQ=255)
                                          cTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGa                                                                                                      >  1:499712/1‑150 (MQ=255)
                                             tAAAAGTTGCGGTATTTAACTCACAACGATTGAACGGATAGTAATAGTATTATCCTCTGTATTATGTGTTATATTCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGAttt                                                                                                   <  1:1241842/150‑1 (MQ=255)
                                              aaaaGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTTTTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGAtttt                                                                                                  >  2:340627/1‑150 (MQ=255)
                                              aaaaGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGAtttt                                                                                                  >  1:1260261/1‑150 (MQ=255)
                                                        gTCTTTCACTCACCCCGATAGAACGGCGCGTAATATTATTATCCTATGTATTTTGTGTTATAGGCGCTTTACTAAAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGc                                                                                         <  2:636390/150‑1 (MQ=255)
                                                            ttaactcacCCCGATAGCACGGCGAGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAACCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGc                                                                                    <  1:353931/147‑1 (MQ=255)
                                                            ttCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGc                                                                                    >  2:607704/1‑150 (MQ=255)
                                                             taactcacCCCGATAGCACTGCGCGTCATATTATTATCCTCTTTATTATGTGTTATAGGCGCTTTACTTAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGcc                                                                                   <  2:1512387/148‑1 (MQ=255)
                                                               actcacCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAg                                                                                >  2:363176/1‑150 (MQ=255)
                                                                ataacCCCTATAGCACGGCGCGTAATATAATTATCCTCTTTATTTTGTTTTATAGGCTCTTTACTAAAAAAAAACCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAg                                                                                <  1:1708672/147‑1 (MQ=255)
                                                                 tcacCCCGATTGCACGGCGCGGCATAGTATTATCCTATGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAgg                                                                               <  1:622324/150‑1 (MQ=255)
                                                                  cacCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGc                                                                              >  1:1484621/1‑150 (MQ=255)
                                                                   acCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGcc                                                                             >  2:575235/1‑150 (MQ=255)
                                                                    ccccGATAGAACGGCGCGTCATAGTATTATCCTCTGTATTATGAGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCAACGGCATCTCCGAGGccc                                                                            >  1:889014/1‑150 (MQ=255)
                                                                     cacGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTTTAGGCGCTTTACTCAAAAATAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCg                                                                           <  2:626296/148‑1 (MQ=255)
                                                                          tAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGaata                                                                      <  2:494851/150‑1 (MQ=255)
                                                                            gCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGaataat                                                                    <  2:701183/150‑1 (MQ=255)
                                                                             cACGGCGCTTCATATTATTATCCAGTGTATTATGTGTAATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATg                                                                   <  1:746808/150‑1 (MQ=255)
                                                                             cACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGaataatt                                                                   >  2:925225/1‑149 (MQ=255)
                                                                             cACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATATCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATg                                                                   <  2:1673031/150‑1 (MQ=255)
                                                                                      tCATAGTATATACCTCTGTATTTTGTGTTATAGGCGCTTAACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGc                                                          <  2:715645/150‑1 (MQ=255)
                                                                                         tAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCg                                                       <  2:749828/150‑1 (MQ=255)
                                                                                           gTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTTGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGtt                                                     >  1:62981/1‑150 (MQ=255)
                                                                                           gTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGtt                                                     >  2:1177145/1‑150 (MQ=255)
                                                                                           gTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGtt                                                     >  2:839939/1‑150 (MQ=255)
                                                                                             attatCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCa                                                   >  2:220008/1‑150 (MQ=255)
                                                                                             attatCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCa                                                   >  2:1010128/1‑150 (MQ=255)
                                                                                               tatCCTCTGTATTATGTGTTATAGGAGCTTTACTCAAAAAAAAGCCGACTCTTAAATTCGTCTTTAATTATTTTTATTCTTTATTTCGTTACAAAGATATTCTCACCTGCATCGCCGAGTCCAGTAATAATGTATACGTGCTCGTTCAGt                                                 >  2:1083328/1‑150 (MQ=255)
                                                                                                  ccTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCACGGAATAATGTATCCGTGCTCGTTCAGTcac                                              >  1:1389772/1‑148 (MQ=255)
                                                                                                   ctctGTATTATGTGTTATAGTCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCt                                             <  2:646280/150‑1 (MQ=255)
                                                                                                   ctctGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCt                                             >  2:982000/1‑150 (MQ=255)
                                                                                                      tGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGAt                                          >  1:83721/1‑150 (MQ=255)
                                                                                                       gTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATc                                         <  2:410193/150‑1 (MQ=255)
                                                                                                         attatGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCaa                                       >  1:1186530/1‑150 (MQ=255)
                                                                                                              gtgtTATAGGCGCTTTACTCAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATg                                 <  1:363176/150‑1 (MQ=255)
                                                                                                               tgtTATAGGCGCTTTACTCAAAAAAACCCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATg                                 <  2:1537803/150‑1 (MQ=255)
                                                                                                               tgtTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATg                                 >  2:1352972/1‑150 (MQ=255)
                                                                                                               tgtTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGc                                <  2:1061016/150‑1 (MQ=255)
                                                                                                               tgtAATAGGCGCTTTACTCAAAAAAAACCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATg                                 <  2:1067494/150‑1 (MQ=255)
                                                                                                                 ttATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCg                               <  1:1060536/150‑1 (MQ=255)
                                                                                                                   ataGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGt                             <  1:957353/150‑1 (MQ=255)
                                                                                                                    taGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGta                            >  1:201872/1‑150 (MQ=255)
                                                                                                                     aGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGtat                           <  2:1605187/150‑1 (MQ=255)
                                                                                                                      ggCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGtata                          >  2:382885/1‑150 (MQ=255)
                                                                                                                      ggCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGtata                          <  1:607704/150‑1 (MQ=255)
                                                                                                                      ggCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGACCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGtata                          >  1:806725/1‑150 (MQ=255)
                                                                                                                      ggCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGATTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGTCATCGCCGAGGCGCGGAAAAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGtata                          >  1:1539400/1‑150 (MQ=255)
                                                                                                                         gcTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAg                       >  2:1008213/1‑150 (MQ=255)
                                                                                                                          cTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTTGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGt                      <  2:62981/150‑1 (MQ=255)
                                                                                                                            ttACTCAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCg                   >  1:988988/1‑150 (MQ=255)
                                                                                                                             ttcctaaaaaaaacCCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTc                    <  1:375810/144‑1 (MQ=255)
                                                                                                                              aCTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCGa                  >  1:1185007/1‑150 (MQ=255)
                                                                                                                               cTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCGAc                 >  1:763824/1‑150 (MQ=255)
                                                                                                                                  cAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGCGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCGCCgg               >  2:463537/1‑149 (MQ=255)
                                                                                                                                     aaaaaaGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCGACGTCCg            <  1:87496/150‑1 (MQ=255)
                                                                                                                                      ttaaaGCCGACTTTAAGAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGGCCCTGAGCAATCGATGCGGTAGACAGGTCGCCGTCCgg           >  2:1119298/3‑150 (MQ=255)
                                                                                                                                       taaaGCCGACTTTAAGAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATTTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCGACGTCCggg          >  1:437603/2‑150 (MQ=255)
                                                                                                                                        aaaGCCGACTTTAAGAGTCGGCTTTAATTATTTTTATTCGTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATAAGTGCTCGTTCAGTCCCTTATAAATCGATGCGGTATACAGTTCGACGTCTGGGt         >  1:170891/1‑150 (MQ=255)
                                                                                                                                         aaGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCGACGTCCGGGTg        >  2:225857/1‑150 (MQ=255)
                                                                                                                                          aGCCGACTTTAAGAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCGACGTCCGGGTgc       >  1:1096412/1‑150 (MQ=255)
                                                                                                                                           gCCGACTTTAAGAGTCGGCTTTAATTATTTTTATTCTTTAATTCGTACCAAAGATTTTGTCACCGGAATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCGACGTCCGGTTgcg      >  2:1018430/1‑150 (MQ=255)
                                                                                                                                              gACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCGACGTCCGGGTGCGCtt   >  1:580046/1‑150 (MQ=255)
                                                                                                                                              gACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCGACGTCCGGGTGCGCtt   >  2:1155938/1‑150 (MQ=255)
                                                                                                                                               aCTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCGACGTCCGGGTGCGCttt  <  2:1260261/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
ACCAAACATGGCGAACAAATGTTGCAAACTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCGACGTCCGGGTGCGCTTT  >  NC_000913/2620071‑2620362

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: