Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,828,980 A→C 5.7% Q12H (CAA→CAC)  xanP → xanthine permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,828,9800AC 5.7% 1475.5 / 6.9 544Q12H (CAA→CAC) xanPxanthine permease
Reads supporting (aligned to +/- strand):  ref base A (277/221);  new base C (7/39);  total (284/260)
Fisher's exact test for biased strand distribution p-value = 1.01e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 566 total aligned reads displayed.

ATATAAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTGGG  >  NC_000913/3828831‑3829129
                                                                                                                                                      |                                                                                                                                                      
 tataAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAAc                                                                                                                                                       >  2:8204607/1‑150 (MQ=255)
  ataAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAAcc                                                                                                                                                      <  2:667493/150‑1 (MQ=255)
  ataAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAAcc                                                                                                                                                      <  1:7305122/150‑1 (MQ=255)
   taAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCg                                                                                                                                                     <  1:7024950/150‑1 (MQ=255)
   taAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCg                                                                                                                                                     >  2:3672441/1‑150 (MQ=255)
   taAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCg                                                                                                                                                     >  2:10956449/1‑150 (MQ=255)
     aaaaaGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGt                                                                                                                                                   >  1:6499454/1‑150 (MQ=255)
     aaaaaGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGt                                                                                                                                                   >  2:2964755/1‑150 (MQ=255)
      aaaaGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCCACCCGtt                                                                                                                                                  >  2:97121/1‑150 (MQ=255)
      aaaaGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGtt                                                                                                                                                  <  2:7464638/150‑1 (MQ=255)
             cATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCatt                                                                                                                                           <  2:9832846/150‑3 (MQ=255)
              aTTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGAc                                                                                                                                          <  2:11139277/150‑1 (MQ=255)
              aTTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGAc                                                                                                                                          <  2:5791418/150‑1 (MQ=255)
              aTTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGAc                                                                                                                                          >  1:9549539/1‑150 (MQ=255)
                  aTCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCaa                                                                                                                                      >  1:10524905/1‑150 (MQ=255)
                        tttAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGc                                                                                                                                >  2:2628984/1‑150 (MQ=255)
                        tttAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAAAAGc                                                                                                                                >  1:4459910/1‑150 (MQ=255)
                           aaaTATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGaa                                                                                                                             <  2:7004758/150‑1 (MQ=255)
                           aaaTATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACACAAAACAGCGac                                                                                                                             >  1:9900086/1‑149 (MQ=255)
                               atTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGACAGAAAATGCGCAACCGGGTGAGCACACTCAAAACAGCGAACTGa                                                                                                                         >  1:5663648/1‑150 (MQ=255)
                                tttAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGAt                                                                                                                        <  2:8335668/150‑1 (MQ=255)
                                tttAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACAGAt                                                                                                                        >  2:4937138/1‑150 (MQ=255)
                                 ttAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGAtt                                                                                                                       <  2:4606699/150‑1 (MQ=255)
                                    gATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCCCGGGAGTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCCGGTGCGCAGACTCAAAACAGCGAACGGATTTAc                                                                                                                    >  2:7582403/1‑150 (MQ=255)
                                       gAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCCCGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCCGTTGCGCAGACTCAAAACAGCGAACTGATTTACCgg                                                                                                                 >  1:805607/1‑149 (MQ=255)
                                       gAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGt                                                                                                                 <  2:9173782/150‑1 (MQ=255)
                                            tgtgGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCCAACGCTTGCGCAGACTCAAAACAGCGAACAGATTTAACGTCttta                                                                                                            >  1:4838582/1‑148 (MQ=255)
                                             gtgGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGaa                                                                                                           <  2:1560303/150‑1 (MQ=255)
                                               ggCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGa                                                                                                         >  2:6564456/1‑150 (MQ=255)
                                               ggCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCCGTTGCGCAGACTCAAAAAAGCGAACATATTTAACGTCATGAAGa                                                                                                         >  1:10028600/1‑150 (MQ=255)
                                                gCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCCACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCGTGAAGAt                                                                                                        >  2:3149679/1‑150 (MQ=255)
                                                gCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGAt                                                                                                        >  2:8601997/1‑150 (MQ=255)
                                                 cAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATc                                                                                                       >  2:93600/1‑150 (MQ=255)
                                                 aaaagttttGCTTTTGACTTTCGGTCACCTAAACGGCCGGCTTTGCTACCGGGGTTTTT‑TTTTCGTTGTCTTTTACCACCCTCGAGTCACAAATTGCGCAACGGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATc                                                                                                       <  1:1219507/144‑1 (MQ=255)
                                                      gTTTGCTTTTCTCTTGCGGTCAGCTAAACTGCCCGCTTTGCTACCTCGGGTTGG‑TTGTCGATGTCTGTTTCCACCCTCTAGACATAAAAGGCGCAACGGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTccg                                                                                                  <  1:6725211/150‑1 (MQ=255)
                                                       tttGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTccgc                                                                                                 >  2:1613395/1‑150 (MQ=255)
                                                           cTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACAGATTTACCGTCTTGAAGAACGTCCGCCGCt                                                                                             >  1:10250371/1‑150 (MQ=255)
                                                              ttCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTcc                                                                                          >  2:1070852/1‑150 (MQ=255)
                                                              ttCTCTTTCGGTCAGATAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCAACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTcc                                                                                          <  2:10400582/150‑1 (MQ=255)
                                                                   tttCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCACCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCaaa                                                                                     <  2:10060563/150‑1 (MQ=255)
                                                                   tttCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCGTCCTCaaa                                                                                     >  1:225763/1‑150 (MQ=255)
                                                                     tCGGTCAGCTAAAATGCCCGCTTTGCTACCTCGGGATTG‑TTTTCTATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAc‑c                                                                                  <  2:206650/150‑1 (MQ=255)
                                                                       ggTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCCACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCt                                                                                >  1:8334758/1‑150 (MQ=255)
                                                                       ggTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCGTCCTCAAAC‑CCt                                                                                >  1:7875490/1‑150 (MQ=255)
                                                                         tCAGCTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCTGt                                                                              >  2:275278/1‑150 (MQ=255)
                                                                         tCAGCTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCAGt                                                                              >  1:1081177/1‑150 (MQ=255)
                                                                           aGCTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCTGttt                                                                            <  2:1943196/150‑1 (MQ=255)
                                                                            gCTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCGTCCTCAAAC‑CCTGTTTg                                                                           >  1:3595077/1‑150 (MQ=255)
                                                                             cTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACAGTCTTGAAGATCGTCCGCCGCTTCATCAAAC‑ACTGTTTgc                                                                          >  2:4462124/1‑150 (MQ=255)
                                                                              taaacggCCCGCTTTGCTACC‑CGGGGTTGTTTTTCGTGGTCTGTGACCACCCTGGAGTCAGAAATTGCGCACCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCTGTTTgcc                                                                         <  1:5496042/144‑1 (MQ=255)
                                                                              tAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCTGTTTgcc                                                                         >  2:8988995/1‑150 (MQ=255)
                                                                               aaaaTGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCTGTTTgccg                                                                        <  2:8015841/150‑1 (MQ=255)
                                                                                aattGCCCGCTTTGCTACCTCGGGTTTG‑TTTTGGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCTGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCTGTTTgccgc                                                                       <  2:394043/147‑1 (MQ=255)
                                                                                  aTGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCACAC‑CCTGTTTGCCGCCt                                                                     <  2:6170606/150‑1 (MQ=255)
                                                                                  aTGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCTGTTTGCCGCCt                                                                     >  1:5996531/1‑150 (MQ=255)
                                                                                     cccGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGACCTGAGTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCGGTGTGCCGCCTGTc                                                                  >  2:1803263/1‑150 (MQ=255)
                                                                                     cccGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCTGTTTGCCGCCTGTc                                                                  <  2:4205842/150‑1 (MQ=255)
                                                                                      ccGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCTGTTTGCCGCCTGTCa                                                                 >  2:6470767/1‑150 (MQ=255)
                                                                                      ccGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAAAAGCGAACTGATTTACCGTCTTGAATATAGTACGCCGCATCCACAAAC‑CCAGGTTTCCGCCAGtaa                                                                 >  1:3188389/1‑148 (MQ=255)
                                                                                       cGCTTTGCTACCAGGGTTTTG‑TTTGCGATGTCTGTTTCCACCCTCGAGTCAGAAATTGCGCACCCGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCTGTTTGCCGCCTGTCAg                                                                <  1:8755266/150‑1 (MQ=255)
                                                                                        gCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCTGTTTGCCGCCTGTCAGc                                                               <  2:4012291/150‑1 (MQ=255)
                                                                                            tGCTACCACGGGTTTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCTGTTTGCCGCCTGTCAGCATct                                                           <  2:241708/150‑1 (MQ=255)
                                                                                              cTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGAGTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCTGTTTGCCGCCTGTCAGCATctgc                                                         >  2:9529966/1‑150 (MQ=255)
                                                                                              cTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCCGTTGCGCACAATCAAAAAAGCGAAAATATTTTCCGTCTTGAAGATCGTCCGCCCCGTCAACAAAA‑CCGGTTTGGCGCCAGTCAGCATcagc                                                         >  1:4685103/1‑150 (MQ=255)
                                                                                              cTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCCGGTGCGCAGACACAAAACACCCAACTGATTTACCCGCGTGAAGATCGTCCCCCGCGTCCACAAAC‑CCCGTTTTCCGCCAGTCAGCATcggc                                                         >  1:2727811/1‑150 (MQ=255)
                                                                                                aCCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCTGTTTGCCGCCTGTCAGCATctgctg                                                       <  2:272553/150‑1 (MQ=255)
                                                                                                 ccACGGTATTT‑TTTTCTATTTCTTTTTCCACCCTCGAGTCATAAAATGCGCAACCGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCTGTTTGCCGCCTGTCAGCATCTGCTgg                                                      <  1:5799558/150‑1 (MQ=255)
                                                                                                  cACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAATTGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCTGTTTGCCGCCTGTCAGCATCTGCTGGc                                                     <  1:7629730/150‑1 (MQ=255)
                                                                                                  cACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGGGCCACCGGGTGCGCAAGAACAAAAAAGCGAACTTTTTTTCCGTCATGAATATCGGCCGCCGCGGCCTCAAAA‑CCGGTTTGCCGCCGGTCAGCAACAGCTCGc                                                     >  1:9258893/1‑150 (MQ=255)
                                                                                                    cGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCGACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCTGTTTGCCGCCTGTCAGCATCTGCTGGCGa                                                   <  2:1309782/150‑1 (MQ=255)
                                                                                                    cGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGGTGCGCAGACTCAAAAAAGCGAACGGATTTACCGACTTGAAGATCGTCCGCCGCGTCCTCAAACACGGGTTTGCC‑CCTGTCAGCATCTTCTGGCGa                                                   >  1:8431194/1‑150 (MQ=255)
                                                                                                       gATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCCACCCGTTGCGCAGACACAAAAAAGCGAACTGATTTAACGTCTTGAAGATCGTCCGCCGCTTCCTCAAAA‑CCAGTTTGCCGCCGGTCAGCATCAGCTGGCGATGt                                                >  1:1820565/1‑150 (MQ=255)
                                                                                                       gATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCACCGTGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGt                                                <  1:7152324/150‑1 (MQ=255)
                                                                                                       gATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAAACGGTTGCGCAGACAAAAAAAAGCGAACTTATTTAACGTCATGAATATAGTCCGCCGCTTCCTCAAAA‑CCAGTTTTCCGCCTGTAAGCATCATCTAGCGATGt                                                >  1:8235903/1‑150 (MQ=255)
                                                                                                         ttG‑TTTTCGATGTCTGTTTCCACCCTCGAGACACAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTc                                              <  1:3407313/150‑1 (MQ=255)
                                                                                                                tCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCCGGTGCGCAGACTCAAAAAAGCGAACAGATGTAACGTCTTGTAGATCGTGCGCCGCGTCCACAAAC‑CCAGTTTGCCGCCTGTAAGCATCTGCTGGCGGGGTTAGTTggg                                        >  1:4361870/1‑148 (MQ=255)
                                                                                                                tCGAGGTCTGTTCCCCCCCTCGAGTCAGAAAATGCGCACCCGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCg                                        <  1:1338470/150‑1 (MQ=255)
                                                                                                                 cGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGGTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCgg                                       <  1:9020049/150‑1 (MQ=255)
                                                                                                                   aTGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTg                                     >  2:9139956/1‑150 (MQ=255)
                                                                                                                    tGTCTGTTACCACCCTCGAGTCAGAAAATGCGCACCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGa                                    <  1:5141644/150‑1 (MQ=255)
                                                                                                                     gTCTGTTTCCACCCTCGAGTCAGAAAATGCACAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGAt                                   <  2:5265672/150‑1 (MQ=255)
                                                                                                                     gTCTGTTTCCACCCTCGAGACAGAAAATGCGCACCCGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGAt                                   <  1:961926/150‑1 (MQ=255)
                                                                                                                      tCTGTTTCCACCCTCGAGTCAGAAAATGCGCCACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATc                                  >  2:1164492/1‑150 (MQ=255)
                                                                                                                      tCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATc                                  >  2:3293247/1‑150 (MQ=255)
                                                                                                                      tCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCCGTTGCGCAGACTCAAAACAGCGAAAGGATTTAACGTAGTTAAGATAGTCCGCCGCTTCCTCAAAC‑CCTGTTTGGCGCCGGTCAGCAACAGCTGGCGAAGTTCGTTGCGGTgtta                                  >  1:6671423/1‑147 (MQ=255)
                                                                                                                            tCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCa                            <  2:5490633/150‑1 (MQ=255)
                                                                                                                             ccccccGCGAGGCAGAAACGGCGCACCCGGGGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAg                           <  1:284090/147‑1 (MQ=255)
                                                                                                                             ccACCCTCCAGCAAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAg                           <  1:3781933/150‑1 (MQ=255)
                                                                                                                               aCCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAgcg                         <  1:5268995/150‑1 (MQ=255)
                                                                                                                                        tCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATc                >  2:7835512/1‑150 (MQ=255)
                                                                                                                                         cAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCt               <  1:2408886/150‑1 (MQ=255)
                                                                                                                                          aGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTg              <  2:7141408/150‑1 (MQ=255)
                                                                                                                                           gAAAATGCGCAACCGGTTGCGCAGACTCAAAAAAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCATCCTCAAAC‑CCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGc             >  2:3826139/1‑150 (MQ=255)
                                                                                                                                             aaaTGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCa           >  2:732292/1‑150 (MQ=255)
                                                                                                                                             aaaTGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCa           <  2:8785233/150‑1 (MQ=255)
                                                                                                                                             aaaTGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCa           >  2:3456871/1‑150 (MQ=255)
                                                                                                                                                    cAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTg    >  1:3979061/1‑150 (MQ=255)
                                                                                                                                                    cAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTg    >  2:5988584/1‑150 (MQ=255)
                                                                                                                                                     aaCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTgg   >  2:4773616/1‑150 (MQ=255)
                                                                                                                                                     aaCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCGTCCTCAAAC‑CCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTgg   <  2:10560156/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
ATATAAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTG‑TTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAC‑CCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTGGG  >  NC_000913/3828831‑3829129

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑