Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,190,776 T→G K404T (AAG→ACG)  purB ← adenylosuccinate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,190,7760TG96.6% 1742.6 / 27.5 468K404T (AAG→ACG) purBadenylosuccinate lyase
Reads supporting (aligned to +/- strand):  ref base T (8/8);  new base G (226/225);  total (235/233)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

Only 100 of 468 total aligned reads displayed.

ATCAACCATCGTGATAGCTCGACCAATATAGTTAGCCGGCGTCATCGCTTTCAGGCGGGCTTTCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCTTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAA  >  NC_000913/1190628‑1190923
                                                                                                                                                    |                                                                                                                                                   
atcaACCATCGTGATAGCTCGACCAATATAGTTAGCCGGCGTCATCGCTTTCAGGCGGGCTTTCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGtc                                                                                                                                                   <  8:772373/151‑1 (MQ=255)
    aCCATCGTGATAGCTCGACCAATATAGTTAGCCGGCGTCATCGCTTTCAGGCGGGCTTTCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGt                                                                                                                                               <  7:1435091/151‑1 (MQ=255)
     ccATCGTGATAGCTCGACCAATATAGTTAGCCGGCGTCATCGCTTTCAGGCGGGCTTTCTCTTCTTCTGGCGAGGGCAGAGCTTCGATAGACTGCTTCATGGCTTCGGCGTCAACGAGCTTACCGCGAGTGAGCTGTTTCAGCGTGTCGt                                                                                                                                               >  4:1351418/1‑150 (MQ=255)
     ccATCGTGATAGCTCGACCAATATAGTTAGCCGGCGTCATCGCTTTCAGGCGGGCTTTCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTa                                                                                                                                              <  2:154605/151‑1 (MQ=255)
      cATCGTGATAGCTCGACCAATATAGTTAGCCGGCGTCATCGCTTTCAGGCGGGCTTTCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTAc                                                                                                                                             <  5:200367/151‑1 (MQ=255)
      cATCGTGATAGCTCGACCAATATAGTTAGCCGGCGTCATCGCTTTCAGGCGGGCTTTCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTAc                                                                                                                                             >  8:630201/1‑151 (MQ=255)
         cGTGATAGCTCGACCAATATAGTTAGCCGGCGTCATCGCTTTCAGGCGGGCTTTCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGt                                                                                                                                          <  2:493/151‑1 (MQ=255)
            gATAGCTCGACCAATATAGTTAGCCGGCGTCATCGCTTTCAGGCGGGCTTTCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGtttt                                                                                                                                       <  4:1462884/151‑1 (MQ=255)
            gATAGCTCGACCAATATAGTTAGCCGGCGTCATCGCTTTCAGGCGGGCTTTCTCTTCTTCTGGCAACGCCAGACCATCCATAAACTCCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGtttt                                                                                                                                       >  6:1484796/1‑151 (MQ=255)
             aTAGCTCGACCAATATAGTTAGCCGGCGTCATCGCTTTCAGGCGGGCTTTCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGttttt                                                                                                                                      <  5:913376/151‑1 (MQ=255)
              tAGCTCGACCAATATAGTTAGCCGGCGTCATCGCTTTCAGGCGGGCTTTCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGttttt                                                                                                                                      >  1:996095/1‑150 (MQ=255)
              tAGCTCGACCAATATAGTTAGCCGGCGTCATCGCTTTCAGGCGGGCTTTCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTc                                                                                                                                     >  5:603481/1‑151 (MQ=255)
               aGCTCGACCAATATAGTTAGCCGGCGTCATCGCTTTCTGGCGGGCTTTCTCTTCTTCTGGCACCGCAACTACATCGATAAACTGCTTCATGACGTCGGCGTCAACGCGGTTACCGCGAGTCAGGTCTTTCAGCGTCTCGTCCGGTTTTTCg                                                                                                                                    >  6:1645659/1‑151 (MQ=255)
               aGCTCGACCAATATAGTTAGCCGGCGTCATCGCTTTCAGGCGGGCTTTCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTc                                                                                                                                     <  7:1593892/150‑1 (MQ=255)
                 cTCGACCAATATAGTTAGCCGGCGTCATCGCTTTCAGGCGGGCTTTCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGAt                                                                                                                                  >  4:535121/1‑151 (MQ=255)
                 cTCGACCAATATAGTTAGCCGGCGTCATCGCTTTCAGGCGGGCTTTCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGAt                                                                                                                                  >  2:1059951/1‑151 (MQ=255)
                  tCGACCAATATAGTTAGCCGGCGTCATCGCTTTCAGGCGGGCTTTCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATg                                                                                                                                 >  4:1655169/1‑151 (MQ=255)
                  gcGACCAATATAGTTAGCCGGCGGCATCGCTGTCAGGCGGGCTTTCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATg                                                                                                                                 <  4:974003/150‑1 (MQ=255)
                      ccAATATAGTTAGCCGGCGTCATCGCTTTCAGGCGGGCTTTCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCAt                                                                                                                             <  1:1003947/151‑1 (MQ=255)
                        aaTATAGTTAGCCGGCGTCATCGCTTTCAGGCGGGCTTTCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCAt                                                                                                                             >  2:295377/1‑149 (MQ=255)
                        aaTATAGTTAGCCGGCGTCATCGCTTTCAGGCGGGCTTTCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCAt                                                                                                                             >  3:1480881/1‑149 (MQ=255)
                         atatAGTTAGCCGGCGTCATCGCTTTCAGGCGGGCTTTCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGc                                                                                                                          >  7:1148405/1‑151 (MQ=255)
                             aGTTAGCCGGCGTCATCGCTTTCAGGCGGGCTTTCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGAc                                                                                                                       <  1:478599/150‑1 (MQ=255)
                             aGTTAGCCGGCGTCATCGCTTTCAGGCGGGCTTTCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACg                                                                                                                      >  4:342295/1‑151 (MQ=255)
                                 aGCCGGCGGCATCGCTTTCAGGCGGGCTTTCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATa                                                                                                                  <  2:1179654/151‑1 (MQ=255)
                                  gCCGGCGTCATCGCTTTCAGGCGGGCTTTCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATaa                                                                                                                 <  2:1113967/151‑1 (MQ=255)
                                   ccGGCGTCATCGCTTTCAGGCGGGCTTTCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATaaa                                                                                                                <  3:1136752/151‑2 (MQ=255)
                                    cGGCGTCATCGCTTTCAGGCGGGCTTTCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACt                                                                                                               >  2:1644918/1‑151 (MQ=255)
                                    cGGCGTCATCGCTTTCAGGCGGGCTTTCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACt                                                                                                               >  1:354059/1‑151 (MQ=255)
                                     ggCGTCATCGCTTTCAGGCGGGCTTTCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTg                                                                                                              >  4:1007727/1‑151 (MQ=255)
                                      gCGTCATCGCTTTCAGGCGGGCTTTCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGt                                                                                                             <  5:791522/151‑1 (MQ=255)
                                       cGTCATCGCTTTCAGGCGGGCTTTCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTc                                                                                                            <  4:1140573/151‑1 (MQ=255)
                                          cATCGCTTTCAGGCGGGCTTTCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTgg                                                                                                         >  8:1034411/1‑151 (MQ=255)
                                                 ttCAGGCGGGCTTTCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCGGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGt                                                                                                   <  8:441812/150‑1 (MQ=255)
                                                 ttCAGGCGGGCTTTCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCTTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGt                                                                                                   <  2:98351/150‑1 (MQ=255)
                                                   cAGGCGGGCTTTCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTca                                                                                                <  5:1261799/151‑1 (MQ=255)
                                                   cAGGCGGGCTTTCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTc                                                                                                 <  1:528786/150‑1 (MQ=255)
                                                    aGGCGGGCTTTCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTcag                                                                                               <  5:567856/151‑1 (MQ=255)
                                                    aGGCGGGCTTTCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTc                                                                                                 >  3:530272/1‑149 (MQ=255)
                                                    aGGCGGGCTTTCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTc                                                                                                 >  5:1130491/1‑149 (MQ=255)
                                                       cGGGCTTTCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTcagcca                                                                                            <  7:704744/151‑1 (MQ=255)
                                                        gggCTTTCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTcagcc                                                                                             >  8:1188735/1‑149 (MQ=255)
                                                        gggCTTTCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTcagcc                                                                                             >  6:133800/1‑149 (MQ=255)
                                                          gCTTTCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTcagccagc                                                                                          >  3:709859/1‑150 (MQ=255)
                                                           cTTTCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTcagccagccc                                                                                        >  4:254532/1‑149 (MQ=255)
                                                            tttCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACt                                                                                       >  3:613905/1‑151 (MQ=255)
                                                             ttCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCAc                                                                                        <  5:410437/149‑1 (MQ=255)
                                                              tctcTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCCGCCAGCACTTc                                                                                     >  7:1198798/1‑151 (MQ=255)
                                                              tctcTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACtt                                                                                      <  7:428095/150‑1 (MQ=255)
                                                               ctctTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTcc                                                                                    >  5:1647799/1‑151 (MQ=255)
                                                                  ttcttctGGCAACGCAAGACCATCGATAAACAGCTTCATGACTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGACATAGCGACGAATAACTGAGTCGATTGGTTCAGCCAGCACTTCCCat                                                                                 >  1:1574503/1‑150 (MQ=255)
                                                                   tcttctGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAg                                                                                 <  3:1347638/150‑1 (MQ=255)
                                                                   tcttctGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAg                                                                                 <  2:597595/150‑1 (MQ=255)
                                                                    cttctGGCAACGCCAGACCATCGATAAAATGCTTCATGACTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCCTCTCGTACGGTTTTTCGATGACATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCTCAGtt                                                                               >  8:1691954/1‑151 (MQ=255)
                                                                      tctGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTtg                                                                              <  3:391209/150‑1 (MQ=255)
                                                                       ctGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTtgtg                                                                            >  8:1270601/1‑151 (MQ=255)
                                                                        tGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCTTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTtgtg                                                                            >  8:1365846/1‑150 (MQ=255)
                                                                         ggCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTtgt                                                                             <  8:148160/148‑1 (MQ=255)
                                                                             aCGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATcc                                                                        >  1:698041/1‑149 (MQ=255)
                                                                             aCGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAg                                                                      <  7:1084615/151‑1 (MQ=255)
                                                                              cGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGt                                                                     >  8:366071/1‑151 (MQ=255)
                                                                               gCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGtt                                                                    <  7:180064/151‑1 (MQ=255)
                                                                               gCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGATCAGCCAGCACTTCCCAGTTGTGATACAGtt                                                                    >  2:938555/1‑151 (MQ=255)
                                                                                  aGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCAt                                                                 >  2:1478168/1‑151 (MQ=255)
                                                                                  aGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCAt                                                                 >  6:40297/1‑151 (MQ=255)
                                                                                   gACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATc                                                                <  3:854276/151‑1 (MQ=255)
                                                                                      cATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCcag                                                             <  6:961767/151‑1 (MQ=255)
                                                                                       aTCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCcagc                                                            >  4:1499465/1‑151 (MQ=255)
                                                                                         cGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCcagcag                                                          <  5:1549324/151‑1 (MQ=255)
                                                                                         cGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCcagca                                                           <  7:614636/150‑1 (MQ=255)
                                                                                           aTAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGAt                                                        <  3:732157/151‑1 (MQ=255)
                                                                                             aaaCTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATgg                                                      >  4:394999/1‑151 (MQ=255)
                                                                                             aaaCTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATg                                                       <  6:734212/150‑1 (MQ=255)
                                                                                                  gCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACg                                                 <  6:1050375/151‑1 (MQ=255)
                                                                                                     tCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCTTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGACTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGt                                               <  4:83712/150‑1 (MQ=255)
                                                                                                      cATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTc                                             >  4:635416/1‑151 (MQ=255)
                                                                                                          ccTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTcccgt                                         >  1:966514/1‑147 (MQ=255)
                                                                                                             tCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTcc                                       >  8:1604650/1‑150 (MQ=255)
                                                                                                                gCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGtt                                   <  5:1405452/151‑1 (MQ=255)
                                                                                                                gCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGtt                                   <  7:36544/151‑1 (MQ=255)
                                                                                                                     aaCGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTc                              <  7:947165/151‑1 (MQ=255)
                                                                                                                      aCGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGATGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCa                             <  4:1220926/151‑1 (MQ=255)
                                                                                                                      aCGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCa                             <  1:784341/151‑1 (MQ=255)
                                                                                                                      aCGCGCTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCa                             <  2:824904/151‑1 (MQ=255)
                                                                                                                          gcTTACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGcc                         >  1:1214862/1‑151 (MQ=255)
                                                                                                                            ttACCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCtt                       >  8:23008/1‑151 (MQ=255)
                                                                                                                              aCCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTc                     >  7:249351/1‑151 (MQ=255)
                                                                                                                              aCCGCGAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTc                     >  3:280726/1‑151 (MQ=255)
                                                                                                                               ccGCGTGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTAGCTGACGCCTTTCa                    <  2:1304432/151‑1 (MQ=255)
                                                                                                                                cgcgAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAg                   <  6:670757/151‑1 (MQ=255)
                                                                                                                                  cgAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAggg                 <  1:1650363/151‑1 (MQ=255)
                                                                                                                                  cgAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAggg                 <  2:1258009/151‑1 (MQ=255)
                                                                                                                                   gAGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGt                <  8:939553/151‑1 (MQ=255)
                                                                                                                                    aGTCAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTg               >  7:821185/1‑151 (MQ=255)
                                                                                                                                       cAGCTCTTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGAt            <  1:1394802/151‑1 (MQ=255)
                                                                                                                                            ctTTCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGa         <  1:874169/149‑1 (MQ=255)
                                                                                                                                              ttCAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACTCCTTTCAGGGTGGATTGata       <  6:245735/149‑1 (MQ=255)
                                                                                                                                                cAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGc    <  6:1149314/150‑1 (MQ=255)
                                                                                                                                                cAGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACGTCCCATTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGc    >  5:679060/1‑150 (MQ=255)
                                                                                                                                                 aGCGTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCa   >  3:1267755/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATCAACCATCGTGATAGCTCGACCAATATAGTTAGCCGGCGTCATCGCTTTCAGGCGGGCTTTCTCTTCTTCTGGCAACGCCAGACCATCGATAAACTGCTTCATGCCTTCGGCGTCAACGCGCTTACCGCGAGTCAGCTCTTTCAGCTTCTCGTACGGTTTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAA  >  NC_000913/1190628‑1190923

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