Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,940,283 G→C 8.3% intergenic (‑90/+30) lpxM ← / ← mepM myristoyl‑acyl carrier protein‑dependent acyltransferase/peptidoglycan DD‑endopeptidase MepM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,940,2830GC 8.3% 258.1 / 10.7 81intergenic (‑90/+30)lpxM/mepMmyristoyl‑acyl carrier protein‑dependent acyltransferase/peptidoglycan DD‑endopeptidase MepM
Reads supporting (aligned to +/- strand):  ref base G (33/41);  new base C (2/5);  total (35/46)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCT  >  NC_000913/1940135‑1940431
                                                                                                                                                    |                                                                                                                                                    
tGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTactga                                                                                                                                                     <  1:898696/150‑1 (MQ=255)
 ggCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTactgac                                                                                                                                                    <  1:930213/150‑1 (MQ=255)
 ggCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTactgaa                                                                                                                                                    >  2:557379/1‑149 (MQ=255)
  gCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTactgact                                                                                                                                                   >  1:697375/1‑150 (MQ=255)
   cGAAAGGATTTATCAAACTCAGGAATGTATTCGCTTTTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTactgactg                                                                                                                                                  <  1:418983/150‑1 (MQ=255)
   cGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCATTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTactgactg                                                                                                                                                  >  2:284140/1‑150 (MQ=255)
    taaaGGATTTATCAAACTCAGGAATGTATTCTCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGc                                                                                                                                                 <  1:693762/149‑1 (MQ=255)
     aaaGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCg                                                                                                                                                <  2:777203/150‑1 (MQ=255)
     aaaGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCg                                                                                                                                                <  1:22373/150‑1 (MQ=255)
      aaGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGt                                                                                                                                               >  2:1630726/1‑150 (MQ=255)
         gATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTAcc                                                                                                                                            <  1:947983/150‑1 (MQ=255)
          aTTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCg                                                                                                                                           <  1:1268338/150‑1 (MQ=255)
             tATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCt                                                                                                                                        <  2:557041/150‑1 (MQ=255)
             tATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCt                                                                                                                                        >  2:1392706/1‑150 (MQ=255)
             tATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACGCAGTCAGTACCGGCt                                                                                                                                        <  1:1544299/150‑1 (MQ=255)
             tATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACCGACTGCGTACCGGCt                                                                                                                                        >  1:460529/1‑150 (MQ=255)
              aTCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTg                                                                                                                                       >  2:1485596/1‑150 (MQ=255)
              aTCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTg                                                                                                                                       >  1:649769/1‑150 (MQ=255)
               tCAAACTCAGGAATGTAGTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATACGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACGCAGTCAGTACCGGCttt                                                                                                                                      <  2:1028760/150‑3 (MQ=255)
                 aaaCTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACGCAGTCAGTACCGGCttttt                                                                                                                                    <  1:1100155/150‑5 (MQ=255)
                 aaaCTCAGGAATATATTCTCTATTATTTTTTTTCGTTTCAATGCTTTTCCAGTTTCTGATAAGGCAAAAATAAATATTGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGa                                                                                                                                    >  2:676087/1‑150 (MQ=255)
                   aCTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACGCAGTCAGTACCGGCttttttt                                                                                                                                  <  2:1324158/150‑7 (MQ=255)
                    cTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATg                                                                                                                                 <  1:1191813/150‑1 (MQ=255)
                      cAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGa                                                                                                                               >  1:806946/1‑150 (MQ=255)
                         gAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGt                                                                                                                            <  1:1188347/150‑1 (MQ=255)
                            tGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGttaa                                                                                                                         <  2:1312174/150‑1 (MQ=255)
                               aTTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGttaatta                                                                                                                      >  2:1619390/1‑150 (MQ=255)
                                ttCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGttaattaa                                                                                                                     <  2:396238/150‑1 (MQ=255)
                                 tCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTtaattaat                                                                                                                    <  2:567127/150‑1 (MQ=255)
                                  cGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACGCAGTCAGTACCGGCttttttt                                                                                                                                  >  2:927325/1‑129 (MQ=255)
                                  cGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACGCAGTCAGTACCGGCttttttt                                                                                                                                  <  1:927325/135‑7 (MQ=255)
                                    cTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCaa                                                                                                                 >  1:980050/1‑150 (MQ=255)
                                    cTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGAGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCaa                                                                                                                 >  2:32307/1‑150 (MQ=255)
                                       ttattTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAAcc                                                                                                              <  1:1561357/150‑1 (MQ=255)
                                         attTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATTGATGTTAATTAATCAAACCGt                                                                                                            <  2:82006/150‑1 (MQ=255)
                                          ttttttttCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTa                                                                                                           <  1:1562679/150‑1 (MQ=255)
                                           tttttttCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAg                                                                                                          >  1:53649/1‑150 (MQ=255)
                                           tttttttCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCATACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAg                                                                                                          <  1:568177/150‑1 (MQ=255)
                                                 tCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCgg                                                                                                    <  1:1382765/150‑1 (MQ=255)
                                                  cGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGc                                                                                                   <  1:444296/150‑1 (MQ=255)
                                                     ttCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGcaca                                                                                                <  2:132039/150‑1 (MQ=255)
                                                       ccATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAACCCGTAGCTGCGGCACAAt                                                                                              >  2:1103441/1‑150 (MQ=255)
                                                       ccATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAAt                                                                                              >  1:1262065/1‑150 (MQ=255)
                                                       caaTGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGAGCCGCAACTTGCACCGCACCAAATAAAAAAAGCCGGTACTGACTCCGTACCGGATGCGAATGGATCTTAATTAATCAAACCGTAGCTCCGGCACAAt                                                                                              >  2:1223173/3‑150 (MQ=255)
                                                        cATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAAtc                                                                                             <  2:408441/150‑1 (MQ=255)
                                                        cATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAAtc                                                                                             <  2:1083772/150‑1 (MQ=255)
                                                            cTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCtt                                                                                         <  1:686014/150‑1 (MQ=255)
                                                            cTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCtt                                                                                         >  1:121051/1‑150 (MQ=255)
                                                             ttttCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCttt                                                                                        <  1:1392706/150‑1 (MQ=255)
                                                             ttttCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCttt                                                                                        <  2:754634/150‑1 (MQ=255)
                                                             ttttCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCttt                                                                                        <  2:1235709/150‑1 (MQ=255)
                                                              tttCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTg                                                                                       <  1:449143/150‑1 (MQ=255)
                                                               ttCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGGTGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTgg                                                                                      >  2:200994/1‑150 (MQ=255)
                                                                     tttCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTgcg                                                                                <  2:215417/150‑1 (MQ=255)
                                                                     tttCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTgcg                                                                                <  1:14296/150‑1 (MQ=255)
                                                                     tttCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTgcg                                                                                <  1:1177299/150‑1 (MQ=255)
                                                                             aaGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAAt                                                                        >  1:692676/1‑150 (MQ=255)
                                                                               ggCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTc                                                                      >  1:1386266/1‑150 (MQ=255)
                                                                               ggCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTc                                                                      <  2:1632529/150‑1 (MQ=255)
                                                                               ggCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTc                                                                      >  2:1413745/1‑150 (MQ=255)
                                                                                 cAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTcgc                                                                    <  1:437417/150‑1 (MQ=255)
                                                                                         aatcTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCg                                                            >  2:1557797/1‑150 (MQ=255)
                                                                                          atcTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCgg                                                           >  2:536480/1‑150 (MQ=255)
                                                                                            cTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAg                                                         <  1:1630726/150‑1 (MQ=255)
                                                                                            cTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAg                                                         >  2:1219407/1‑150 (MQ=255)
                                                                                              ggTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGcc                                                       >  1:1211186/1‑150 (MQ=255)
                                                                                                    aGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAg                                                 <  1:1413745/150‑1 (MQ=255)
                                                                                                        tAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCt                                             >  2:111949/1‑150 (MQ=255)
                                                                                                        gtgcggggCAAGTTGCGCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGATGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCt                                             >  2:1195186/3‑150 (MQ=255)
                                                                                                         tgcggggCAAGTTGCGCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCtt                                            >  2:1268611/2‑150 (MQ=255)
                                                                                                         aGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCtt                                            <  2:1386266/150‑1 (MQ=255)
                                                                                                           cggcgCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCg                                          >  1:35203/1‑150 (MQ=255)
                                                                                                              cgcAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGAGTACCGGATGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTa                                       >  1:290280/1‑150 (MQ=255)
                                                                                                               gcAACTTGCCCCGCACCAAATAAAAAAAGCCGGAACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGAACAATCTCTATGGCCTGCCCCAGGAATTCGCGACGATCGGAGCCGGGCAGCCCTTAGGTAc                                      >  2:71851/1‑150 (MQ=255)
                                                                                                                 aaCTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGTATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTAcgc                                    <  2:1567159/150‑1 (MQ=255)
                                                                                                                  aCTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTAcgcg                                   <  2:1182535/150‑1 (MQ=255)
                                                                                                                    ttGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGc                                 <  2:160963/150‑1 (MQ=255)
                                                                                                                        cccGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGtt                             >  2:1142874/1‑150 (MQ=255)
                                                                                                                         ccGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGttt                            <  1:1485596/150‑1 (MQ=255)
                                                                                                                           gCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTg                          <  1:557379/150‑1 (MQ=255)
                                                                                                                              ccAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCg                       >  2:113986/1‑150 (MQ=255)
                                                                                                                                aaaTAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTc                     >  1:301698/1‑150 (MQ=255)
                                                                                                                                         caGCCGGTACGCAGTCAGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGttt            >  1:900892/2‑150 (MQ=255)
                                                                                                                                          aGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTa           >  2:11132/1‑150 (MQ=255)
                                                                                                                                            ccGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACg         <  2:485711/150‑1 (MQ=255)
                                                                                                                                            ccGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACg         <  2:200708/150‑1 (MQ=255)
                                                                                                                                               gTACTGACTGCGTACGGGATTCGAATGTATGTTAATTATTCAAACTGTAGCTGCTGCACAATCTCTTTGGCCTGCTCCAGGAATTCGCGACGATGGTAGCCGTTGAGCCGTTCGTTAGGCGTCAGTTTTGCGGTCAGCGGGTTTACGGcc      <  2:290280/150‑1 (MQ=255)
                                                                                                                                                  ctgactgCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCctgc   >  2:1311744/1‑150 (MQ=255)
                                                                                                                                                   tgactgCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCctgct  >  1:344771/1‑150 (MQ=255)
                                                                                                                                                    gactgCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGtt             >  1:926214/1‑138 (MQ=255)
                                                                                                                                                    gactgCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGtt             <  2:926214/138‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCT  >  NC_000913/1940135‑1940431

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: