Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,542,928 C→T 6.0% A201T (GCG→ACG)  cysP ← thiosulfate/sulfate ABC transporter periplasmic binding protein CysP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,542,9280CT 6.0% 362.9 / 6.9 116A201T (GCG→ACG) cysPthiosulfate/sulfate ABC transporter periplasmic binding protein CysP
Reads supporting (aligned to +/- strand):  ref base C (49/59);  new base T (2/5);  total (51/65)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 120 total aligned reads displayed.

CGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGG  >  NC_000913/2542782‑2543076
                                                                                                                                                  |                                                                                                                                                    
cGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGcc                                                                                                                                                   <  1:1042268/150‑1 (MQ=255)
 gCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCa                                                                                                                                                  <  1:97866/150‑1 (MQ=255)
 gCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCa                                                                                                                                                  <  2:1096660/150‑1 (MQ=255)
  ccACCGGGAATTCCGACAGAATGTTCGATTTCGTAATCACCACTTCAAACCCCTGCTCTTCATACTGTTTACGGATGTTGTTAACTTCCGATTCGAAGCTAATCAGCACATCGCCAAGGCCGCGCTCGCCAAAAGTCGTGGTCGCGCCAc                                                                                                                                                 >  2:902710/1‑150 (MQ=255)
   cACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTCACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACg                                                                                                                                                >  1:177759/1‑150 (MQ=255)
     ccGGGAATTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGAc                                                                                                                                              >  1:1067085/1‑150 (MQ=255)
      cGGGAATTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGAcc                                                                                                                                             <  1:1480519/150‑1 (MQ=255)
          aaTTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCa                                                                                                                                         <  2:69375/150‑1 (MQ=255)
          aaTTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGCGGTCGTGGCGCGAccacca                                                                                                                                         <  1:1475724/150‑1 (MQ=255)
           aTTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGCGGTCGTGGCGCGAccaccac                                                                                                                                        <  1:22207/150‑2 (MQ=255)
           aTTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGCGGTCGTGGCGCGAccaccac                                                                                                                                        <  1:389813/150‑2 (MQ=255)
              ccGCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTAt                                                                                                                                     <  1:938167/150‑1 (MQ=255)
               cGCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATc                                                                                                                                    <  2:373869/150‑1 (MQ=255)
                gCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCg                                                                                                                                   <  2:1493094/150‑1 (MQ=255)
                gCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCg                                                                                                                                   <  1:1028174/150‑1 (MQ=255)
                 ccAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGa                                                                                                                                  <  1:1348522/150‑1 (MQ=255)
                 ccAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGa                                                                                                                                  <  1:1410329/150‑1 (MQ=255)
                     aaTGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAacac                                                                                                                              <  1:544434/150‑1 (MQ=255)
                     aaTGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAacac                                                                                                                              >  2:217126/1‑150 (MQ=255)
                       tGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACtt                                                                                                                            <  1:488947/150‑1 (MQ=255)
                        gTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTc                                                                                                                           >  2:898922/1‑150 (MQ=255)
                          tGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCaa                                                                                                                         >  2:902349/1‑150 (MQ=255)
                              tttCGGAATCAACACTTCAAAGCCCTGCGCTTCATACTGTTTACGGAAGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAGCACTTCAACGtt                                                                                                                     >  2:484752/1‑150 (MQ=255)
                                tCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGTGGCGCGAccaccact                                                                                                                                       <  2:132401/130‑3 (MQ=255)
                                tCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGtttt                                                                                                                   <  1:239251/150‑1 (MQ=255)
                                tCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGtttt                                                                                                                   >  2:362927/1‑150 (MQ=255)
                                tCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGtttt                                                                                                                   <  2:570381/150‑1 (MQ=255)
                                  ggAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTc                                                                                                                 <  2:672806/150‑1 (MQ=255)
                                  ggAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGGCGCGCCACGACCGCCAGTATCGAACACTTCATCGTTTTTc                                                                                                                 <  2:521313/150‑1 (MQ=255)
                                  ggAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTGGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTc                                                                                                                 <  1:1324088/150‑1 (MQ=255)
                                    aaTCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAg                                                                                                               <  1:1376252/150‑1 (MQ=255)
                                    aaTCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAg                                                                                                               <  1:1045650/150‑1 (MQ=255)
                                        accacTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAAc                                                                                                           >  2:1127897/1‑150 (MQ=255)
                                        accacTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAAc                                                                                                           >  2:1122415/1‑150 (MQ=255)
                                        accacTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCAAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAAc                                                                                                           >  2:939201/1‑150 (MQ=255)
                                        accacTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCAAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAAc                                                                                                           >  2:944925/1‑150 (MQ=255)
                                         ccacTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACt                                                                                                          >  2:1285844/1‑150 (MQ=255)
                                         ccacTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACt                                                                                                          >  1:703376/1‑150 (MQ=255)
                                             ttCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGt                                                                                                      <  1:1186719/150‑1 (MQ=255)
                                              tCACAGCCCTGCGCGTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTATTCTGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGCGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTc                                                                                                     <  2:1377862/150‑1 (MQ=255)
                                              tCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTc                                                                                                     <  2:1067085/150‑1 (MQ=255)
                                                    cccTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAAc                                                                                               <  1:1403385/150‑1 (MQ=255)
                                                        gcgcTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGtt                                                                                           >  2:1301434/1‑150 (MQ=255)
                                                               aTACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGtttt                                                                                    >  2:1031939/1‑150 (MQ=255)
                                                                  cTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGcc                                                                                 <  1:950093/150‑1 (MQ=255)
                                                                      ttACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTg                                                                             <  2:434490/150‑1 (MQ=255)
                                                                      ttACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTg                                                                             <  1:902349/150‑1 (MQ=255)
                                                                      ttACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTg                                                                             <  2:60430/150‑1 (MQ=255)
                                                                       tACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATAGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGACACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGt                                                                            >  1:432587/1‑150 (MQ=255)
                                                                        aCGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTc                                                                           <  2:78286/150‑1 (MQ=255)
                                                                           gATGTTGTTCACTTCCGATTCGAAGCTAATCGGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTcacc                                                                        <  1:12128/150‑1 (MQ=255)
                                                                             agttgttCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTATTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTcaccac                                                                      >  1:844006/2‑150 (MQ=255)
                                                                              gttgttCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGGACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCaccacc                                                                     >  1:9080/1‑150 (MQ=255)
                                                                                  ttCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCCCCACCGTcc                                                                 >  2:64484/1‑149 (MQ=255)
                                                                                  ttCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCa                                                                 >  1:1208715/1‑150 (MQ=255)
                                                                                   tCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAg                                                                >  2:1084585/1‑150 (MQ=255)
                                                                                   tCAATTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAg                                                                >  2:766687/1‑150 (MQ=255)
                                                                                     aCTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAACAGTGGTGGTCGCGCCACGGCCGCCAGTATCGGACACTGCAGCGTTGTTCGGGAGCTGGGTCATAAACGGTTCGGGTTTGCCTTTGTGGGCGCCGTCAggt                                                              >  2:1537087/1‑148 (MQ=255)
                                                                                     aCTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACt                                                                                              >  1:843850/1‑118 (MQ=255)
                                                                                     aCTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACt                                                                                              <  2:843850/118‑1 (MQ=255)
                                                                                      cTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACTCTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCtt                                                             >  2:630658/1‑150 (MQ=255)
                                                                                        tCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTa                                                           <  1:903420/150‑1 (MQ=255)
                                                                                          cGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATc                                                         >  1:602095/1‑150 (MQ=255)
                                                                                               cGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTg                                                    >  2:290815/1‑150 (MQ=255)
                                                                                               cGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTg                                                    >  2:627165/1‑150 (MQ=255)
                                                                                                gAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTgc                                                   >  2:1384102/1‑150 (MQ=255)
                                                                                                    cTAATCATCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGccc                                               <  1:885416/150‑1 (MQ=255)
                                                                                                       aTCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCAt                                            >  2:273994/1‑150 (MQ=255)
                                                                                                        tCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGAATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATg                                           <  2:177759/150‑1 (MQ=255)
                                                                                                             acaTCGCAGATGCCGAGCTCGGCTATTGTGGTTGTCGATCCACGAGGGGGAGTATCTAACTATTCTAAGTTTTTCAGGAACTTTTTAATTAACTGTTCGGTTTTGCCTTTCTCACCAACGTCAGCTTTTTCCTCTGCGGCCCTTgccggc                                      <  2:495966/150‑3 (MQ=255)
                                                                                                              caTCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCa                                     >  2:1163382/1‑150 (MQ=255)
                                                                                                               aTCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAg                                    <  1:79041/150‑1 (MQ=255)
                                                                                                                tCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGa                                   >  2:1239858/1‑150 (MQ=255)
                                                                                                                  gCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATa                                 <  1:939201/150‑1 (MQ=255)
                                                                                                                  gCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATa                                 <  1:944925/150‑1 (MQ=255)
                                                                                                                    ccAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAgg                               <  1:217126/150‑1 (MQ=255)
                                                                                                                     cAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAgg                               >  2:18983/1‑149 (MQ=255)
                                                                                                                     cAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAgg                               <  1:18983/149‑1 (MQ=255)
                                                                                                                     cAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGt                              <  1:1084585/150‑1 (MQ=255)
                                                                                                                      aGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGta                             <  1:362927/150‑1 (MQ=255)
                                                                                                                       ggCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCTCATGCCGCCAGATAGGtat                            >  1:1569810/1‑150 (MQ=255)
                                                                                                                       ggCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGtat                            <  1:1572290/150‑1 (MQ=255)
                                                                                                                       ggCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGtat                            >  1:432285/1‑150 (MQ=255)
                                                                                                                       ggCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGtat                            >  1:779517/1‑150 (MQ=255)
                                                                                                                         ccGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGtat                            <  2:343256/148‑1 (MQ=255)
                                                                                                                         ccGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGtat                            >  1:343256/1‑148 (MQ=255)
                                                                                                                            cgcTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAAcgc                       >  2:392690/1‑150 (MQ=255)
                                                                                                                              cTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAAcgcgc                     <  1:290815/150‑1 (MQ=255)
                                                                                                                                  gCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTa                 >  1:56956/1‑150 (MQ=255)
                                                                                                                                   cAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTAc                >  1:342268/1‑150 (MQ=255)
                                                                                                                                     aaaGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTAccc              <  1:1031939/150‑1 (MQ=255)
                                                                                                                                       aGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGa            <  1:201083/150‑1 (MQ=255)
                                                                                                                                         tggtggTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACg          <  1:578967/150‑1 (MQ=255)
                                                                                                                                          ggcggtCGTGGCGCGACCACCTGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGt         >  2:1025157/4‑150 (MQ=255)
                                                                                                                                             ggtCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGtttt      >  1:1357273/1‑150 (MQ=255)
                                                                                                                                               tCGCGCCACGACCGCCAGTATCGGACACTTCAACGTTTTTCAGTAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTttcg    <  2:9080/150‑1 (MQ=255)
                                                                                                                                                cgcgCCACGCCCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCGTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTttcgg   >  2:121662/1‑150 (MQ=255)
                                                                                                                                                cgcgCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTttcgg   <  1:1301434/150‑1 (MQ=255)
                                                                                                                                                 gcgcCACGACCGCCAGTATCGTACACTTCAACGTTTTTCAGGAACTGGGTCATAATCTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGGCAGATTGGTATAACGCGCGTTACCCGACGTTttcggg  <  2:155272/150‑1 (MQ=255)
                                                                                                                                                 gcgcCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTttcggg  <  2:376746/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGG  >  NC_000913/2542782‑2543076

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: