Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 894,796 C→A 5.5% P338H (CCT→CAT) ‡ potF → putrescine ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913894,7960CA 5.5% 311.5 / 5.6 91P338H (CCT→CAT) ‡potFputrescine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (48/38);  new base A (3/2);  total (51/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

CAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAAT  >  NC_000913/894647‑894942
                                                                                                                                                     |                                                                                                                                                  
cAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGcc                                                                                                                                                    >  1:900087/1‑150 (MQ=255)
cAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGcc                                                                                                                                                    >  2:1416995/1‑150 (MQ=255)
cAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGcc                                                                                                                                                    >  2:826093/1‑150 (MQ=255)
   aGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGc                                                                                                                                                 >  2:1239922/1‑150 (MQ=255)
   aGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGc                                                                                                                                                 <  1:37935/150‑1 (MQ=255)
      cGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGa                                                                                                                                              >  1:654414/1‑150 (MQ=255)
        aaGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCAGGCGGAta                                                                                                                                            <  1:417069/150‑2 (MQ=255)
         aGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGt                                                                                                                                           <  2:911657/150‑1 (MQ=255)
            cTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCg                                                                                                                                        >  1:1261306/1‑150 (MQ=255)
               tCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGc                                                                                                                                     <  2:905324/150‑1 (MQ=255)
               tCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGc                                                                                                                                     <  2:1541795/150‑1 (MQ=255)
                cAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCg                                                                                                                                    <  2:45812/150‑1 (MQ=255)
                  gTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGaa                                                                                                                                  >  2:989485/1‑150 (MQ=255)
                    tCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGc                                                                                                                                >  2:494212/1‑150 (MQ=255)
                      cTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTg                                                                                                                              >  1:1148153/1‑150 (MQ=255)
                          aTTACCTGCTTATCCCGTATGTATTAGCGCATATTTCCGACAATGTGTTCTATGCAAAATCCAATAAAGCAGCCACGACTCTGGTGAGTGCGGAAGTCCGTTAGAACACAGGTAATTAACCGCCTGCGGATGTTCGTGAGAAGCTGtaaa                                                                                                                          >  2:630486/1‑147 (MQ=255)
                          aTTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCa                                                                                                                          >  2:21949/1‑150 (MQ=255)
                          aTTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCa                                                                                                                          >  2:197718/1‑150 (MQ=255)
                          aTTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCa                                                                                                                          >  1:267808/1‑150 (MQ=255)
                           ttACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCAc                                                                                                                         <  1:826093/150‑1 (MQ=255)
                             aCCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCActc                                                                                                                       >  2:914481/1‑150 (MQ=255)
                             aCCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCActc                                                                                                                       <  1:782483/150‑1 (MQ=255)
                                    gcgcCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGt                                                                                                                <  1:128574/150‑1 (MQ=255)
                                       cccGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACACCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCa                                                                                                             >  2:93759/1‑150 (MQ=255)
                                        ccGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAg                                                                                                            <  1:494212/150‑1 (MQ=255)
                                        ccGGATCTAGTAGCGCATATTTCCTATCATTTGTTCTATGCCAACTCCAATAAATCGTCGACTCCGCTGGTGAGTGCTGAAGTCCGTTATAACGAAGTTATTTATATGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAg                                                                                                            <  2:1008651/150‑1 (MQ=255)
                                             tGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATcc                                                                                                       >  2:1010352/1‑150 (MQ=255)
                                               tagtagCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGa                                                                                                     >  2:1633601/1‑150 (MQ=255)
                                                agtagCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGaa                                                                                                    <  1:329916/150‑1 (MQ=255)
                                                   agCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGAATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAAt                                                                                                 >  1:1352302/1‑150 (MQ=255)
                                                        atatTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGAcc                                                                                            <  1:197718/150‑1 (MQ=255)
                                                           tttCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCgtg                                                                                         <  1:434304/150‑1 (MQ=255)
                                                              ccGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGa                                                                                                     >  2:924176/1‑135 (MQ=255)
                                                              ccGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGa                                                                                                     <  1:924176/135‑1 (MQ=255)
                                                              ccGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTgc                                                                                      >  1:1814008/1‑150 (MQ=255)
                                                                      gtgtTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCcgc                                                                              >  2:1102744/1‑150 (MQ=255)
                                                                       tgtTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCcgcg                                                                             <  1:1416995/150‑1 (MQ=255)
                                                                           cTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTg                                                                         >  2:1198704/1‑150 (MQ=255)
                                                                           cTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTg                                                                         <  2:543582/150‑1 (MQ=255)
                                                                            tATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTgg                                                                        <  2:1148153/150‑1 (MQ=255)
                                                                             aTGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGa                                                                       >  2:1035851/1‑150 (MQ=255)
                                                                             aTGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGa                                                                       >  2:1029008/1‑150 (MQ=255)
                                                                              tGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGAc                                                                      >  1:28581/1‑150 (MQ=255)
                                                                              tGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGATCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGAc                                                                      <  2:654414/150‑1 (MQ=255)
                                                                               gCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGAcc                                                                     <  2:1261306/150‑1 (MQ=255)
                                                                                           aaaGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGc                                                         >  2:496660/1‑150 (MQ=255)
                                                                                               cagcCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGaa                                                     <  1:1633601/150‑1 (MQ=255)
                                                                                                  ccACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAaat                                                  >  2:1975187/1‑150 (MQ=255)
                                                                                                  ccACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAaat                                                  >  2:1441071/1‑150 (MQ=255)
                                                                                                      gccgcTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATc                                              <  1:989485/150‑1 (MQ=255)
                                                                                                            ggTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAg                                         >  1:77621/1‑149 (MQ=255)
                                                                                                            ggTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAg                                         <  2:77621/149‑1 (MQ=255)
                                                                                                            ggTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGt                                        <  1:1463818/150‑1 (MQ=255)
                                                                                                              tGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCg                                      <  1:554526/150‑1 (MQ=255)
                                                                                                              tGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCg                                      <  2:1700970/150‑1 (MQ=255)
                                                                                                                tgTGCGGAAGTCCGTGAGAACCCAGTTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTa                                    <  2:1945433/149‑1 (MQ=255)
                                                                                                                aGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTa                                    >  1:587442/1‑150 (MQ=255)
                                                                                                                 gTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAg                                   <  2:28581/150‑1 (MQ=255)
                                                                                                                 gTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAg                                   >  2:641645/1‑150 (MQ=255)
                                                                                                                 gTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAg                                   >  1:1959224/1‑150 (MQ=255)
                                                                                                                  tGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGa                                  >  2:1656581/1‑150 (MQ=255)
                                                                                                                    cGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATg                                <  1:1102744/150‑1 (MQ=255)
                                                                                                                      gAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGcc                              >  2:143757/1‑150 (MQ=255)
                                                                                                                        aGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCgg                            <  2:1106583/150‑1 (MQ=255)
                                                                                                                        aGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCgg                            >  1:21592/1‑150 (MQ=255)
                                                                                                                          tCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAg                          >  2:717693/1‑150 (MQ=255)
                                                                                                                            cGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAggg                        <  2:1814008/150‑1 (MQ=255)
                                                                                                                            cGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAggg                        >  2:1760885/1‑150 (MQ=255)
                                                                                                                              tGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGgc                      >  2:86594/1‑150 (MQ=255)
                                                                                                                              tGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGgc                      >  2:235935/1‑150 (MQ=255)
                                                                                                                              tGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGgc                      >  2:1079865/1‑150 (MQ=255)
                                                                                                                               gagaACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCATGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGgcg                     >  2:1925849/1‑150 (MQ=255)
                                                                                                                                agaACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGgggc                    <  1:1239922/150‑1 (MQ=255)
                                                                                                                                agaACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGgcgc                    >  2:694514/1‑150 (MQ=255)
                                                                                                                                   aCCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGcacc                 <  2:1065624/150‑1 (MQ=255)
                                                                                                                                   aCCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGcacc                 >  1:1649703/1‑150 (MQ=255)
                                                                                                                                    cccAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGcacca                >  2:1483342/1‑150 (MQ=255)
                                                                                                                                      cAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGcaccacc              >  2:898005/1‑149 (MQ=255)
                                                                                                                                      cAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCCGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACcaca              <  1:1312248/150‑1 (MQ=255)
                                                                                                                                        cacgaacaTCCGCAGGCGGATGTTCGTGCGAAGCTGTTCCCTCTGAAAGTGCAGGATCCGAAAATCGCCCGTGTGCGCCCCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCCCCACAcc            >  1:757918/8‑150 (MQ=255)
                                                                                                                                        cacgaacaTCCGCAGGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGcac                  >  1:1335472/8‑144 (MQ=255)
                                                                                                                                        cacgaacaTCCGCAGGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGcac                  <  2:1335472/137‑1 (MQ=255)
                                                                                                                                        cacgaacaTCCGCAGGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACAcc            >  1:988137/8‑150 (MQ=255)
                                                                                                                                            tttATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACccgcc        >  2:1229770/1‑150 (MQ=255)
                                                                                                                                              tATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCgg      <  1:309096/150‑1 (MQ=255)
                                                                                                                                              tATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCgg      >  2:241541/1‑150 (MQ=255)
                                                                                                                                              tATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCgg      >  1:1670834/1‑150 (MQ=255)
                                                                                                                                               aTCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGc     <  1:1035851/150‑1 (MQ=255)
                                                                                                                                               aTCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGc     <  1:1029008/150‑1 (MQ=255)
                                                                                                                                               aTCCGCCTGCGAATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGc     <  2:1352302/150‑1 (MQ=255)
                                                                                                                                                tCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCa    >  1:920443/1‑150 (MQ=255)
                                                                                                                                                tCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCa    >  2:617252/1‑150 (MQ=255)
                                                                                                                                                   gCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAAt  >  2:1236942/1‑149 (MQ=255)
                                                                                                                                                   gCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAAt  <  1:1236942/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAAT  >  NC_000913/894647‑894942

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: