Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,620,218 T→A 7.1% intergenic (‑58/+28) 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,2180TA 7.1% 319.3 / 12.1 112intergenic (‑58/+28)uraA/uppuracil permease/uracil phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (59/45);  new base A (3/5);  total (62/50)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.38e-01

Only 100 of 116 total aligned reads displayed.

GCACCAAACATGGCGAACAAATGTTGCAAACTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCGACGTCCGGGTGCGCTTTTTCC  >  NC_000913/2620069‑2620366
                                                                                                                                                      |                                                                                                                                                    
gCACCAAACATGGCGAACAAATGTTGCAAACTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCtt                                                                                                                                                      >  2:1775337/1‑150 (MQ=255)
 cACCAAACATGGCGAACAAATGTTGCAAACTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCAGACTCTTa                                                                                                                                                     >  2:1478567/1‑150 (MQ=255)
   ccAAACATGGCGAACAAATGTTGCAAACTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAAGCCGACTCTTaa                                                                                                                                                    >  1:574246/1‑150 (MQ=255)
    cAAACATGGCGAACAAATGTTGCAAACTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAg                                                                                                                                                  >  2:1640322/1‑150 (MQ=255)
      aaCATGGCGAACAAATGTTGCTAACTAAGCGTGATTGTCTGTAAAATTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTTTTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTc                                                                                                                                                <  2:1111272/150‑1 (MQ=255)
      aaCATGGCGAACAAATGTTGCAAAATAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTAACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTc                                                                                                                                                <  2:1105489/150‑1 (MQ=255)
        cATGGCGAACAAATGTTGCAAACTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCgg                                                                                                                                              >  1:1560787/1‑150 (MQ=255)
         aTGGCGAACAAATGTTGCAAACTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTTTAAGAGTCGGc                                                                                                                                             <  1:1519755/150‑1 (MQ=255)
           ggCGAACAAATGTTGCAAACTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCtt                                                                                                                                           >  2:71265/1‑150 (MQ=255)
              gAACAAATGTTGCAAACTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTaa                                                                                                                                        <  2:61333/150‑1 (MQ=255)
                aCAAATGTTGCAAAATAAGCGGGATTGTCTGTAAAAGTGGCGTGCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTTTAAGAGTCGGCttttttt                                                                                                                                      <  2:1160749/150‑5 (MQ=255)
                aCAAACGTTCCAAACTAAGCGGGATTGGCAGAAAAAGTGGCGGTCTTTCACTCACCCAGATAGCACGGCGAGTCATAGTAATATCCTCTGTATTCTGTGGTATAGGCGCTGTACTCAAAAATAAGCCGACTTTAAGAGTCGGCttttttt                                                                                                                                      <  1:271796/150‑5 (MQ=255)
                  aaaTGTTGCAAACTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTTTAAGAGTCGGCttttttttg                                                                                                                                    <  1:453648/150‑7 (MQ=255)
                   aaTGTTGCATATTTCGGGGTATTTTTTGTAAATGTGGTTGTATTTCACTCACCCCGATACCCCGGCGCGTCATTGTATTATCCTTTTTATTATGTGTTATTGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAattatt                                                                                                                                   <  2:156369/150‑1 (MQ=255)
                   aaTGTTGCAAACTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTTTAAGAGTCGGCttttttttga                                                                                                                                   <  1:892834/150‑8 (MQ=255)
                       ttGCAAACTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTa                                                                                                                               >  1:1530688/1‑150 (MQ=255)
                       ttGCAAACTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTAACTCTCGCAGATAGCTCGGAGCGTGATAGTATTTTCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTGTTAAAGTCGGGTTTAATTAtttttt                                                                                                                               >  1:104036/1‑149 (MQ=255)
                              cTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTa                                                                                                                        >  2:780150/1‑150 (MQ=255)
                              cTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTa                                                                                                                        >  2:1450648/1‑150 (MQ=255)
                                   cgtgATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCAAGTCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCTCTTTACTCAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGt                                                                                                                  <  1:144025/147‑1 (MQ=255)
                                   cGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATATCACGGCGTGTAATAGTATTATCATCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCg                                                                                                                   <  1:118310/150‑1 (MQ=255)
                                          gTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCaaa                                                                                                            >  1:727016/1‑150 (MQ=255)
                                            cTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGa                                                                                                          >  1:1622639/1‑150 (MQ=255)
                                            cTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGa                                                                                                          >  1:552529/1‑150 (MQ=255)
                                            cTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGa                                                                                                          >  1:600402/1‑150 (MQ=255)
                                            cTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGa                                                                                                          >  2:513980/1‑150 (MQ=255)
                                            cTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGa                                                                                                          <  2:1608066/150‑1 (MQ=255)
                                            cTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTATTTATTTCGTACCAAAGa                                                                                                          >  2:285087/1‑150 (MQ=255)
                                             tGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGAt                                                                                                         >  1:1168747/1‑150 (MQ=255)
                                              gTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGAtt                                                                                                        >  1:992327/1‑150 (MQ=255)
                                               tAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGAtt                                                                                                        >  2:1262215/1‑150 (MQ=255)
                                                 aaaGTGGCGGTCTTTCACTCACCCCTATAGCACGGGGTGTCATATTATTATCCTCTTTATTTTGTGTTATAGGCGCTTTAATAAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTg                                                                                                     <  2:1653469/150‑1 (MQ=255)
                                                  atgTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGt                                                                                                    <  2:1642007/148‑1 (MQ=255)
                                                   aGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTc                                                                                                   >  2:712344/1‑150 (MQ=255)
                                                            cTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATc                                                                                          >  1:80385/1‑150 (MQ=255)
                                                                cactcacCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACAAAAGATTTTGTCACCGGCATCGccc                                                                                      >  2:545356/1‑149 (MQ=255)
                                                                 actcacCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACGATTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGa                                                                                     >  2:293411/1‑150 (MQ=255)
                                                                  ctcacCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAgg                                                                                   >  2:1707357/1‑150 (MQ=255)
                                                                   tcacCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAgg                                                                                   <  2:1301372/150‑1 (MQ=255)
                                                                    cacCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGc                                                                                  >  1:126295/1‑150 (MQ=255)
                                                                     acCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGcc                                                                                 <  1:656264/150‑1 (MQ=255)
                                                                      ccccGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGGTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGccc                                                                                >  1:1139572/1‑150 (MQ=255)
                                                                       cccGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCg                                                                               >  2:212608/1‑150 (MQ=255)
                                                                       cccGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCg                                                                               <  1:212608/150‑1 (MQ=255)
                                                                         cGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGa                                                                             <  1:848725/150‑1 (MQ=255)
                                                                         cGATAGCACGGCGCGTCATAGTATTATCCTCTGAATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGa                                                                             <  2:1279360/150‑1 (MQ=255)
                                                                           ttAGCACGGCGCGTCATCTTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGaat                                                                           <  2:731584/149‑1 (MQ=255)
                                                                              gCACGGCTCGTCATAGTTTTATCCTCTGTATTATTTTTAATAGGGGCTTTACTAAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGaataat                                                                        <  2:870445/150‑1 (MQ=255)
                                                                              gCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGaataat                                                                        >  1:1053561/1‑150 (MQ=255)
                                                                              gCACGGCGCGTCATAGTATTATCCTCTGTATTAAGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGaataat                                                                        >  1:222495/1‑150 (MQ=255)
                                                                               cACGGCGCTTAATAGTATTATCCTATTTATTATGTGTAAAAGGCTCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTAATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATg                                                                       <  1:60970/150‑1 (MQ=255)
                                                                               cACGGCGCGTCATAGTATTATCCTCTGTATTATTTTTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATg                                                                       <  1:1667432/150‑1 (MQ=255)
                                                                               cACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATg                                                                       >  1:586874/1‑150 (MQ=255)
                                                                               cACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATg                                                                       <  1:1774206/150‑1 (MQ=255)
                                                                                aCGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGt                                                                      >  1:837275/1‑150 (MQ=255)
                                                                                   gcgcGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATc                                                                   <  2:815888/150‑1 (MQ=255)
                                                                                       gTCATAGTATTATCCTCTGTATTATGTGTTATAGTCGCTTTACTCAAAAAAAATCCGACTCTTAAAGTCTGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCTTg                                                               <  1:1180711/150‑1 (MQ=255)
                                                                                                  atCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTAACAAAGATTTTGTAACCGGAATCGACGAGGCCAGGAATAATGTATACGTGCTCTTTAAGTc                                                    >  1:1043954/1‑150 (MQ=255)
                                                                                                    ccTCTGTATTATGTGTTATCGGCGCTTTACACAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTccc                                                  <  2:1560601/150‑1 (MQ=255)
                                                                                                     ctctGTTTTATGTGTTATTGGCGCTTAACTCAAAAAAAAGCCGACTCTTAAAGGCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCt                                                 <  1:255384/150‑1 (MQ=255)
                                                                                                     ctctGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCt                                                 >  1:1726541/1‑150 (MQ=255)
                                                                                                      tctGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTg                                                >  1:1479854/1‑150 (MQ=255)
                                                                                                      tctGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTg                                                >  1:1249304/1‑150 (MQ=255)
                                                                                                         gTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATc                                             <  2:1050769/150‑1 (MQ=255)
                                                                                                          tattatGCGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAACGTATCCGTGCTCGTTCAGTCCCTGATCa                                            <  1:1449454/150‑1 (MQ=255)
                                                                                                             tatGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATc                                         <  1:102736/150‑1 (MQ=255)
                                                                                                              atGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTAGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCg                                        >  1:1409265/1‑150 (MQ=255)
                                                                                                                gtgtTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGAt                                      >  1:1422585/1‑150 (MQ=255)
                                                                                                                gtgtTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGAt                                      >  2:196702/1‑150 (MQ=255)
                                                                                                                 tgtTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATg                                     <  2:48894/150‑1 (MQ=255)
                                                                                                                 tgtTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATg                                     <  1:503523/150‑1 (MQ=255)
                                                                                                                         gcgcTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATAc                             <  2:727016/150‑1 (MQ=255)
                                                                                                                         gcgaTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTAATTCTTTATTTCGTACCAAAGTATTTGTCACCGGCATCGCCGAGGCCCGGAATATTGTATCCGAGGTCGTTCAGTCCCTGATCAATCGATGCGGTATAc                             <  1:318973/150‑1 (MQ=255)
                                                                                                                           gcTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAg                           >  1:1020014/1‑150 (MQ=255)
                                                                                                                           gcTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAg                           <  1:343584/150‑1 (MQ=255)
                                                                                                                            cTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGt                          <  1:1014239/150‑1 (MQ=255)
                                                                                                                              ttACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTc                        <  2:1464899/150‑1 (MQ=255)
                                                                                                                              ttACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTc                        <  2:1314025/150‑1 (MQ=255)
                                                                                                                               tACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCg                       <  2:1530688/150‑1 (MQ=255)
                                                                                                                               tACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCg                       <  2:837275/150‑1 (MQ=255)
                                                                                                                                aatCAAAAAAACGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCGa                      <  2:1792346/148‑1 (MQ=255)
                                                                                                                                 cTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCGAc                     >  1:440478/1‑150 (MQ=255)
                                                                                                                                 cTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTAGACAGTTCGAc                     >  1:1399190/1‑150 (MQ=255)
                                                                                                                                 cTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCATGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCGAc                     >  1:1533186/1‑150 (MQ=255)
                                                                                                                                 acgaaaaaaaatGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCGa                      <  1:341424/146‑1 (MQ=255)
                                                                                                                                   cAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCAGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCGACGt                   >  2:1674972/1‑150 (MQ=255)
                                                                                                                                     aaaaaaaaGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCGACGTc                  >  2:1187311/1‑150 (MQ=255)
                                                                                                                                       attaaaGCCGACTTTAAGAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCAGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCGACGTCCg                >  2:93816/4‑150 (MQ=255)
                                                                                                                                       aaaaaaGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCGACGTCCg                <  1:712344/150‑1 (MQ=255)
                                                                                                                                          aaaGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCGACGTCCGGGt             <  2:600402/150‑1 (MQ=255)
                                                                                                                                           aaGCCGACTTTAAGAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCGACGTCCGGGTg            >  2:714058/1‑150 (MQ=255)
                                                                                                                                           aaGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCGACGTCCGGGTg            >  2:498614/1‑150 (MQ=255)
                                                                                                                                           aaGCCGACTCTTAAAGTCGGCATTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCTAGGCCCGGAATAATGTATCCGTGCTAGTTCAGTCCCTGATCAATCGATGCGGTATACAGGTCGACGTCCGGGTg            >  2:309197/1‑150 (MQ=255)
                                                                                                                                             gCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCGACGTCCGGGTgcg          >  2:1166546/1‑150 (MQ=255)
                                                                                                                                               cGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGAAGCGGTATACAGTTCGACGTCCGGGTGCGCt        >  2:1373648/1‑150 (MQ=255)
                                                                                                                                                   tctTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCGACGTCCGGGTGCGCttttt    >  1:806843/1‑150 (MQ=255)
                                                                                                                                                    ctTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCGACGTCCGGGTGCGCTTTTTc   >  1:1281199/1‑150 (MQ=255)
                                                                                                                                                    ctTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCGACGTACGGGTGCGCTTTTTc   >  2:944151/1‑150 (MQ=255)
                                                                                                                                                     ttAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCGACGTCCGGGTGCGCTTTTTcc  <  2:984400/150‑1 (MQ=255)
                                                                                                                                                     ttAAAGTCGGCTTTAATTATTTTTATGCGTTATTTCGTACCAAATATTTTGTCACAGGCATCGCCGAGGCCCCGAATAAAGTATCAGTGCTCTTTCAGTCCCTGATCACTCGATGCGGTATACAGTTCTACGTCCGTGTCCGCTTTTTcc  >  1:1083672/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
GCACCAAACATGGCGAACAAATGTTGCAAACTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCGACGTCCGGGTGCGCTTTTTCC  >  NC_000913/2620069‑2620366

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: