Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,661,948 A→C 6.7% K465Q (AAG→CAG)  ynfF → putative selenate reductase YnfF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,661,9480AC 6.7% 273.0 / 8.6 87K465Q (AAG→CAG) ynfFputative selenate reductase YnfF
Reads supporting (aligned to +/- strand):  ref base A (45/36);  new base C (3/3);  total (48/39)
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.00e+00

ATCTGGGGGTAGAATGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATC  >  NC_000913/1661799‑1662097
                                                                                                                                                     |                                                                                                                                                     
aTCTGGGGGTAGAATGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCa                                                                                                                                                       <  1:371483/150‑1 (MQ=255)
    gggggTAGAATGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGtt                                                                                                                                                   >  2:922900/1‑150 (MQ=255)
    gggggTAGAATGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGtt                                                                                                                                                   >  2:1476816/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:601932/150‑1 (MQ=255)
           gAATGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATgg                                                                                                                                            >  1:145877/1‑150 (MQ=255)
           gAATGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATgg                                                                                                                                            >  1:979410/1‑150 (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:1550653/1‑150 (MQ=255)
                                          aaaCGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAAc                                                                                                             >  1:1000052/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:63433/1‑150 (MQ=255)
                                                                                 aCCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATg                                                                      >  1:1010445/1‑150 (MQ=255)
                                                                                     tGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAg                                                                  <  2:869573/150‑1 (MQ=255)
                                                                                       gTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTTCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCa                                                                >  2:905999/1‑150 (MQ=255)
                                                                                        tACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCaa                                                               >  1:812728/1‑150 (MQ=255)
                                                                                        tACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCaa                                                               >  2:947333/1‑150 (MQ=255)
                                                                                         aCGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAg                                                              >  2:1073189/1‑150 (MQ=255)
                                                                                         aCGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAg                                                              >  2:1056081/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:1400189/150‑1 (MQ=255)
                                                                                                           cGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGc                                            >  1:566433/1‑150 (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                                <  1:509908/150‑1 (MQ=255)
                                                                                                                       cGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTc                                >  2:1305317/1‑150 (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                        >  2:905525/1‑150 (MQ=255)
                                                                                                                               agaaaaACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGc                        >  1:799272/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)
                                                                                                                                               cccATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGc        <  2:1657344/150‑1 (MQ=255)
                                                                                                                                                ccATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCg       <  2:142705/150‑1 (MQ=255)
                                                                                                                                                    cAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGatat   <  1:1491928/150‑1 (MQ=255)
                                                                                                                                                     aaGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATc  <  1:22219/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATCTGGGGGTAGAATGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGTTTTTATGGTGCTACGCCAGTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATC  >  NC_000913/1661799‑1662097

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: