Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,225,009 A→G 6.3% intergenic (‑139/+35) pbpG ← / ← yohC peptidoglycan DD‑endopeptidase PbpG/putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,225,0090AG 6.3% 283.6 / 7.4 95intergenic (‑139/+35)pbpG/yohCpeptidoglycan DD‑endopeptidase PbpG/putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (46/43);  new base G (2/4);  total (48/47)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.72e-01

CATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTT  >  NC_000913/2224868‑2225157
                                                                                                                                             |                                                                                                                                                    
cATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACg                                                                                                                                              <  1:683925/150‑1 (MQ=255)
  tgatgaGCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGcc                                                                                                                                            <  1:43556/150‑1 (MQ=255)
        gCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAAt                                                                                                                                      >  1:639709/1‑150 (MQ=255)
         cATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGGCGTAAATGTCAGCAATg                                                                                                                                     <  1:1100966/150‑1 (MQ=255)
         cATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGGCGTAAATGTCAGCAATg                                                                                                                                     >  2:1150877/1‑150 (MQ=255)
         cATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATg                                                                                                                                     >  1:1161428/1‑150 (MQ=255)
         cATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATg                                                                                                                                     <  2:747795/150‑1 (MQ=255)
          aTTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGGCGTAAATGTCAGCAATGc                                                                                                                                    <  2:941993/150‑1 (MQ=255)
           ttCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCa                                                                                                                                   >  2:1360988/1‑150 (MQ=255)
            tCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGGCGTAAATGTCAGCAATGCag                                                                                                                                  <  2:93093/150‑2 (MQ=255)
             cAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGGCGTAAATGTCAGCAATGCagc                                                                                                                                 <  1:50381/150‑3 (MQ=255)
             cAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAAc                                                                                                                                 <  1:1010490/150‑1 (MQ=255)
             cAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAAc                                                                                                                                 <  1:735968/150‑1 (MQ=255)
              aGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACg                                                                                                                                <  2:771725/150‑1 (MQ=255)
              aGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGAATTGCTGACATTTACGCCAGCAATGCAACg                                                                                                                                >  1:533388/1‑150 (MQ=255)
               gATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGt                                                                                                                               <  1:619175/150‑1 (MQ=255)
                aTAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTc                                                                                                                              <  2:1415870/150‑1 (MQ=255)
                   gTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCaaa                                                                                                                           <  1:183765/150‑1 (MQ=255)
                   gTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCaaa                                                                                                                           <  1:189365/150‑1 (MQ=255)
                    tGGTTCACGCTCACGGGTTGCGCACCGCCGTATTAAGGATTTACTGATGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAg                                                                                                                          <  2:1350866/150‑1 (MQ=255)
                    tGGTTCACGCGCACTGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAg                                                                                                                          <  1:1211740/150‑1 (MQ=255)
                    tGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAg                                                                                                                          >  1:533261/1‑150 (MQ=255)
                     ggTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGa                                                                                                                         >  1:273330/1‑150 (MQ=255)
                     ggTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGa                                                                                                                         >  1:457674/1‑150 (MQ=255)
                     ggTTCAAGCGCACTGGTGGCGGACCGCCGGAGTTAGGATTTACTGATGCTAGCGAAGCCATCATAACGAGGAAAATGTGCGTTCTTCATAGGAGAATCGTGAGGTAATGCTGCATTGCTGACATTTTCGCCAGCAATGCAACGTaaaata                                                                                                                         <  2:1406955/150‑3 (MQ=255)
                      gTTCACGCGGACGGGTTGCGCACCGGCGTAGTAATGATTTACTGAGGCTATCGTCGGCATCATAACGATCAAAAACTTCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGaa                                                                                                                        <  2:533388/150‑1 (MQ=255)
                       ttCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAAc                                                                                                                       <  1:1199421/150‑1 (MQ=255)
                       ttCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAAc                                                                                                                       >  2:245914/1‑150 (MQ=255)
                         cACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACtt                                                                                                                     >  2:316303/1‑150 (MQ=255)
                           cgcgCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTc                                                                                                                   >  1:933292/1‑150 (MQ=255)
                           cgcgCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTc                                                                                                                   <  2:895187/150‑1 (MQ=255)
                            gcgcACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCt                                                                                                                  <  2:156339/150‑1 (MQ=255)
                            gcgcACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCt                                                                                                                  <  1:461366/150‑1 (MQ=255)
                            gcgcACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCt                                                                                                                  <  1:1483954/150‑1 (MQ=255)
                               cACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGa                                                                                                               <  1:802292/150‑1 (MQ=255)
                                 cGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAAc                                                                                                             >  2:1194785/1‑150 (MQ=255)
                                 cGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAAc                                                                                                             <  2:719798/150‑1 (MQ=255)
                                   ggTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACaa                                                                                                           >  1:563057/1‑150 (MQ=255)
                                   ggTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACaa                                                                                                           <  1:1482221/150‑1 (MQ=255)
                                    gTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACaaa                                                                                                          >  1:894694/1‑150 (MQ=255)
                                      tGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACg                                                                                                        >  1:1292343/1‑150 (MQ=255)
                                        cgcACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGGAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGAt                                                                                                      >  2:448088/1‑150 (MQ=255)
                                        cgcACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGAt                                                                                                      >  1:483399/1‑150 (MQ=255)
                                        cgcACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACATTCTAGAACAAACGAt                                                                                                      >  1:1153039/1‑150 (MQ=255)
                                             cgccgGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCg                                                                                                 <  2:483399/150‑1 (MQ=255)
                                              gccgGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGt                                                                                                >  1:560765/1‑150 (MQ=255)
                                              gccgGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGt                                                                                                >  1:1226998/1‑150 (MQ=255)
                                                   aGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGccc                                                                                           <  1:331279/150‑1 (MQ=255)
                                                       aGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATa                                                                                       <  1:1360988/150‑1 (MQ=255)
                                                       aGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATa                                                                                       >  1:1050839/1‑150 (MQ=255)
                                                        ggATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAg                                                                                      >  1:930878/1‑150 (MQ=255)
                                                            ttACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGtata                                                                                  <  2:457674/150‑1 (MQ=255)
                                                             tACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTGTAGAACAAACGATAACCGGAGCCCCATAGTATaa                                                                                 >  1:647347/1‑150 (MQ=255)
                                                             tACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATaa                                                                                 >  1:1035683/1‑150 (MQ=255)
                                                             tACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATaa                                                                                 >  1:1469505/1‑150 (MQ=255)
                                                              aCTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAAc                                                                                <  1:215052/150‑1 (MQ=255)
                                                                tGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACgg                                                                              <  2:508863/150‑1 (MQ=255)
                                                                tGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACgg                                                                              <  2:1469505/150‑1 (MQ=255)
                                                                tGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACgg                                                                              >  1:1178037/1‑150 (MQ=255)
                                                                 gAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGt                                                                             >  2:1270119/1‑150 (MQ=255)
                                                                   ggCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTa                                                                           >  1:370156/1‑150 (MQ=255)
                                                                     cTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAgg                                                                         <  2:563057/150‑1 (MQ=255)
                                                                       aGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGc                                                                       <  1:1281498/150‑1 (MQ=255)
                                                                          gACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGag                                                                    >  1:194317/1‑150 (MQ=255)
                                                                            cGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGagca                                                                  >  1:651684/1‑150 (MQ=255)
                                                                               catcatAACGAGCAAAAAGTGCGATCGTCAAACGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAgcagca                                                               <  2:346283/150‑1 (MQ=255)
                                                                                    tAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGGAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTcc                                                          <  1:448088/150‑1 (MQ=255)
                                                                                      aCGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAg                                                        <  1:179430/150‑1 (MQ=255)
                                                                                         aGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTAc                                                     >  1:120319/1‑150 (MQ=255)
                                                                                          gCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTAcc                                                    <  2:1178037/150‑1 (MQ=255)
                                                                                          gCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTAcc                                                    <  2:1015380/150‑1 (MQ=255)
                                                                                          gCAAAAAGTGCGATCGGCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTAcc                                                    <  1:1498337/150‑1 (MQ=255)
                                                                                            aaaaaGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAg                                                  >  1:227309/1‑150 (MQ=255)
                                                                                                 gTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGc                                             <  1:513590/150‑1 (MQ=255)
                                                                                                 gTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGc                                             >  1:808071/1‑150 (MQ=255)
                                                                                                  tGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGcc                                            <  1:1270119/150‑1 (MQ=255)
                                                                                                     gATAGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAAAGATAACCGTAGCCCAATAGTATAACGGTTAGGGAGAGCAGCACTTCAAGTACCAGCACGCCAAt                                         >  2:257914/1‑150 (MQ=255)
                                                                                                      aTCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATc                                        >  1:928839/1‑150 (MQ=255)
                                                                                                          tCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCCCGCCAATCGCCa                                    >  2:256844/1‑150 (MQ=255)
                                                                                                          tCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCa                                    >  1:505673/1‑150 (MQ=255)
                                                                                                           cAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAg                                   >  2:268570/1‑150 (MQ=255)
                                                                                                           cAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTACCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAg                                   >  2:963921/1‑150 (MQ=255)
                                                                                                              aGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGt                                >  1:487659/1‑150 (MQ=255)
                                                                                                                gagaATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCg                              <  2:1008622/150‑1 (MQ=255)
                                                                                                                        tGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTagag                      >  2:74186/1‑150 (MQ=255)
                                                                                                                           ggAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAg                   >  2:725035/1‑150 (MQ=255)
                                                                                                                               aTGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTaa               <  1:316303/150‑1 (MQ=255)
                                                                                                                                   tgcATTGCTGGCGTAAATGTCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGcc           >  2:1409513/1‑150 (MQ=255)
                                                                                                                                    gcATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCt          >  1:609868/1‑150 (MQ=255)
                                                                                                                                       ttGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTcc       <  2:505673/150‑1 (MQ=255)
                                                                                                                                        tGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCct      <  1:466186/150‑1 (MQ=255)
                                                                                                                                         gCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCctc     >  2:658442/1‑150 (MQ=255)
                                                                                                                                         gCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAAGCCTTCctc     >  1:543795/1‑150 (MQ=255)
                                                                                                                                           tGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCtt   >  2:545079/1‑150 (MQ=255)
                                                                                                                                            gACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCttt  >  2:1251460/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
CATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTT  >  NC_000913/2224868‑2225157

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