Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,263,502 T→C 5.7% A53A (GCA→GCG tdcD ← propionate kinase/acetate kinase C, anaerobic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,263,5020TC 5.7% 239.6 / 4.9 89A53A (GCA→GCGtdcDpropionate kinase/acetate kinase C, anaerobic
Reads supporting (aligned to +/- strand):  ref base T (27/53);  new base C (3/2);  total (34/55)
Fisher's exact test for biased strand distribution p-value = 3.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.86e-01

ATCGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAA  >  NC_000913/3263355‑3263651
                                                                                                                                                   |                                                                                                                                                      
aTCGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGttc                                                                                                                                                      >  1:1132201/1‑148 (MQ=255)
aTCGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCAAGCGGTgg                                                                                                                                                      >  1:1738373/1‑149 (MQ=255)
 tCGGTAATAATGGCGGACTCGGTAAAAATCCTGACTCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTATGGTGAGCCAGCGGTGCt                                                                                                                                                     <  2:113317/150‑1 (MQ=255)
 tCGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCt                                                                                                                                                     <  1:1028030/150‑1 (MQ=255)
 tCGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTACCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCt                                                                                                                                                     <  1:250517/150‑1 (MQ=255)
  cGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCACTTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTg                                                                                                                                                    <  1:1365005/150‑1 (MQ=255)
   ggTAATAATGGCGGACTCGGTAATAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGATCAAATGCAATTGCTTTCATTGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTgg                                                                                                                                                   <  1:1449903/150‑1 (MQ=255)
    gTAATTATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTCTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCACCGCTGGc                                                                                                                                                  <  2:183710/150‑1 (MQ=255)
    gTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGc                                                                                                                                                  <  1:487896/150‑1 (MQ=255)
    gTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGc                                                                                                                                                  <  1:1087785/150‑1 (MQ=255)
    gTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGc                                                                                                                                                  <  2:47362/150‑1 (MQ=255)
     taataaTGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGct                                                                                                                                                 >  1:1110674/1‑150 (MQ=255)
      aataatGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCACCGCTGGctc                                                                                                                                                <  1:1416102/150‑1 (MQ=255)
         aatGGCGGACTCGTTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTcc                                                                                                                                            <  2:784417/150‑1 (MQ=255)
               ggACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCAtt                                                                                                                                      <  1:1582835/150‑1 (MQ=255)
                gACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCAttt                                                                                                                                     <  2:604106/150‑1 (MQ=255)
                 aCTCGGTAAAAATAGTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCAACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTa                                                                                                                                    <  2:616293/150‑1 (MQ=255)
                  cTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGTCAGCGGTGCTGGCTCTCCCCCATTTAc                                                                                                                                   <  1:1375684/150‑1 (MQ=255)
                   tCGGTAAATATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACg                                                                                                                                  <  1:290742/150‑1 (MQ=255)
                    cGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAACTTCCGTTTTTCCACTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACgg                                                                                                                                 <  1:484958/150‑1 (MQ=255)
                    cGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACgg                                                                                                                                 <  1:35095/150‑1 (MQ=255)
                     ggTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGGGCCAGCGGGGCGGGCTCTCCCCCCTTTACggt                                                                                                                                >  1:1806329/1‑149 (MQ=255)
                       tAAAAATACGGCAGCCGTGATCGATGCGGTGGCAATTTAAGGCAACTCTGTAATATAAATTCCGTTTTTACAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGGTGTGGTGCGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATa                                                                                                                              <  2:1629112/150‑1 (MQ=255)
                         aaaaTACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCCCCGGGGCTGGGGCTACCCCCTTTTAGGATaat                                                                                                                            >  2:1579551/1‑149 (MQ=255)
                         aaaaTACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAg                                                                                                                            <  2:526281/150‑1 (MQ=255)
                             tACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATg                                                                                                                        <  2:1296977/150‑1 (MQ=255)
                               cTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCa                                                                                                                      >  1:261464/1‑150 (MQ=255)
                                tGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCAt                                                                                                                     <  2:382645/150‑1 (MQ=255)
                                tGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCAt                                                                                                                     >  2:970837/1‑150 (MQ=255)
                                 gccgccGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCCGCGGTGGTGGGTCGCCCCCCATTTCGGATAAGAATGCaat                                                                                                                    >  2:1856252/1‑148 (MQ=255)
                                       gTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGa                                                                                                              <  1:1305905/150‑1 (MQ=255)
                                        ggAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAg                                                                                                             <  2:1601279/149‑1 (MQ=255)
                                            cGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTaa                                                                                                         <  1:1597817/150‑1 (MQ=255)
                                             gATGCGGTGGCCAATTAAGGCCCCACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAAt                                                                                                        >  1:239586/1‑150 (MQ=255)
                                              aTGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAACTGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATa                                                                                                       <  1:903799/150‑1 (MQ=255)
                                                gCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCTTTTTTCCAGTTCAAATGCAATTGCCTTAAATGCACCTACGTAGCTGTGGTGAGCCAGCGGTTCTGGCTCTCCCCCATTTACGCATAAAAATGCATTTTACGAGTTAATaac                                                                                                     >  2:960559/1‑148 (MQ=255)
                                                gCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATAcc                                                                                                     <  2:1296232/150‑1 (MQ=255)
                                                gCGGTGGACAATTAAGGCCACACTGTAATTTAATTTCCGTTTTTCCAGTTCTAATGCAATTGCCTTCAATGCTCCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATAcc                                                                                                     <  1:1439720/150‑1 (MQ=255)
                                                   gTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCGCTCCCCCCTTTAAGGTGTAGAATGCATTTTACGAGTTAATACCGTc                                                                                                  >  2:886443/1‑150 (MQ=255)
                                                      gCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGc                                                                                               <  2:1874051/150‑1 (MQ=255)
                                                      gCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTGCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGc                                                                                               <  2:1657248/150‑1 (MQ=255)
                                                       ccAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGGGCCAGCGGTGCTGGCTCAACCACCATTAAGGATAAAAAATCAATTTAACAGTTAATAACCTCGGCa                                                                                              >  1:1146909/1‑150 (MQ=255)
                                                       ccAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCa                                                                                              <  1:49064/150‑1 (MQ=255)
                                                       ccAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCa                                                                                              >  1:255935/1‑150 (MQ=255)
                                                       ccAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAACCGTGCTGGGTAGCCCCCAGTTAGGGATAAAAATGCATTTTACGCGTTAATACACTCGGcc                                                                                              >  1:721132/1‑149 (MQ=255)
                                                        cAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGGGGTCCCCCAGTTTAGGAGAAGAATGCATTTTGCGAGTTTATAACGTGGGCaa                                                                                             >  2:1440727/1‑150 (MQ=255)
                                                        cAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGGGGGGGGTGTCCCCCCATTTAGGTTAAAAATGGATTTTGAGAGTTAATAAAGTCGGGaa                                                                                             >  1:1613170/1‑150 (MQ=255)
                                                             aaGGCCAAACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCATTCAATGCACCTTCGTAGCTGTGGTGATCCAGCGGTGCTGGCTCTCCCCCAATTACGGAATAGAATTCAATTTCCGAGTTAATACCGTCCGCAATAcca                                                                                        >  2:1025436/1‑149 (MQ=255)
                                                              aGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGTTAAGAATGGATTTTACGAGTTAATAACGTCGGCAATaacag                                                                                       >  2:1861080/1‑146 (MQ=255)
                                                               ggCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTACGAGTTAATACCGTCGGCAATACCTGa                                                                                      >  2:1020407/1‑150 (MQ=255)
                                                               ggCCACACTGTCATTTAAATTCCGTTTTTACAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTCGCTCTCCCCCATTTACGGATATGAATGCATTTTCCGAGTTAATACCGTCGGCAATGCCTgc                                                                                      >  2:222022/1‑149 (MQ=255)
                                                                gCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGAc                                                                                     <  1:438780/150‑1 (MQ=255)
                                                                   acacTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCACCGGTGCTGGCTCTCCCCCAATTACGGATAAAAATGCATTTTACGAGTTAATACCCTCGGCAATACCAGaaaat                                                                                  >  1:1836274/1‑146 (MQ=255)
                                                                   acaaTGTCATTTAAATTCCGTTTTTCCAGTTCTAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTACCCCATTTACGGATATGAATGCTTTTTCCGAGTTAATACCGTCGGCAATAGCTGACAtt                                                                                  <  1:960559/150‑1 (MQ=255)
                                                                    cacTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGGGCCAGCGGGGGTGGCTCTGCCCCCATTTAGGTTTAGAAATCATTTTTCGAGTTAAAAAAGGCGGCAATAACAGACAttt                                                                                 >  2:1524854/1‑149 (MQ=255)
                                                                      cTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAAt                                                                               <  2:1337300/150‑1 (MQ=255)
                                                                      cTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTACGAGTTAATACCGTCGGCAATACCTGACATTAAt                                                                               >  1:966551/1‑150 (MQ=255)
                                                                      atGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAAt                                                                               <  1:1045806/149‑1 (MQ=255)
                                                                          cATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACtt                                                                           >  1:1342310/1‑150 (MQ=255)
                                                                           aTTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTAGACCCCATTTACGGTTAAGAATGCATTTTACGAGTTAATACCGTCGGCAATACCTGACATTAATACTTc                                                                          >  2:1730299/1‑150 (MQ=255)
                                                                             ttAAATTCCGTTTTTCCATTTCAAATCGAATTGCCTTCAATTCACATTCGTAGCTGTGGTTAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTcac                                                                        <  1:1124129/150‑1 (MQ=255)
                                                                             ttAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTcac                                                                        >  1:1051784/1‑150 (MQ=255)
                                                                                      gTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTgg                                                               <  2:1497783/150‑1 (MQ=255)
                                                                                           tCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCg                                                          <  2:647335/150‑1 (MQ=255)
                                                                                           tCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGGGGGGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCg                                                          <  2:1015798/150‑1 (MQ=255)
                                                                                            caaGTTCAAATGCAATTGCCTTCAATGCTCCGTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGa                                                         <  2:1199429/148‑1 (MQ=255)
                                                                                             cAGTACAACTGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAg                                                        <  2:1760660/150‑1 (MQ=255)
                                                                                              aGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGc                                                       <  1:1388092/150‑1 (MQ=255)
                                                                                                       gCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAAc                                              <  1:1519167/150‑1 (MQ=255)
                                                                                                        cAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCCGCGGTGCTGGCTCTCCCCCCTTTAAGGATAAGAATTCATTTTCAGAGTTAATACCGTCGGCAATTACTGACATTTATACTTCACAGGAGCGGGCATCGAGCCCGGaaaaaa                                             >  1:351915/1‑148 (MQ=255)
                                                                                                          aTTGCCTTCAATGAACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTa                                           <  2:1869616/150‑1 (MQ=255)
                                                                                                            tGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAAt                                         <  2:1065263/150‑1 (MQ=255)
                                                                                                             gCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATc                                        <  1:1452759/150‑1 (MQ=255)
                                                                                                             gCCTTCAATGCACCTTCGTAGCTGTGGGGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATc                                        <  2:1284689/150‑1 (MQ=255)
                                                                                                                  cAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTACGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGa                                   >  1:4433/1‑150 (MQ=255)
                                                                                                                   aaTGCACCTTCGTAGCTGTGGTGAGCCATCGGTGCTGTGTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGGTAATAACGTCGGCAATACCTGTCATTAATACTTCACAGTCGCTGGCATCTAGCACGGAAAACTTAATCGAAGAc                                  <  2:1291286/150‑1 (MQ=255)
                                                                                                                      gCACCTTCGTAGCTGTGGTGGGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGaa                               >  2:1338131/1‑150 (MQ=255)
                                                                                                                        aCCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAAcc                             <  2:239586/150‑1 (MQ=255)
                                                                                                                        aCCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAAcc                             <  1:1080204/150‑1 (MQ=255)
                                                                                                                          cTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACcac                           <  2:261464/150‑1 (MQ=255)
                                                                                                                             cGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGt                        <  1:1793703/150‑1 (MQ=255)
                                                                                                                               tAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTta                      <  2:878946/150‑1 (MQ=255)
                                                                                                                                aGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTtaa                     <  1:1524854/150‑1 (MQ=255)
                                                                                                                                     gggagAGCCAGCACCGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAAcc                >  1:823336/5‑150 (MQ=255)
                                                                                                                                      ggagAGCCAGCACCGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCAAAGTTAATAACCa               >  2:1867647/4‑150 (MQ=255)
                                                                                                                                       gtgAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGCCGAACCACAGTTAATAACCaa              >  1:848399/1‑150 (MQ=255)
                                                                                                                                       gagAGCCAGCACCGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAACACTTCACAGTCGCTGGCATCGAGCACGGAAAACATAATCGAAGACGAACCACAGTtaataaacca              >  1:331575/3‑146 (MQ=255)
                                                                                                                                           gCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCTGCAATACCTGACCTTAATACTTCACAGTCGCTGGCATCGATCACGGAAAACTTAATAGAAGACGATCCACAGTTAATACCCAAcaca          >  2:1121790/1‑150 (MQ=255)
                                                                                                                                           gCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGACAACTTAATCGAAGACGACCCACACTTAATAACCAAaaca          >  2:165897/1‑150 (MQ=255)
                                                                                                                                               gCGGTGCTGGCTCTCCCACATTTACGGATAAGAATGCATTTTCCGAGTTAATACAGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCg      <  2:1182319/150‑1 (MQ=255)
                                                                                                                                                   tGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGaaa  >  2:981552/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
ATCGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAA  >  NC_000913/3263355‑3263651

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