Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,828,980 A→C 7.5% 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 7.5% 2248.5 / 27.9 842Q12H (CAA→CACxanPxanthine permease
Reads supporting (aligned to +/- strand):  ref base A (424/331);  new base C (14/71);  total (439/403)
Fisher's exact test for biased strand distribution p-value = 9.70e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 867 total aligned reads displayed.

ATATAAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTGGG  >  NC_000913/3828831‑3829129
                                                                                                                                                     |                                                                                                                                                     
atataAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCaa                                                                                                                                                       <  2:8718070/150‑1 (MQ=255)
  ataAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAAcc                                                                                                                                                     <  2:370291/150‑1 (MQ=255)
   taAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCg                                                                                                                                                    >  1:1576686/1‑150 (MQ=255)
   taAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCg                                                                                                                                                    >  1:3916748/1‑150 (MQ=255)
   taAAAAAGTGCATTTATCTCAGTTAAATATTTAGATGAGGATGTGTCAAACGTTTGCTTTTATATTTAGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCTATGTCTGTTTCCACACTCGAGTCAGAAAATGCGCAACCg                                                                                                                                                    >  1:4992737/1‑150 (MQ=255)
       aaaGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTg                                                                                                                                                <  1:2093529/150‑1 (MQ=255)
       aaaGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTg                                                                                                                                                <  2:1413730/150‑1 (MQ=255)
        aaGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTgc                                                                                                                                               <  1:3671326/150‑1 (MQ=255)
         aGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTgcg                                                                                                                                              >  2:1546554/1‑150 (MQ=255)
          gTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTgcgc                                                                                                                                             <  2:6765926/150‑1 (MQ=255)
             cATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGa                                                                                                                                          >  2:4816314/1‑150 (MQ=255)
              aTTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGAc                                                                                                                                         >  2:3931707/1‑150 (MQ=255)
              aTTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGAc                                                                                                                                         >  1:1349054/1‑150 (MQ=255)
              aTTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCCGTTGCGCAGAc                                                                                                                                         >  1:6374350/1‑150 (MQ=255)
                   tctcATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCaaa                                                                                                                                    >  1:3376426/1‑150 (MQ=255)
                     tcATTTAAATATATAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAAc                                                                                                                                  >  2:7067096/1‑150 (MQ=255)
                         ttAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCg                                                                                                                              >  2:6942158/1‑150 (MQ=255)
                         ttAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCCGTTGCGCAGACTCAAAACAGCg                                                                                                                              >  1:177825/1‑150 (MQ=255)
                            aaTATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAAc                                                                                                                           <  2:309976/150‑1 (MQ=255)
                                tttAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGTGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGAt                                                                                                                       <  2:5960496/150‑1 (MQ=255)
                                 ttAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCTAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGAtt                                                                                                                      <  2:8390516/150‑1 (MQ=255)
                                 ttAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGAtt                                                                                                                      <  2:7687152/150‑1 (MQ=255)
                                   aGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAAGGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTa                                                                                                                    <  2:1865560/150‑1 (MQ=255)
                                   aGATGAGGATGTGGCAAACGTTTGCTGTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCCACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTa                                                                                                                    >  2:7333652/1‑150 (MQ=255)
                                     aTGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCCGACTCAAAACAGCCAACAGATTTaac                                                                                                                  >  1:1020948/1‑148 (MQ=255)
                                       gAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACACAAAACAGCGAACAGATTTACCGt                                                                                                                >  1:7763419/1‑150 (MQ=255)
                                          tttgtgGCAATCGTTTGCTTTTATCTTTCGGTCAGCTAAAATGCCCGCTTTTCTACCACGGTATTGTTTTGTTTGTCTGTTTCCACACTCGTGTCAGAAAATGCGCATCGGGTTTCGCAGACTCAAAACAGCGAACTGATTTACCGTCtt                                                                                                             <  2:7935186/148‑1 (MQ=255)
                                          gATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCtt                                                                                                             <  2:2279076/150‑1 (MQ=255)
                                          gATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCCCAACCGGTTGCGCAGACTCAAAAAAGCGAACAGATTTAACGTCtt                                                                                                             >  1:5979604/1‑150 (MQ=255)
                                          gATGTGGCAAACGTTTGCTTTTCTCTTACGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGACGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCtt                                                                                                             <  2:8560606/150‑1 (MQ=255)
                                            tgtgGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGa                                                                                                           <  2:6070673/150‑1 (MQ=255)
                                              tgGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATATACCGTCTTGAAg                                                                                                         >  2:5733146/1‑150 (MQ=255)
                                               ggCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGa                                                                                                        <  2:5270878/150‑1 (MQ=255)
                                                gCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAAAAGCGAACTGATTTACCGTCTTGAAGAt                                                                                                       >  2:5544127/1‑150 (MQ=255)
                                                gCAAACGTTTGCTTTTCTCTTTCGGTCAGCAAAAATGCCCGCTTTGCTACCACGGGATTGTTTTAGATGTCTGTTTCCACCCTCGAGTCATAAAAAGCGCAAACCGTTGCGCAGACACAAAAAAACCAACAGATTTACACGCTTGAAGAt                                                                                                       >  1:8109959/1‑150 (MQ=255)
                                                 cAACCGTTTGCTTTTCTATTTCGGTCAGCTAAAATGCCCGCTTTGCTACCAGGGTTTTGTTTTCGCTGTCTTTGTCCACCCTCGAGACAGAAAAGGCGCAACCGGTTGCCCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATc                                                                                                      <  1:4201397/150‑1 (MQ=255)
                                                 cAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATc                                                                                                      <  2:1091819/150‑1 (MQ=255)
                                                  aaaCGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCg                                                                                                     >  2:3654684/1‑150 (MQ=255)
                                                  aaaCGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACGGATTTACCGTCTTGAAGAtag                                                                                                     >  1:5364493/1‑148 (MQ=255)
                                                    aCGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGAACGTc                                                                                                   >  1:6917507/1‑150 (MQ=255)
                                                     tttttGCTTTTCTCTGTGGGTCTGCTAAACGGCCCTCTTTGCTACCTGGGGTTTGTTTTCGTTGTGTGTTTCCACCCTCGAGTCAGAAAAGGCGCACCGGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTcc                                                                                                  <  1:7618764/148‑1 (MQ=255)
                                                       tttGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACAGATTTACCGTCTTGAAGATAGTccgc                                                                                                >  2:3202348/1‑150 (MQ=255)
                                                           cTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCt                                                                                            <  2:3554127/150‑1 (MQ=255)
                                                            ttttCTCTTTCGGTCAGCTAAAATGACCGCTATGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCtt                                                                                           >  2:2441389/1‑150 (MQ=255)
                                                             tttCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTc                                                                                          >  2:2774803/1‑150 (MQ=255)
                                                             tttCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTc                                                                                          >  2:3987662/1‑150 (MQ=255)
                                                                 tctTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCa                                                                                      >  1:3036250/1‑150 (MQ=255)
                                                                   tttCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTTCGCAGACTCAAAACAGCGAACTGATTTAACGTCTTGAAGATATTCCGCCGCTTACTCaaa                                                                                    >  1:8833679/1‑150 (MQ=255)
                                                                    ttCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAc                                                                                   <  2:5411264/150‑1 (MQ=255)
                                                                       ggTCAGCTAAAATGCCCGCTTTGCTACCACGGTATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCt                                                                                <  2:1680364/150‑1 (MQ=255)
                                                                       ggTCAGCTAAAATGCCCGCTTTGCTACCACGGTATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCt                                                                                <  2:7832578/150‑1 (MQ=255)
                                                                       ggTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCCGTTGCGCAGACTCAAAACAGCGAACAGATTTACCGGCTTGAAGATCGTCCGCCGCTTCCACAAACCCt                                                                                >  1:4032403/1‑150 (MQ=255)
                                                                        gTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTg                                                                               >  2:9002287/1‑150 (MQ=255)
                                                                         tCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACACTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCCAACGGATTTAACGTAGTGAATATAGTCCGCCGCTTCCTCAAAACCtat                                                                              >  1:6683863/1‑148 (MQ=255)
                                                                            gatAAACTGCCCGCTTTGCTACCTGGTGTTTTTTTTCTTTGTCTGTTCCCACCCTCGAGTCAGAAAATGCGAAACGGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTg                                                                           <  1:8603492/148‑1 (MQ=255)
                                                                                     cccGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTc                                                                  >  2:7486279/1‑150 (MQ=255)
                                                                                      ccGCTTTGCTACCCCGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCa                                                                 <  2:4080947/150‑1 (MQ=255)
                                                                                      ccGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCa                                                                 <  2:2162546/150‑1 (MQ=255)
                                                                                      ccGCTTGGCTACCAGGGTTTTGTTTTCGATGTGTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCa                                                                 <  1:4539735/150‑1 (MQ=255)
                                                                                       cGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCACCCGGTTGCGCAGACACAAAAAAGCGAACTGATTTAACGTCTTGAAGATCGTCCGCCGCGTCCTCAAACCCAGTTTGCCGCCGGTCAg                                                                >  1:2898812/1‑150 (MQ=255)
                                                                                       cGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGGTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCCCTTCCTCAAAACCGGTTTTGCGCCTGTCAg                                                                >  1:958730/1‑150 (MQ=255)
                                                                                       cGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGAGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAg                                                                >  2:4473601/1‑150 (MQ=255)
                                                                                        gCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAAAAGCGAACAGATTTACCGTATTGAATATAGTCCGCCGCTTCCTCAAACCCAGTTTGCCGCCTGTCAGc                                                               >  1:3533315/1‑150 (MQ=255)
                                                                                        gCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGAATCAAAACAGCGAACAGATTTAACGTCATGAAGATAGGCCGCCCCATCCTCAAACCCAGTTTGCCGCCAGTccac                                                               >  1:8952254/1‑147 (MQ=255)
                                                                                         ttttGCTACCAGGGTTTTTTTTTGGTTTTCTGTTACAACCCTCTTGTCATAATTTGCGAAACCGTTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCa                                                              <  1:229994/149‑1 (MQ=255)
                                                                                            tGCTACCTCGGTTTTGTTTTCGTTGTCTGTTTCCACCCTCGAGTCATAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATct                                                           <  1:5535766/150‑1 (MQ=255)
                                                                                             gCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATctg                                                          >  2:5090621/1‑150 (MQ=255)
                                                                                              cTACCAGGGGTTTGTTTTCTATTTCTGTTTCCACCCTCGAGCCAGAAAATGCGCACCGGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATctgc                                                         <  1:4306088/150‑1 (MQ=255)
                                                                                              cTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATctgc                                                         <  2:1792262/150‑1 (MQ=255)
                                                                                              cTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATctgc                                                         >  2:6258720/1‑150 (MQ=255)
                                                                                               tACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCCGTTGCGCAGACTCAAAAAAGCCAACGGATTTACCGTCTTGAAGATCGTCCGCCCCTTCCTCAAACCCAGTTTTCCGCCAGTCAGCATctgct                                                        >  1:8425017/1‑150 (MQ=255)
                                                                                                aCCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATctgctg                                                       >  2:8177187/1‑150 (MQ=255)
                                                                                                  cccGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAACATGCGCCACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGc                                                     <  1:6504721/148‑1 (MQ=255)
                                                                                                    ggggTTTTTTTTCGATGTCTGTTTCCACCCTCGAGTAATAATATGCTCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGa                                                   <  1:7300955/149‑1 (MQ=255)
                                                                                                     gggATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACGGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGAt                                                  <  1:4776354/150‑1 (MQ=255)
                                                                                                        aTTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCAATGtt                                               >  2:1212047/1‑150 (MQ=255)
                                                                                                            ttttCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGtt                                           >  2:2589301/1‑150 (MQ=255)
                                                                                                               tCGATGTCTGTTCCCACCCTCGAGCCAGAAAAGGCGCACCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCg                                        <  1:2577822/150‑1 (MQ=255)
                                                                                                                cGATGTCTGTTTCCACCCTCGAGTCAGAAATTGCGCACCGGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCgg                                       <  1:3202348/150‑1 (MQ=255)
                                                                                                                cGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCACCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCgg                                       <  2:4202094/150‑1 (MQ=255)
                                                                                                                    gTCTGTTGCCACCCTCGCGTCAGAAAATGCGCCCCGGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGAt                                   <  1:5190179/150‑1 (MQ=255)
                                                                                                                       tGTTTCCACCCTCGAGGCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCAc                                >  2:6577251/1‑150 (MQ=255)
                                                                                                                           tCCACCCTCGAGTCAGAAAATGCGCACCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCa                            <  2:4688773/150‑1 (MQ=255)
                                                                                                                           accACCCTCGCGTCAGAAAATGCGCAACCGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCa                            <  1:7934979/149‑1 (MQ=255)
                                                                                                                            ccACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAg                           >  2:9179872/1‑150 (MQ=255)
                                                                                                                             cACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTACTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAgc                          >  2:809464/1‑150 (MQ=255)
                                                                                                                              aCCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAgcg                         >  2:7228078/1‑150 (MQ=255)
                                                                                                                              aCCCTCGAGTCAGAAAAGGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAgcg                         <  1:2761901/150‑1 (MQ=255)
                                                                                                                                  tCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTAt                     <  2:7616569/150‑1 (MQ=255)
                                                                                                                                      gTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAAt                 >  2:6463116/1‑150 (MQ=255)
                                                                                                                                          gAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGc             >  2:2073497/1‑150 (MQ=255)
                                                                                                                                           aaaaTGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGCTTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGcc            <  1:6157921/150‑1 (MQ=255)
                                                                                                                                              aTGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAgg         <  1:8310310/150‑1 (MQ=255)
                                                                                                                                              aTGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAgg         >  1:5826410/1‑150 (MQ=255)
                                                                                                                                                gcgcAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGgcg       >  2:4870890/1‑150 (MQ=255)
                                                                                                                                                gcgcAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGgcg       <  2:3369105/150‑1 (MQ=255)
                                                                                                                                                 cgcAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGgcgc      >  2:299715/1‑150 (MQ=255)
                                                                                                                                                  gcAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCt     >  1:7209181/1‑150 (MQ=255)
                                                                                                                                                   cAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTg    <  2:5664585/150‑1 (MQ=255)
                                                                                                                                                    aaCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTgg   <  2:4260781/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATATAAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTGGG  >  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: