Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,263,518 T→A 9.5% N48I (AAT→ATT)  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,5180TA 9.5% 1156.5 / 15.7 393N48I (AAT→ATT) tdcDpropionate kinase/acetate kinase C, anaerobic
Reads supporting (aligned to +/- strand):  ref base T (116/216);  new base A (47/2);  total (174/219)
Fisher's exact test for biased strand distribution p-value = 2.90e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 395 total aligned reads displayed.

CGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTCATTCATCTCTTTT  >  NC_000913/3263369‑3263667
                                                                                                                                                     |                                                                                                                                                     
     tCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACg                                                                                                                                                  <  1:1264260/150‑1 (MQ=255)
     tCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACg                                                                                                                                                  <  1:1425946/150‑1 (MQ=255)
       ggTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGa                                                                                                                                                <  2:6778086/150‑1 (MQ=255)
        gTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGAt                                                                                                                                               <  1:6346626/150‑1 (MQ=255)
        gTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCAATTAAGGAt                                                                                                                                               >  2:168915/1‑150 (MQ=255)
          aaaaaTACTGCCGCCGTGAGCGATGCGGTGGCCATTTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATaa                                                                                                                                             <  1:4566284/150‑1 (MQ=255)
              aTACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAAt                                                                                                                                         <  2:4391069/150‑1 (MQ=255)
              aTACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAAt                                                                                                                                         <  1:5614638/150‑1 (MQ=255)
              aTACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAAt                                                                                                                                         <  1:5850840/150‑1 (MQ=255)
                aCTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGc                                                                                                                                       <  2:376334/150‑1 (MQ=255)
                   gccgccGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGTGCCCGCGGGGCTGGCGGGCCCCCAATTACGGCTAAGAAAGCAtt                                                                                                                                    >  1:5504402/1‑150 (MQ=255)
                    ccgccgTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGTGCCAGCGGTGCTGGCTCTCCCCCCGTTACGGATAAGAATGCAatt                                                                                                                                   >  1:4167255/1‑150 (MQ=255)
                    ccgccgTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCCCTGCTTGCTCTCCCCCCTTTACGGATAATAAAGCAttt                                                                                                                                   >  2:5671591/1‑150 (MQ=255)
                    ccgccgTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCACCGGTGCTGGGGCTACCCCAATTAAGGATAAGAATTCAatt                                                                                                                                   >  1:3337784/1‑150 (MQ=255)
                     cgccgTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGGTCTCCCCCAAATACGGATAAGAATGCAtttt                                                                                                                                  >  1:373484/1‑150 (MQ=255)
                     cgccgTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGGGCTGGCTCTCCCCCCTTTACGGATAAGAATGCAtttt                                                                                                                                  >  1:2420930/1‑150 (MQ=255)
                      gccgTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTc                                                                                                                                 >  2:4210329/1‑150 (MQ=255)
                        cgTGAGCGATGCGGTGGCCAATTCAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCg                                                                                                                               <  2:4877782/150‑1 (MQ=255)
                         gTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCTCCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGa                                                                                                                              <  2:102832/150‑1 (MQ=255)
                         gTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCTATTTACGGATAAGAATGCATTTTCCGa                                                                                                                              <  2:5473544/150‑1 (MQ=255)
                          tGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGGGCCCGCGGTGCGGGCTCACCCCCCGTTACGGAGAACAATGCAATTTCCCAg                                                                                                                             >  1:2766376/1‑150 (MQ=255)
                          tGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCAATTAAGGATAAGAATGCAttttaagag                                                                                                                             >  1:5810963/1‑145 (MQ=255)
                           gAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGt                                                                                                                            <  2:3011230/150‑1 (MQ=255)
                              cGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTaa                                                                                                                         <  1:5071498/150‑1 (MQ=255)
                                    ggTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTATTTCCAGTTCAAATGCAATTGCCTTCAATGCACCATCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGt                                                                                                                   <  2:4700422/150‑1 (MQ=255)
                                      tGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCCCCGGGGCTGGGTCTCCCCCCATTAAGGGTAAGAATTGAATTTCCGAGTTAATaaccacg                                                                                                                 >  1:2234072/1‑144 (MQ=255)
                                       ggCCAATTAAGGCCACACTGTCATTTAAAGTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATAACGTCgg                                                                                                                <  2:115871/150‑1 (MQ=255)
                                        gCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGc                                                                                                               <  2:6436950/150‑1 (MQ=255)
                                        gCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGc                                                                                                               <  1:3462668/150‑1 (MQ=255)
                                        gCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCCGTGCGGGCACTCCCCCCAATACGGATAAGAAAGCAATTTCCCAGTTAATACCGACGGc                                                                                                               >  2:7053515/1‑150 (MQ=255)
                                          cAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGGTCCAGCGGTGCTTGCTCTCCCCCAATTAAGGATAAGAAAGCATTTTACGAGTTTATAACGACGGCaa                                                                                                             >  2:269170/1‑150 (MQ=255)
                                          cAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCaa                                                                                                             <  2:5851040/150‑1 (MQ=255)
                                          cAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCaa                                                                                                             <  2:2248924/150‑1 (MQ=255)
                                            aTTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATa                                                                                                           >  2:6007748/1‑150 (MQ=255)
                                             ttAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATAc                                                                                                          <  2:3507330/150‑1 (MQ=255)
                                                 ggCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGa                                                                                                      <  1:6716462/150‑1 (MQ=255)
                                                 ggCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGa                                                                                                      <  2:6731105/150‑1 (MQ=255)
                                                 ggCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGa                                                                                                      <  1:4372110/150‑1 (MQ=255)
                                                    cacaCTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTAAGGATAAGAATGCATTTTACGAGTTAATAACGTCGGCAAAACCAGACAt                                                                                                   >  2:4946420/1‑150 (MQ=255)
                                                       acTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGGGCCAGCGGGGCTGGGGCTGCCCCAAATAAGGAAAATAAAGCATATTGCGAGTTAAAAACGTCGGCAAAAACAGACaaaaa                                                                                                >  1:317284/1‑146 (MQ=255)
                                                       acTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTaa                                                                                                <  1:5931266/150‑1 (MQ=255)
                                                              tttAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCCGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTGCGAGTTAATACCGTCGGCAATACCAGACAATAATACttaa                                                                                         >  2:7003175/1‑148 (MQ=255)
                                                              tttAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCACCGGGGCTGGCTCTCCCCCCATTACGGGTAAGAATTCAGTTTACGAATTAATAACGTCGGCAAAAACAGACATTAATACttaa                                                                                         >  2:1274596/1‑148 (MQ=255)
                                                                tAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTcaca                                                                                       <  2:1481488/150‑1 (MQ=255)
                                                                   aTTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTc                                                                                    <  1:6863569/150‑1 (MQ=255)
                                                                     tCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGc                                                                                  <  1:6646738/150‑1 (MQ=255)
                                                                     tCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGc                                                                                  <  1:2312615/150‑1 (MQ=255)
                                                                      ccGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCt                                                                                 <  2:5195908/150‑1 (MQ=255)
                                                                      ccGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCt                                                                                 <  1:4483797/150‑1 (MQ=255)
                                                                        gTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGGGCCCGCGGGGCGGGCGCGACCCCAAATAAGGATAAGAAAGGAATTTCCGAGTTAATAACGGCGGGAAAAACGATCATTAAAAAATCACAGTCGCggg                                                                               >  2:2501960/1‑150 (MQ=255)
                                                                        gTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCCGCGGGGCTTGCAGTCCCCCATTTAAGGATAAGAATGCAATTTACGATTTAATAACCTCGGCAAAAACTGAAAATAATAATTAAAAGTAGCTgg                                                                               >  2:852786/1‑150 (MQ=255)
                                                                        gTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTgg                                                                               <  1:6962768/150‑1 (MQ=255)
                                                                        gTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTgg                                                                               <  1:4930421/150‑1 (MQ=255)
                                                                        gTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTgg                                                                               <  2:3424884/150‑1 (MQ=255)
                                                                          ttttCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCa                                                                             <  2:361892/150‑1 (MQ=255)
                                                                            ttCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATc                                                                           <  2:3917432/150‑1 (MQ=255)
                                                                              ccAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGa                                                                         <  1:1793195/150‑1 (MQ=255)
                                                                               cAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAg                                                                        <  1:5300390/150‑1 (MQ=255)
                                                                               cAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAg                                                                        <  2:3258520/150‑1 (MQ=255)
                                                                                aGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGc                                                                       <  1:6494402/150‑1 (MQ=255)
                                                                                 gTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCa                                                                      <  2:6549521/150‑1 (MQ=255)
                                                                                  ttCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCCATACCTGACATTAATACTTCACAGTCGCGGGCATCGAGccc                                                                     >  2:5985374/1‑148 (MQ=255)
                                                                                   tCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACg                                                                    <  1:3100561/150‑1 (MQ=255)
                                                                                   tCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACg                                                                    <  2:4098518/150‑1 (MQ=255)
                                                                                             ttGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGGGCTGGCTCTCCCCCATTTAACGAGAAGAATGCATTTTCCGAGGTAATACCCTCGGCAAAACCTGACAATAATACATAACAGTCGCGGGCCTCGAGCACCGAAAACTTaa                                                          >  2:7131968/1‑150 (MQ=255)
                                                                                               gCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATc                                                        <  2:2068540/150‑1 (MQ=255)
                                                                                                ccTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCg                                                       <  2:2193110/150‑1 (MQ=255)
                                                                                                  ttCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATAACTGACATTAATAATTCACAGTCGCTGGCATCGAGCCCGGAAAAAATAATCGaa                                                     >  2:2376312/1‑150 (MQ=255)
                                                                                                      aTGCACCTTCGTAGCTGTGGTGGGCCAGCGGTGGTGGGTGTACCCCATTTACGGATAATAATTCATTTTGCGTGTTAATACCGTGGGCAATACCTGTCATTTATTATTGTCATTGGGTGGGATTGAGGTCGGAAAAATTTATAGAAGACg                                                 >  2:2772624/1‑150 (MQ=255)
                                                                                                        gCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGaa                                               >  2:5326822/1‑150 (MQ=255)
                                                                                                             ttCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACcaca                                          <  2:6303208/150‑1 (MQ=255)
                                                                                                                gTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGtt                                       <  2:4321447/150‑1 (MQ=255)
                                                                                                                   gCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTtaat                                    <  2:2959249/150‑1 (MQ=255)
                                                                                                                     tgtgGTGCGCCAGCGGTGCGGGCTCTCCCCCATTTACGGATAAGAAAGCATTTTCCGAGTTAATACCGGCGGCAATACCTGACATTAATACTTCACAGTCGCGGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGGgaataa                                  >  2:744600/1‑150 (MQ=255)
                                                                                                                      gtggtgAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAAc                                 <  2:4481197/150‑1 (MQ=255)
                                                                                                                      gtggtgAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAAc                                 <  1:497118/150‑1 (MQ=255)
                                                                                                                        ggagAGCCAGCACCGCTGGCTCTCCCCCATTTACGGATACGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGCAGACGACCCACAGTTAATAAccc                               >  2:4812477/4‑149 (MQ=255)
                                                                                                                        ggagAGCCAGCACCGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCa                               >  2:6064610/4‑150 (MQ=255)
                                                                                                                        ggagAGCCAGCACCGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCa                               >  2:6154625/4‑150 (MQ=255)
                                                                                                                          tgAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCaaa                             <  1:5117464/150‑1 (MQ=255)
                                                                                                                          tgAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGAATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCaaa                             <  1:968294/150‑1 (MQ=255)
                                                                                                                           gAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCaaaa                            <  1:4643886/150‑1 (MQ=255)
                                                                                                                           gAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCaaaa                            <  2:663337/150‑1 (MQ=255)
                                                                                                                             gCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAaaca                          >  1:1944645/1‑150 (MQ=255)
                                                                                                                                aGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAAcc                       >  2:1110372/1‑150 (MQ=255)
                                                                                                                                    gTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGaa                   >  2:2805659/1‑150 (MQ=255)
                                                                                                                                     tggtGGCTCTCCCCATTTTAGGGATATGAATGCTTTTTCCGAGTTAATACCGTCGGCATTACCTGAATTTATTAGTTCACAGTCGGTGGAATGGTGAACGGAAAACTTATTCGAATACGAACAACTGTTAATAACCAAAACAACCGGaaa                  <  1:5435030/147‑1 (MQ=255)
                                                                                                                                      gCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAat                 <  2:1375843/150‑1 (MQ=255)
                                                                                                                                        tGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAattc               >  1:1028741/1‑150 (MQ=255)
                                                                                                                                        tGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAattc               >  1:2985113/1‑150 (MQ=255)
                                                                                                                                        tGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAattc               >  1:6273898/1‑150 (MQ=255)
                                                                                                                                         ggCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGACattca              >  1:2008950/1‑150 (MQ=255)
                                                                                                                                         ggCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAattca              <  2:5104926/150‑1 (MQ=255)
                                                                                                                                          gCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAattcat             <  2:3038810/150‑1 (MQ=255)
                                                                                                                                            tctcCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAattcattc           >  2:6891179/1‑150 (MQ=255)
                                                                                                                                            tctcCCCCATTTACGGATAAGAATGCATTTTCAGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAaatcattc           >  1:3960613/1‑150 (MQ=255)
                                                                                                                                             ctcCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAAttcattca          >  2:2560276/1‑150 (MQ=255)
                                                                                                                                                ccccATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTCATTCAtct       >  1:1521678/1‑150 (MQ=255)
                                                                                                                                                 cccATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACAGTAAACTTAATCGTAGACGTACCACAGTTAATTACCAAAACAACCGGAAATTCATTCAtgtc      >  2:1930866/1‑150 (MQ=255)
                                                                                                                                                     tttACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACATTTAATAACCAAAACAACCGGAAATTCATTCATCTCtttt  <  1:7164075/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTCATTCATCTCTTTT  >  NC_000913/3263369‑3263667

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: