Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,828,980 A→C 8.8% 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 8.8% 1074.0 / 15.3 435Q12H (CAA→CACxanPxanthine permease
Reads supporting (aligned to +/- strand):  ref base A (229/156);  new base C (12/38);  total (241/194)
Fisher's exact test for biased strand distribution p-value = 2.97e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 443 total aligned reads displayed.

ATATAAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTGGG  >  NC_000913/3828831‑3829129
                                                                                                                                                     |                                                                                                                                                     
 tataAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAAc                                                                                                                                                      >  1:6853286/1‑150 (MQ=255)
   taAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCg                                                                                                                                                    >  1:4527965/1‑150 (MQ=255)
    aaaaaaGTGCATTTATCTCATTTAAATCTTTAGATGAGGATGTGGCAAACGTTTGCTGTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGAGTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAACATGCGCACCCgg                                                                                                                                                   >  2:6745510/1‑150 (MQ=255)
     aaaaaGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGt                                                                                                                                                  >  1:5981360/1‑150 (MQ=255)
     aaaaaGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGt                                                                                                                                                  >  2:532951/1‑150 (MQ=255)
      aaaaGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCCGtt                                                                                                                                                 >  1:8216802/1‑150 (MQ=255)
           tGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCCCAACCGGTTGCGCa                                                                                                                                            >  2:1514828/1‑150 (MQ=255)
             cATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCACTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCatt                                                                                                                                          <  2:3894914/150‑3 (MQ=255)
              aTTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGAc                                                                                                                                         >  1:1015179/1‑150 (MQ=255)
                    ctcATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGTTTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCaaaa                                                                                                                                   <  2:4207228/150‑1 (MQ=255)
                    ctcATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCaaaa                                                                                                                                   <  2:240505/150‑1 (MQ=255)
                     tcctttAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAACATGCGCAACCGGTTGCGCAGACTCAAAAc                                                                                                                                  >  2:5423281/4‑150 (MQ=255)
                     tcATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACACaaaaa                                                                                                                                  >  1:2312268/1‑149 (MQ=255)
                        tttAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTTCGCAGACTCAAAACAGc                                                                                                                               >  2:7098256/1‑150 (MQ=255)
                        tttAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGc                                                                                                                               >  1:6207812/1‑150 (MQ=255)
                           aaaTATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGaa                                                                                                                            <  2:1952318/150‑1 (MQ=255)
                             atatTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCCACCGGTTTCGCAGACTCAAAAAAGCGAAca                                                                                                                          >  1:2965753/1‑149 (MQ=255)
                             atatTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCACCCGGTTGCGCAGACTCAAAACAGCGAACt                                                                                                                          >  2:4997662/1‑150 (MQ=255)
                                 ttAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACAGAtt                                                                                                                      >  1:1283014/1‑150 (MQ=255)
                                     aTGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTAcc                                                                                                                  <  2:3374122/150‑1 (MQ=255)
                                      tGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAAAGCCCGCTTTGCTACCACGGGATTGTTTTCGATGACTGTTTCCACCCTCGAGTCAGAAAATGAGCAAACCGGTGCGCAGACTCAAAAAACCGAAATGATTTAccc                                                                                                                 >  1:5573289/1‑149 (MQ=255)
                                       gAGGATGTGGCACACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGt                                                                                                                <  2:6230012/150‑1 (MQ=255)
                                         ggATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCt                                                                                                              >  2:7461568/1‑150 (MQ=255)
                                             gtgGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGaa                                                                                                          <  2:4952518/150‑1 (MQ=255)
                                              tgGAAAATGTTTGATTTTTTCTTTGGGTCAGATAAAATGCCCGCTTTGCTACCAGGGTATTTTTTTCTATGTCTGTTTCCACCCTCGAGTCATAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAg                                                                                                         <  1:6092468/150‑1 (MQ=255)
                                               ggCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGaact                                                                                                        >  2:6538563/1‑148 (MQ=255)
                                               ggCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGaaca                                                                                                        >  2:1123187/1‑148 (MQ=255)
                                                 aaaagttttGATTTTCTCTTTGGTTCCGATAAAATGCCCGTTTTGCTACCCGGTTTTTTTTTTCTTTTTGTGTTTCCACCCTCGAGACAGAAAAGGCGCACCGGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATc                                                                                                      <  1:1161268/144‑1 (MQ=255)
                                                    aCGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTc                                                                                                   <  2:106461/150‑1 (MQ=255)
                                                        ttGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTccgcc                                                                                               >  2:6908823/1‑150 (MQ=255)
                                                         tGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCCGTAGCGCCCACACAAAACACCCAACAGATTTAACCTCTTGAACATCGTaccccg                                                                                              >  1:3020991/1‑150 (MQ=255)
                                                           cTTTTCTCTTTCGGTCAGCTAATTGGCCCGCTTTGCTACCTGGGGTTTGTTTTCTTTGTGTGTTTCCCCCCTCTAGTCATAAAATGCGCACCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCt                                                                                            <  1:5423163/150‑1 (MQ=255)
                                                               tctcTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCt                                                                                        <  2:6981170/150‑1 (MQ=255)
                                                               tctcTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCt                                                                                        >  2:8118491/1‑150 (MQ=255)
                                                                 tctTTCGGTCAGCTAAAATGCCCGCTTAGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCCCCCGGTTGCGCAGACACAAAACAGCGAACAGATTTACACTCTAGAAGATAGACCGCCGCATACTCa                                                                                      >  2:3993575/1‑150 (MQ=255)
                                                                  ctTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCaa                                                                                     >  2:645239/1‑150 (MQ=255)
                                                                   tttCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGGCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCaaa                                                                                    >  2:3223482/1‑150 (MQ=255)
                                                                        gTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGTTCGTCCGCCGCTTCCTCAAACCCTg                                                                               >  1:2598373/1‑150 (MQ=255)
                                                                         tCAGATAAAATGCCCTCTTTGCTACCTGGTTATTTTTTTCGATGTCTGTTCCCACCCTCTAGTAAGAAAATGCGCACCGGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGt                                                                              <  1:4997902/150‑1 (MQ=255)
                                                                           aGCTAAAATGCCGGCTTTGCTACCAGGGGTTTTTTTTCGATTTCTTTTTCCACCCTCGAGTCACAATTTGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCATCAAACCCTGttt                                                                            <  1:903654/150‑1 (MQ=255)
                                                                           aGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCACCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGttt                                                                            >  1:316121/1‑150 (MQ=255)
                                                                           aGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATAGTCCGCCGCTTCCTCAAACCCGGttt                                                                            >  1:6527453/1‑150 (MQ=255)
                                                                            gCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAAACGGTTGCGCAGAATCAAAAAAGCGAACTGATTTACCGTCTTGAAGAACGTCCGCCGCGTCCACAAACCCTGTTTg                                                                           >  1:637589/1‑150 (MQ=255)
                                                                              tAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTgcc                                                                         >  2:1564983/1‑150 (MQ=255)
                                                                                aaaTGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGGTGCGCAGACTCAAAACAGCGAACTGATTTAACGTCTTGTAGATCGTCCGCCGCGTCCACAAACCCAGTTTgccgc                                                                       >  1:726826/1‑150 (MQ=255)
                                                                                  tggCCCGCTTTGCTACCCCGGGTTTGTTTTCGATGTCTGTTTCCCCCCTCGTGTCCGAAAATGCGCACCGGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCt                                                                     <  1:5368484/148‑1 (MQ=255)
                                                                                  aTGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCt                                                                     <  2:4770374/150‑1 (MQ=255)
                                                                                       cGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCCGTTGCGCAGACACAAAACAGCGAACAGATTTACCGTCTTTAAGATAGTCCGCCCCATACACAAACCCAGTTTGCCGCCAGTCAg                                                                >  1:3606104/1‑150 (MQ=255)
                                                                                        gCTTTGCTACCACGGGTTTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGc                                                               <  2:4444202/150‑1 (MQ=255)
                                                                                        gCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGc                                                               <  2:3275929/150‑1 (MQ=255)
                                                                                           ttGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGGCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATc                                                            <  2:2472624/150‑1 (MQ=255)
                                                                                            tGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATct                                                           >  2:8245665/1‑150 (MQ=255)
                                                                                                aCCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAAGGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATctgctg                                                       <  1:6993312/150‑1 (MQ=255)
                                                                                                 ccACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTgg                                                      <  1:532951/150‑1 (MQ=255)
                                                                                                  cACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCACCCGGTTGCGCAGACTCAAAACAGCGAACTGAATTACCGTCTTGAAGATCGTCCGCCGCGTCCGCAAACCCTGGTTGCCGCCTGTCAGCGTCTGCTGGc                                                     >  2:6134553/1‑150 (MQ=255)
                                                                                                  cACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGc                                                     <  2:7375925/150‑1 (MQ=255)
                                                                                                      ggtttGTTTTCGTTGTCTGTTTCCCCCCTCGAGTCAGAAACGGCGCACCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATg                                                 <  1:4724743/147‑1 (MQ=255)
                                                                                                        aTTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCCCCCGGGTGCGCAGACACAAAAAAGCGAACAGATTTAACGGCGTGAAGATCGTCCGCCGCGTCCACAAAACCGGGTTGCCGCCTGTCAGCAAATGCTGGCGCTGtt                                               >  1:1436759/1‑150 (MQ=255)
                                                                                                             tttCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTg                                          <  2:1852139/150‑1 (MQ=255)
                                                                                                              ttCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGc                                         >  2:2264902/1‑150 (MQ=255)
                                                                                                                 gATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGt                                      >  2:5744605/1‑150 (MQ=255)
                                                                                                                 gATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAAAAGCGAACTGATTTACAGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGt                                      >  2:4297076/1‑150 (MQ=255)
                                                                                                                  aTGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCCACCGGTTTCGCACACTAAAAACAGCGAAATGTTTTAACGTCTTTAAGATCGTCCGCCGCTTACTAAAACCCAGTTTTGCGCCGGTCAGCATATGCTGGCGATGTTGGTTGCGGTg                                     >  1:5752845/1‑150 (MQ=255)
                                                                                                                    gTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGAt                                   <  1:32729/150‑1 (MQ=255)
                                                                                                                     tCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATc                                  >  2:840092/1‑150 (MQ=255)
                                                                                                                     tCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCGTCCTCAAACCCTGTTTGCCGCCGGTCAGCATCTGCTGGCGATGTTCGTTGCGGGGATc                                  >  2:7774828/1‑150 (MQ=255)
                                                                                                                      cTGTTTCCACCCTCGAGTCAGAAAATGCGCACCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCa                                 <  2:5617547/150‑1 (MQ=255)
                                                                                                                      cTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCa                                 >  2:4551331/1‑150 (MQ=255)
                                                                                                                       tGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCAc                                >  2:2874246/1‑150 (MQ=255)
                                                                                                                       tGTTGCCACCCTCGAGTAATAAAATGCGAAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCAc                                <  1:8359484/150‑1 (MQ=255)
                                                                                                                         tttCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGc                              <  2:2965753/150‑1 (MQ=255)
                                                                                                                          ttCCACCCTCGAGTCAGAAAATGCGCACCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGcc                             <  1:8099083/150‑1 (MQ=255)
                                                                                                                           tCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCGTGAAGATCGTCCGCCGCGTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCa                            >  2:920463/1‑150 (MQ=255)
                                                                                                                            ccACCCTCGAGTCATAAAATGCGCACCGGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAg                           <  1:4558182/150‑1 (MQ=255)
                                                                                                                              aCCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAgcg                         <  2:5709913/150‑1 (MQ=255)
                                                                                                                               cccTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAgcgc                        <  2:2594728/150‑1 (MQ=255)
                                                                                                                                ccTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCt                       <  1:2120636/150‑1 (MQ=255)
                                                                                                                                 cTCGAGTCAGAAATTGCGCACCGGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTa                      <  1:3998043/150‑1 (MQ=255)
                                                                                                                                 cTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTa                      >  2:283354/1‑150 (MQ=255)
                                                                                                                                  tCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTAt                     >  2:7602356/1‑150 (MQ=255)
                                                                                                                                    gAGTCAGAAAATGCGCCACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTAACGTCTTGAAGATCGTCCGCCGCTTCCACAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCCGCGCTATTa                   >  1:5254823/1‑150 (MQ=255)
                                                                                                                                     aGTCAGAAAATGCGCACCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTaa                  <  1:6769051/150‑1 (MQ=255)
                                                                                                                                     aGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTaa                  <  1:3845701/150‑1 (MQ=255)
                                                                                                                                      gTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAAt                 <  1:1795237/150‑1 (MQ=255)
                                                                                                                                       tCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATc                >  2:3441047/1‑150 (MQ=255)
                                                                                                                                        cAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCt               >  2:2113764/1‑150 (MQ=255)
                                                                                                                                        cAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCt               <  1:4962342/150‑1 (MQ=255)
                                                                                                                                          gAAAATGCGCAACCGGTTGCGCAGACTAAAAACAGCGAACTGATTTAACGTCTTGAAGATCGTCCGCCGCATCCTCAAAACCTGTTTTCCGCCAGTCAGCAAATGCTGGCGATGTTCGTTGCGGTGATAACGCCAGCGCTATTAATCTGc             >  2:3628857/1‑150 (MQ=255)
                                                                                                                                           aaatTGCTCAACGGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCTTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGca            <  2:3621476/150‑2 (MQ=255)
                                                                                                                                            aaaTGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCa           <  2:2035913/150‑1 (MQ=255)
                                                                                                                                            aaaTGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCa           <  2:5930775/150‑1 (MQ=255)
                                                                                                                                             tctGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAg          >  2:1727596/3‑150 (MQ=255)
                                                                                                                                               tGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGgc        >  2:3550422/1‑150 (MQ=255)
                                                                                                                                                 cgcAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGgcgc      >  2:7276483/1‑150 (MQ=255)
                                                                                                                                                 cgcAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGgcgc      >  2:4122843/1‑150 (MQ=255)
                                                                                                                                                 cgcAACCGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGgcgc      <  1:272687/150‑1 (MQ=255)
                                                                                                                                                   cAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTg    >  1:7641155/1‑150 (MQ=255)
                                                                                                                                                   cAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTg    <  2:1314453/150‑1 (MQ=255)
                                                                                                                                                   cAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTg    <  1:4274928/150‑1 (MQ=255)
                                                                                                                                                    aaCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTgg   <  2:5380962/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: