Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,661,941 C→G 8.0% V462V (GTC→GTG) ‡ ynfF → putative selenate reductase YnfF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,661,9410CG 8.0% 279.2 / 10.8 85V462V (GTC→GTG) ‡ynfFputative selenate reductase YnfF
Reads supporting (aligned to +/- strand):  ref base C (43/35);  new base G (3/3);  total (47/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.26e-01

AGCTGGGATCTGGGGGTAGAATGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTG  >  NC_000913/1661792‑1662090
                                                                                                                                                     |                                                                                                                                                     
aGCTGGGATCTGGGGGTAGAATGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTc                                                                                                                                                       >  1:1535154/1‑150 (MQ=255)
 gCTGGGATCTGGGGGTAGAATGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTcc                                                                                                                                                      >  2:1480010/1‑150 (MQ=255)
    gggATCTGGGGGTAGAATGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCa                                                                                                                                                   <  1:505832/150‑1 (MQ=255)
     ggATCTGGGGGTAGAATGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCAt                                                                                                                                                  <  2:155455/150‑1 (MQ=255)
      gATCTGGGGGTAGAATGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATc                                                                                                                                                 >  1:869573/1‑150 (MQ=255)
       aTCTGGGGGTAGAATGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCa                                                                                                                                                <  1:371483/150‑1 (MQ=255)
           gggggTAGAATGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGtt                                                                                                                                            >  2:1476816/1‑150 (MQ=255)
           gggggTAGAATGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGtt                                                                                                                                            >  2:922900/1‑150 (MQ=255)
             gggTAGAATGGTTCCCTATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTTGATGGGGACATCCAGtttt                                                                                                                                          <  1:1178071/150‑1 (MQ=255)
               gTAGAATGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACTTAAATGACCGCGACCCGTTATTGTGTTCGTGGAAAAGAAAAACTTGATGGGGACATCCAGtttttc                                                                                                                                        <  2:922693/150‑2 (MQ=255)
                tAGAATGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTTGATGGGGACATCCAGtttttct                                                                                                                                       <  1:1099395/150‑3 (MQ=255)
                 aGAATGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATg                                                                                                                                      <  1:922900/150‑1 (MQ=255)
                 aGAATGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGAAGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTCTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTATTTATg                                                                                                                                      >  2:278378/1‑150 (MQ=255)
                  gAATGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATgg                                                                                                                                     >  1:979410/1‑150 (MQ=255)
                  gAATGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATgg                                                                                                                                     >  1:145877/1‑150 (MQ=255)
                  gAATGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATgg                                                                                                                                     <  1:601932/150‑1 (MQ=255)
                   aaTGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGt                                                                                                                                    >  2:35672/1‑150 (MQ=255)
                     tGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGc                                                                                                                                  <  2:627005/150‑1 (MQ=255)
                           ccGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGcc                                                                                                                            <  2:145877/150‑1 (MQ=255)
                             gaagCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAg                                                                                                                          <  1:511775/147‑1 (MQ=255)
                               tGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTa                                                                                                                        <  1:819566/150‑1 (MQ=255)
                                           cTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGatca                                                                                                            >  1:1113907/1‑150 (MQ=255)
                                            tGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGatcaa                                                                                                           >  1:313645/1‑150 (MQ=255)
                                             gTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGatcaat                                                                                                          >  1:461332/1‑150 (MQ=255)
                                             gTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGatcaat                                                                                                          <  2:1428057/150‑1 (MQ=255)
                                               tAAAACGCATATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGAACGAGACACGTGATGGTGTTAGTGGAAAAGAAAAAATGGATTTCACCATCAAGTTTTTATGGTGCTACGCCAGTAAAAAATTGataaatca                                                                                                        >  1:1410594/1‑150 (MQ=255)
                                                aaaaCGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGAtcaatcaa                                                                                                       >  1:1066095/1‑150 (MQ=255)
                                                 aaaCGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAAc                                                                                                      >  1:1000052/1‑150 (MQ=255)
                                                 aaaCGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAAc                                                                                                      >  1:1550653/1‑150 (MQ=255)
                                                  aaCGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACa                                                                                                     >  1:798098/1‑150 (MQ=255)
                                                                  tttACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATcac                                                                                     >  2:826471/1‑150 (MQ=255)
                                                                    tACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACatc                                                                                           >  2:1149949/1‑142 (MQ=255)
                                                                    tACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACatc                                                                                           <  1:1149949/142‑1 (MQ=255)
                                                                    tACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCacac                                                                                   >  2:1526925/1‑150 (MQ=255)
                                                                     aCATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACAcc                                                                                  <  2:700620/150‑1 (MQ=255)
                                                                      cATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACAccc                                                                                 >  1:1657344/1‑150 (MQ=255)
                                                                       aTGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCa                                                                                <  1:35672/150‑1 (MQ=255)
                                                                         ggACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATg                                                                              >  2:747671/1‑150 (MQ=255)
                                                                             aGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACTTCAATCACACCCATGAGGt                                                                          <  2:486076/150‑1 (MQ=255)
                                                                                tGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCt                                                                       <  1:1476816/150‑1 (MQ=255)
                                                                                  ccATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTc                                                                     >  2:436421/1‑150 (MQ=255)
                                                                                    aTCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAg                                                                   <  2:355017/150‑1 (MQ=255)
                                                                                     tCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAgg                                                                  <  1:435224/150‑1 (MQ=255)
                                                                                     tCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAgg                                                                  <  2:1000052/150‑1 (MQ=255)
                                                                                       gagcATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGAt                                                                <  2:382274/147‑1 (MQ=255)
                                                                                        aCCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATg                                                               <  2:313645/150‑1 (MQ=255)
                                                                                        aCCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATg                                                               >  1:1010445/1‑150 (MQ=255)
                                                                                        aCCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATg                                                               >  1:63433/1‑150 (MQ=255)
                                                                                            tGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAg                                                           <  2:869573/150‑1 (MQ=255)
                                                                                              gTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTTCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCa                                                         >  2:905999/1‑150 (MQ=255)
                                                                                               tACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCaa                                                        >  2:947333/1‑150 (MQ=255)
                                                                                               tACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCaa                                                        >  1:812728/1‑150 (MQ=255)
                                                                                                aCGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAg                                                       >  2:1056081/1‑150 (MQ=255)
                                                                                                aCGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAg                                                       >  2:1073189/1‑150 (MQ=255)
                                                                                                aCGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAg                                                       <  1:1480010/150‑1 (MQ=255)
                                                                                                  ggAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTg                                                     <  2:893675/150‑1 (MQ=255)
                                                                                                   gAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGc                                                    >  2:1227239/1‑150 (MQ=255)
                                                                                                         aCCGCGACGCGTGATGGTGTTCGTGGAAAAGAATAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAgatg                                              <  1:278378/150‑1 (MQ=255)
                                                                                                         aCCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAgatg                                              <  2:1432375/150‑1 (MQ=255)
                                                                                                              gACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGAttgt                                         >  1:781299/1‑150 (MQ=255)
                                                                                                                cccGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGAttgttg                                       >  1:142705/1‑150 (MQ=255)
                                                                                                                 ccGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTgg                                      <  1:504095/150‑1 (MQ=255)
                                                                                                                  cGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGc                                     >  1:566433/1‑150 (MQ=255)
                                                                                                                  cGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGc                                     <  1:1400189/150‑1 (MQ=255)
                                                                                                                   gTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCa                                    >  2:921407/1‑150 (MQ=255)
                                                                                                                    tGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCAt                                   >  1:1005319/1‑150 (MQ=255)
                                                                                                                    tGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCAt                                   >  1:399060/1‑150 (MQ=255)
                                                                                                                    tGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCAt                                   <  2:1066095/150‑1 (MQ=255)
                                                                                                                       tGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGa                                <  2:979410/150‑1 (MQ=255)
                                                                                                                              cGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTc                         >  2:1305317/1‑150 (MQ=255)
                                                                                                                              cGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTc                         <  1:509908/150‑1 (MQ=255)
                                                                                                                               ggtgaaaagaaaaACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCa                        <  1:1352726/147‑1 (MQ=255)
                                                                                                                                 ggAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATg                      >  1:59487/1‑150 (MQ=255)
                                                                                                                                   ataagaataACTGGTTGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTTATCAATCTACTTGGCGACATCAATCACACCGATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGAc                    <  1:1612413/148‑1 (MQ=255)
                                                                                                                                   accataaaaaCTTGATGGGGACATCCAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGAc                    >  2:387783/6‑150 (MQ=255)
                                                                                                                                   aaaagaaaaACTGGATGTCCCCATAAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGAc                    >  2:1491928/1‑150 (MQ=255)
                                                                                                                                    ccataaaaaCTTGATGGGGACATCCAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACg                   >  1:1309108/5‑150 (MQ=255)
                                                                                                                                      ataaaaaCTTGATGGGGACATCCAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGc                 >  1:1066427/3‑150 (MQ=255)
                                                                                                                                      agaaaaACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGc                 >  1:799272/1‑150 (MQ=255)
                                                                                                                                      agaaaaACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGc                 >  2:905525/1‑150 (MQ=255)
                                                                                                                                       gaaaaACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGcc                >  2:912918/1‑150 (MQ=255)
                                                                                                                                           aaCTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCgg            <  1:1227239/150‑1 (MQ=255)
                                                                                                                                            aCTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGc           <  2:799272/150‑1 (MQ=255)
                                                                                                                                              tGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTa         <  2:1550653/150‑1 (MQ=255)
                                                                                                                                                   gTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTAt    >  2:652143/1‑150 (MQ=255)
                                                                                                                                                    tCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTAtt   >  2:1476672/1‑150 (MQ=255)
                                                                                                                                                    tCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTAtt   >  2:1482745/1‑150 (MQ=255)
                                                                                                                                                     ccccATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTc              <  1:668096/138‑1 (MQ=255)
                                                                                                                                                     ccccATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTc              >  2:668096/1‑138 (MQ=255)
                                                                                                                                                     ccccATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTg  >  2:493042/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGCTGGGATCTGGGGGTAGAATGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTG  >  NC_000913/1661792‑1662090

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