Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,263,501 G→C 5.2% P54A (CCG→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,5010GC 5.2% 279.9 / 5.2 90P54A (CCG→GCG) tdcDpropionate kinase/acetate kinase C, anaerobic
Reads supporting (aligned to +/- strand):  ref base G (32/53);  new base C (3/2);  total (35/55)
Fisher's exact test for biased strand distribution p-value = 3.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.02e-01

ATCGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCG  >  NC_000913/3263355‑3263647
                                                                                                                                                  |                                                                                                                                                   
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:1087785/150‑1 (MQ=255)
    gTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGc                                                                                                                                              <  1:487896/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)
         aatGGCGGACTCGTTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑cc                                                                                                                                        <  2:784417/150‑1 (MQ=255)
               ggACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCAtt                                                                                                                                  <  1:1582835/150‑1 (MQ=255)
                gACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCAttt                                                                                                                                 <  2:604106/150‑1 (MQ=255)
                 aCTCGGTAAAAATAGTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCAACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTa                                                                                                                                <  2:616293/150‑1 (MQ=255)
                  cTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGTCAGCGGTGCTGGCTCT‑CCCCCATTTAc                                                                                                                               <  1:1375684/150‑1 (MQ=255)
                   tCGGTAAATATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACg                                                                                                                              <  1:290742/150‑1 (MQ=255)
                    cGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAACTTCCGTTTTTCCACTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACgg                                                                                                                             <  1:484958/150‑1 (MQ=255)
                    cGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACgg                                                                                                                             <  1:35095/150‑1 (MQ=255)
                     ggTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGGGCCAGCGGGGCGGGCTCT‑CCCCCCTTTACggt                                                                                                                            >  1:1806329/1‑149 (MQ=255)
                       tAAAAATACGGCAGCCGTGATCGATGCGGTGGCAATTTAAGGCAACTCTGTAATATAAATTCCGTTTTTACAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGGTGTGGTGCGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATa                                                                                                                          <  2:1629112/150‑1 (MQ=255)
                         aaaaTACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCCCCGGGGCTGGGGCT‑ACCCCCTTTTAGGATaat                                                                                                                        >  2:1579551/1‑149 (MQ=255)
                         aaaaTACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAg                                                                                                                        <  2:526281/150‑1 (MQ=255)
                             tACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATg                                                                                                                    <  2:1296977/150‑1 (MQ=255)
                               cTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCa                                                                                                                  >  1:261464/1‑150 (MQ=255)
                                tGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCAt                                                                                                                 >  2:970837/1‑150 (MQ=255)
                                tGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCAt                                                                                                                 <  2:382645/150‑1 (MQ=255)
                                 gccgccGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCCGCGGTGGTGGGTCG‑CCCCCCATTTCGGATAAGAATGCaat                                                                                                                >  2:1856252/1‑148 (MQ=255)
                                       gTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTCCGa                                                                                                          <  1:1305905/150‑1 (MQ=255)
                                        ggAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTCCGAg                                                                                                         <  2:1601279/149‑1 (MQ=255)
                                            cGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTCCGAGTTaa                                                                                                     <  1:1597817/150‑1 (MQ=255)
                                             gATGCGGTGGCCAATTAAGGCCCCACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAAt                                                                                                    >  1:239586/1‑150 (MQ=255)
                                              aTGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAACTGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATa                                                                                                   <  1:903799/150‑1 (MQ=255)
                                                gCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCTTTTTTCCAGTTCAAATGCAATTGCCTTAAATGCACCTACGTAGCTGTGGTGAGCCAGCGGTTCTGGCTCT‑CCCCCATTTACGCATAAAAATGCATTTTACGAGTTAATaac                                                                                                 >  2:960559/1‑148 (MQ=255)
                                                gCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATAcc                                                                                                 <  2:1296232/150‑1 (MQ=255)
                                                gCGGTGGACAATTAAGGCCACACTGTAATTTAATTTCCGTTTTTCCAGTTCTAATGCAATTGCCTTCAATGCTCCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATAcc                                                                                                 <  1:1439720/150‑1 (MQ=255)
                                                   gTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCGCT‑CCCCCCTTTAAGGTGTAGAATGCATTTTACGAGTTAATACCGTc                                                                                              >  2:886443/1‑150 (MQ=255)
                                                      gCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGc                                                                                           <  2:1874051/150‑1 (MQ=255)
                                                      gCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTGCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGc                                                                                           <  2:1657248/150‑1 (MQ=255)
                                                       ccAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGGGCCAGCGGTGCTGGCTCA‑ACCACCATTAAGGATAAAAAATCAATTTAACAGTTAATAACCTCGGCa                                                                                          >  1:1146909/1‑150 (MQ=255)
                                                       ccAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCa                                                                                          <  1:49064/150‑1 (MQ=255)
                                                       ccAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCa                                                                                          >  1:255935/1‑150 (MQ=255)
                                                       ccAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAACCGTGCTGGGTAG‑CCCCCAGTTAGGGATAAAAATGCATTTTACGCGTTAATACACTCGGcc                                                                                          >  1:721132/1‑149 (MQ=255)
                                                        cAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGGGGT‑CCCCCAGTTTAGGAGAAGAATGCATTTTGCGAGTTTATAACGTGGGCaa                                                                                         >  2:1440727/1‑150 (MQ=255)
                                                        cAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGGGGGGGGTGTCCCCCCATTTA‑GGTTAAAAATGGATTTTGAGAGTTAATAAAGTCGGGaa                                                                                         >  1:1613170/1‑150 (MQ=255)
                                                             aaGGCCAAACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCATTCAATGCACCTTCGTAGCTGTGGTGATCCAGCGGTGCTGGCTCT‑CCCCCAATTACGGAATAGAATTCAATTTCCGAGTTAATACCGTCCGCAATAcca                                                                                    >  2:1025436/1‑149 (MQ=255)
                                                              aGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGTTAAGAATGGATTTTACGAGTTAATAACGTCGGCAATaacag                                                                                   >  2:1861080/1‑146 (MQ=255)
                                                               ggCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTACGAGTTAATACCGTCGGCAATACCTGa                                                                                  >  2:1020407/1‑150 (MQ=255)
                                                               ggCCACACTGTCATTTAAATTCCGTTTTTACAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTCGCTCT‑CCCCCATTTACGGATATGAATGCATTTTCCGAGTTAATACCGTCGGCAATGCCTgc                                                                                  >  2:222022/1‑149 (MQ=255)
                                                                gCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGAc                                                                                 <  1:438780/150‑1 (MQ=255)
                                                                   acacTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCACCGGTGCTGGCTCT‑CCCCCAATTACGGATAAAAATGCATTTTACGAGTTAATACCCTCGGCAATACCAGaaaat                                                                              >  1:1836274/1‑146 (MQ=255)
                                                                   acaaTGTCATTTAAATTCCGTTTTTCCAGTTCTAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑ACCCCATTTACGGATATGAATGCTTTTTCCGAGTTAATACCGTCGGCAATAGCTGACAtt                                                                              <  1:960559/150‑1 (MQ=255)
                                                                    cacTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGGGCCAGCGGGGGTGGCTCTGCCCCCATTTA‑GGTTTAGAAATCATTTTTCGAGTTAAAAAAGGCGGCAATAACAGACAttt                                                                             >  2:1524854/1‑149 (MQ=255)
                                                                      cTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAAt                                                                           <  2:1337300/150‑1 (MQ=255)
                                                                      cTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTACGAGTTAATACCGTCGGCAATACCTGACATTAAt                                                                           >  1:966551/1‑150 (MQ=255)
                                                                      atGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAAt                                                                           <  1:1045806/149‑1 (MQ=255)
                                                                          cATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACtt                                                                       >  1:1342310/1‑150 (MQ=255)
                                                                           aTTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTAG‑ACCCCATTTACGGTTAAGAATGCATTTTACGAGTTAATACCGTCGGCAATACCTGACATTAATACTTc                                                                      >  2:1730299/1‑150 (MQ=255)
                                                                             ttAAATTCCGTTTTTCCATTTCAAATCGAATTGCCTTCAATTCACATTCGTAGCTGTGGTTAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTcac                                                                    <  1:1124129/150‑1 (MQ=255)
                                                                             ttAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTcac                                                                    >  1:1051784/1‑150 (MQ=255)
                                                                                      gTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTgg                                                           <  2:1497783/150‑1 (MQ=255)
                                                                                           tCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCg                                                      <  2:647335/150‑1 (MQ=255)
                                                                                           tCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGGGGGGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCg                                                      <  2:1015798/150‑1 (MQ=255)
                                                                                            caaGTTCAAATGCAATTGCCTTCAATGCTCCGTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGa                                                     <  2:1199429/148‑1 (MQ=255)
                                                                                             cAGTACAACTGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAg                                                    <  2:1760660/150‑1 (MQ=255)
                                                                                              aGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGc                                                   <  1:1388092/150‑1 (MQ=255)
                                                                                                       gCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAAc                                          <  1:1519167/150‑1 (MQ=255)
                                                                                                        cAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCCGCGGTGCTGGCTCT‑CCCCCCTTTAAGGATAAGAATTCATTTTCAGAGTTAATACCGTCGGCAATTACTGACATTTATACTTCACAGGAGCGGGCATCGAGCCCGGaaaaaa                                         >  1:351915/1‑148 (MQ=255)
                                                                                                          aTTGCCTTCAATGAACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTa                                       <  2:1869616/150‑1 (MQ=255)
                                                                                                            tGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAAt                                     <  2:1065263/150‑1 (MQ=255)
                                                                                                             gCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATc                                    <  1:1452759/150‑1 (MQ=255)
                                                                                                             gCCTTCAATGCACCTTCGTAGCTGTGGGGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATc                                    <  2:1284689/150‑1 (MQ=255)
                                                                                                                  cAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTACGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGa                               >  1:4433/1‑150 (MQ=255)
                                                                                                                   aaTGCACCTTCGTAGCTGTGGTGAGCCATCGGTGCTGTGTCT‑CCCCCATTTACGGATAAGAATGCATTTTCCGAGGTAATAACGTCGGCAATACCTGTCATTAATACTTCACAGTCGCTGGCATCTAGCACGGAAAACTTAATCGAAGAc                              <  2:1291286/150‑1 (MQ=255)
                                                                                                                      gCACCTTCGTAGCTGTGGTGGGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGaa                           >  2:1338131/1‑150 (MQ=255)
                                                                                                                        aCCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAAcc                         <  2:239586/150‑1 (MQ=255)
                                                                                                                        aCCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAAcc                         <  1:1080204/150‑1 (MQ=255)
                                                                                                                          cTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACcac                       <  2:261464/150‑1 (MQ=255)
                                                                                                                             cGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGt                    <  1:1793703/150‑1 (MQ=255)
                                                                                                                               tAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTta                  <  2:878946/150‑1 (MQ=255)
                                                                                                                                aGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTtaa                 <  1:1524854/150‑1 (MQ=255)
                                                                                                                                     gggagAGCCAGCACCGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAAcc            >  1:823336/5‑150 (MQ=255)
                                                                                                                                      ggagAGCCAGCACCGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCAAAGTTAATAACCa           >  2:1867647/4‑150 (MQ=255)
                                                                                                                                       gtgAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGCCGAACCACAGTTAATAACCaa          >  1:848399/1‑150 (MQ=255)
                                                                                                                                       gagAGCCAGCACCGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAACACTTCACAGTCGCTGGCATCGAGCACGGAAAACATAATCGAAGACGAACCACAGTtaataaacca          >  1:331575/3‑146 (MQ=255)
                                                                                                                                           gCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCTGCAATACCTGACCTTAATACTTCACAGTCGCTGGCATCGATCACGGAAAACTTAATAGAAGACGATCCACAGTTAATACCCAAcaca      >  2:1121790/1‑150 (MQ=255)
                                                                                                                                           gCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGACAACTTAATCGAAGACGACCCACACTTAATAACCAAaaca      >  2:165897/1‑150 (MQ=255)
                                                                                                                                               gCGGTGCTGGCTCT‑CCCACATTTACGGATAAGAATGCATTTTCCGAGTTAATACAGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCg  <  2:1182319/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
ATCGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCT‑CCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCG  >  NC_000913/3263355‑3263647

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: ‑