Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,604,991 T→A 8.3% L61M (TTG→ATG) ‡ hyfD → hydrogenase 4 component D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,604,9910TA 8.3% 336.2 / 17.2 121L61M (TTG→ATG) ‡hyfDhydrogenase 4 component D
Reads supporting (aligned to +/- strand):  ref base T (52/59);  new base A (5/5);  total (57/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.63e-01

Only 100 of 126 total aligned reads displayed.

TGCCTTTTATCGGCGCACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAAC  >  NC_000913/2604842‑2605140
                                                                                                                                                     |                                                                                                                                                     
tgcCTTTTATCGGCGCACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCacgttgacgt                                                                                                                                                       <  1:1110914/150‑1 (MQ=255)
   cTTTTATCGGCGCACTTGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTgg                                                                                                                                                    <  2:1037687/150‑1 (MQ=255)
   cTTTTATCGGCGCACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTgg                                                                                                                                                    <  1:1012912/150‑1 (MQ=255)
    ttttATCGGCGCACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGt                                                                                                                                                   <  2:1993678/150‑1 (MQ=255)
        aTCTGCGCACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAAc                                                                                                                                               >  1:2115354/1‑150 (MQ=255)
         tCGGCGCACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACg                                                                                                                                              >  1:1962691/1‑150 (MQ=255)
           ggCGCACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACCAACGTCAACGTg                                                                                                                                            <  1:1759501/150‑1 (MQ=255)
            gcgcACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTgg                                                                                                                                           >  2:2152433/1‑150 (MQ=255)
             ggctcTGGTCGTTTCGTTTTCGGGACTAGGGGGGGCGGCCGAAGGGGGGGTTTTGTGCGCCGCGCTGACCTCTCTGTGCTTGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTAACTTGGACGAACGTAAACGTgac                                                                                                                                          <  2:1677844/146‑3 (MQ=255)
             cgcACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTggg                                                                                                                                          <  1:972266/150‑1 (MQ=255)
             cgcACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTggg                                                                                                                                          <  2:369563/150‑1 (MQ=255)
               cACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTggggg                                                                                                                                        <  2:2381390/150‑1 (MQ=255)
               cACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTggggg                                                                                                                                        >  2:1660854/1‑150 (MQ=255)
               cACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACCAACGTCAACGTgacgg                                                                                                                                        <  1:1461862/150‑5 (MQ=255)
               cACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACCAACGTCAACGTgacgg                                                                                                                                        <  1:1166255/150‑5 (MQ=255)
                 cTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTCGGTCAACGTGGGGGAt                                                                                                                                      >  1:1577307/1‑150 (MQ=255)
                    gTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACAACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGAtgtg                                                                                                                                   >  2:1953684/1‑150 (MQ=255)
                        tttcgtttTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCgt                                                                                                                               >  1:358565/1‑150 (MQ=255)
                        tttcgtttTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCgt                                                                                                                               <  1:318815/150‑1 (MQ=255)
                        tttcgtttTCGCCACAACGTCGGGCCGCAGAATGCGCGCTTTTCTTCACCGCGCTTACCACGCTGTTCATGTTGTCACTGATCTCCCCCTATTATCAGCCCGATAAACATCCCGTCACGTTGACGTTTCTCAACGTGCGGGAtgtcgagt                                                                                                                               >  2:62789/1‑145 (MQ=255)
                          tcgtttTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCgttg                                                                                                                             <  2:951963/150‑1 (MQ=255)
                           cgtttTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACCAACGTCAACGTgac                                                                                                                                          <  1:702185/136‑3 (MQ=255)
                           cgtttTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACCAACGTCAACGTgac                                                                                                                                          >  2:702185/1‑134 (MQ=255)
                            gtttTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGAtgtg                                                                                                                                   <  1:2115777/142‑1 (MQ=255)
                            gtttTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGAtgtg                                                                                                                                   >  2:2115777/1‑142 (MQ=255)
                             ttttTGCCACAACGTCGGGCCGCCGTATGGGGGGTTTTGTTTGCCCCGGTGTCCAAGCTTTTCATGTTGTCACTGTTATTCTTTTTTTATCAGGCCGTTAAAGTTGCCGTTACGTTGACGTTGGTCAACGTTGGGTATGTGGCGTTGttt                                                                                                                          <  2:1636740/150‑1 (MQ=255)
                              tttCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTg                                                                                                                         <  2:1449135/150‑1 (MQ=255)
                               ttCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTgg                                                                                                                        <  2:614398/150‑1 (MQ=255)
                               ttCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGGCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTgg                                                                                                                        <  2:2383788/150‑1 (MQ=255)
                                tCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGc                                                                                                                       >  1:835854/1‑150 (MQ=255)
                                   ccACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTg                                                                                                                    >  1:207258/1‑150 (MQ=255)
                                   ccACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTg                                                                                                                    <  2:19366/150‑1 (MQ=255)
                                    cacaACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTgg                                                                                                                   <  1:1455962/150‑1 (MQ=255)
                                    cacaAAGTCGGGCCGCGGAATGGGGGGTTATGTTCGCCGCGCTGAGCACGCTGTGCCTGTTGTCTCTGATCACCGAGTGTTATCAGGCGGTTAAAGTTGCCGTCTCGTTGAGGTTGGTCAACGTGGGGGATGTGGCGTTGTGTGGCCTgg                                                                                                                   <  2:2336526/150‑1 (MQ=255)
                                             gggCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATc                                                                                                          >  1:519452/1‑150 (MQ=255)
                                                  gccgAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATGGCGTg                                                                                                     <  1:1266839/150‑1 (MQ=255)
                                                  gccgAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTg                                                                                                     <  1:2335370/150‑1 (MQ=255)
                                                  gccgAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTg                                                                                                     >  2:460646/1‑150 (MQ=255)
                                                     gAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGt                                                                                                  >  1:419620/1‑150 (MQ=255)
                                                     gAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGt                                                                                                  <  1:1962427/150‑1 (MQ=255)
                                                       aTGGGGGGGGTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTAc                                                                                                <  2:1988957/150‑1 (MQ=255)
                                                        tGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACg                                                                                               <  2:2124646/150‑1 (MQ=255)
                                                               ttttGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTGCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTc                                                                                        <  1:1350139/150‑1 (MQ=255)
                                                               ttttGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTc                                                                                        <  1:1076130/150‑1 (MQ=255)
                                                               ttttGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTc                                                                                        >  2:2393760/1‑150 (MQ=255)
                                                                 ttGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTg                                                                                      >  2:1713916/1‑150 (MQ=255)
                                                                   gTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGtt                                                                                    <  1:1904091/150‑1 (MQ=255)
                                                                       gccgcGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTtgtg                                                                                >  1:48975/1‑150 (MQ=255)
                                                                           cgcTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTgtggtgg                                                                            <  2:309600/150‑1 (MQ=255)
                                                                           cgcTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTgtggtgg                                                                            <  2:2115354/150‑1 (MQ=255)
                                                                                aCCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGttt                                                                       >  2:441367/1‑150 (MQ=255)
                                                                                   aCGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTtctc                                                                    <  1:460646/150‑1 (MQ=255)
                                                                                    cGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTTGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATAGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCg                                                                   >  1:279116/1‑150 (MQ=255)
                                                                                    cGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCg                                                                   >  1:1032208/1‑150 (MQ=255)
                                                                                      cTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGt                                                                 >  2:208864/1‑150 (MQ=255)
                                                                                          gCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGc                                                             >  2:457539/1‑150 (MQ=255)
                                                                                           cATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCt                                                            <  1:841619/150‑1 (MQ=255)
                                                                                           cATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCt                                                            >  2:607524/1‑150 (MQ=255)
                                                                                           cATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCt                                                            <  2:2093816/150‑1 (MQ=255)
                                                                                            aTGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTg                                                           <  2:960879/150‑1 (MQ=255)
                                                                                             tgttgtCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTgg                                                          <  1:187215/150‑1 (MQ=255)
                                                                                                 gtCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCAc                                                      <  1:133345/150‑1 (MQ=255)
                                                                                                   cACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGa                                                    >  1:2089343/1‑150 (MQ=255)
                                                                                                     cTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATc                                                  <  1:62789/150‑1 (MQ=255)
                                                                                                     cTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATc                                                  <  1:1457692/150‑1 (MQ=255)
                                                                                                       gATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTa                                                >  1:734902/1‑150 (MQ=255)
                                                                                                       gATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTa                                                <  1:1094889/150‑1 (MQ=255)
                                                                                                         tctcCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACt                                              <  2:2268425/150‑1 (MQ=255)
                                                                                                          ctcCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTc                                             <  2:1496131/150‑1 (MQ=255)
                                                                                                            ggcaactttATCGGCCTGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTATGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTccc                                           >  2:294547/7‑149 (MQ=255)
                                                                                                             cgcgTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCAc                                          >  1:1188890/1‑150 (MQ=255)
                                                                                                              gcgTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACg                                         <  2:1962691/150‑1 (MQ=255)
                                                                                                              caactttATCGGCCTGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACg                                         >  2:1794227/5‑150 (MQ=255)
                                                                                                               cgTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACgg                                        <  2:680855/150‑1 (MQ=255)
                                                                                                               cgTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACgg                                        <  2:1429698/150‑1 (MQ=255)
                                                                                                               cgTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACgg                                        <  1:441367/150‑1 (MQ=255)
                                                                                                                 ttttATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGt                                      >  1:1325272/1‑150 (MQ=255)
                                                                                                                 ttttATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGt                                      <  1:1765276/150‑1 (MQ=255)
                                                                                                                  tttATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGtt                                     <  2:1508864/150‑1 (MQ=255)
                                                                                                                      tCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCt                                 <  2:1887181/150‑1 (MQ=255)
                                                                                                                      tCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCt                                 <  1:614754/150‑1 (MQ=255)
                                                                                                                          gCCGATAATGTTGCCGTCACGTTGTGGTTGGTCAAGGTGGGGGTTGTGGCGTTGTTTGGCCTGGTCATTGATCGAGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTATGTTTGATTTTCACGATCTAGTCCACGTGTTATCTGACg                             <  2:966231/150‑1 (MQ=255)
                                                                                                                          gCCGATAAAGTTGCCGTCACGTTGACGTCGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACg                             <  2:1577307/150‑1 (MQ=255)
                                                                                                                            cGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGa                           >  1:886063/1‑150 (MQ=255)
                                                                                                                             gATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGAt                          >  2:337171/1‑150 (MQ=255)
                                                                                                                                 aaGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATaaaa                      >  2:2389311/1‑150 (MQ=255)
                                                                                                                                 aaGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATaaaa                      >  1:363198/1‑150 (MQ=255)
                                                                                                                                   gTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAAt                    >  1:1153902/1‑150 (MQ=255)
                                                                                                                                     tGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATcg                  <  1:8745/150‑1 (MQ=255)
                                                                                                                                     tGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATcg                  >  2:696437/1‑150 (MQ=255)
                                                                                                                                       ccGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATcgcg                <  1:208864/150‑1 (MQ=255)
                                                                                                                                       ccGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATcgcg                >  1:1911662/1‑150 (MQ=255)
                                                                                                                                        ccccACGTTGACCAACGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGa               >  2:2293452/4‑150 (MQ=255)
                                                                                                                                        cGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGa               <  2:1589201/150‑1 (MQ=255)
                                                                                                                                           cACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAaca            >  2:1779035/1‑150 (MQ=255)
                                                                                                                                           cACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAaca            <  2:734902/150‑1 (MQ=255)
                                                                                                                                           cACGTTGACCAACGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAaca            >  1:1770442/1‑150 (MQ=255)
                                                                                                                                               ttgacgttgGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCg        <  2:419620/150‑1 (MQ=255)
                                                                                                                                                 gacgttgGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCa      <  1:2363426/150‑1 (MQ=255)
                                                                                                                                                     ttgGTAAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAAc  >  2:82887/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGCCTTTTATCGGCGCACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAAC  >  NC_000913/2604842‑2605140

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