Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,525,563 T→A 5.5% intergenic (+140/+26) mrcA → / ← nudE peptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcA/ADP‑sugar diphosphatase NudE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,525,5630TA 5.5% 534.2 / 13.7 163intergenic (+140/+26)mrcA/nudEpeptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcA/ADP‑sugar diphosphatase NudE
Reads supporting (aligned to +/- strand):  ref base T (69/84);  new base A (3/6);  total (73/90)
Fisher's exact test for biased strand distribution p-value = 7.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 164 total aligned reads displayed.

TTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCCATCATATGCGCCAGCGGCCAGCGCACCT  >  NC_000913/3525415‑3525710
                                                                                                                                                    |                                                                                                                                                   
  gTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATtgt                                                                                                                                                  >  2:1656960/1‑150 (MQ=255)
   ttCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATtgtg                                                                                                                                                 <  1:1838063/150‑1 (MQ=255)
   ttCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACAGtgtg                                                                                                                                                 >  2:425976/1‑150 (MQ=255)
    tCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCTTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGc                                                                                                                                                >  1:2074978/1‑150 (MQ=255)
     caaactgaaaaGGCGCCGAAGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCa                                                                                                                                               >  2:1471474/8‑150 (MQ=255)
     cTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCa                                                                                                                                               <  1:1218463/150‑1 (MQ=255)
     cTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCa                                                                                                                                               >  1:1084098/1‑150 (MQ=255)
      tGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCAt                                                                                                                                              >  1:1088801/1‑150 (MQ=255)
       gATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATc                                                                                                                                             >  1:530074/1‑150 (MQ=255)
         ttAAAAAGGCGCTTCGGCGCCTTTTCAGTTTACTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAg                                                                                                                                           >  2:764301/1‑150 (MQ=255)
          tAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAgg                                                                                                                                          <  1:924431/150‑1 (MQ=255)
           aaaaaGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATGCACAATGTTTATCAGGc                                                                                                                                         <  1:810224/150‑1 (MQ=255)
            aaaaGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATGCACAATGTTTATCAGGct                                                                                                                                        <  2:451559/150‑2 (MQ=255)
            aaaaGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCa                                                                                                                                        >  2:1691320/1‑150 (MQ=255)
             aaaGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATGCACAATGTTTATCAGGctt                                                                                                                                       <  1:1475697/150‑3 (MQ=255)
                ggCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACt                                                                                                                                    >  1:437265/1‑150 (MQ=255)
                ggCGCTTAGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACt                                                                                                                                    >  2:1641147/1‑150 (MQ=255)
                    cTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCAc                                                                                                                                <  1:167492/150‑1 (MQ=255)
                    cTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTATGCATCAGGCAAACTTCAc                                                                                                                                <  2:282353/150‑1 (MQ=255)
                         gcgcCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCAtt                                                                                                                           <  2:1525436/150‑1 (MQ=255)
                         gcgcCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCAtt                                                                                                                           >  1:1625543/1‑150 (MQ=255)
                          cgcCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCAttt                                                                                                                          <  2:2093121/150‑1 (MQ=255)
                           gcCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTa                                                                                                                         <  2:1521161/150‑1 (MQ=255)
                            ccTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTac                                                                                                                        <  2:2372934/150‑1 (MQ=255)
                             cTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTaca                                                                                                                       <  2:1782382/150‑1 (MQ=255)
                               tttCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACt                                                                                                                     >  1:640801/1‑150 (MQ=255)
                                  cAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGc                                                                                                                  >  1:965069/1‑150 (MQ=255)
                                           gACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCt                                                                                                         <  2:2274579/150‑1 (MQ=255)
                                           gACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCt                                                                                                         >  1:2194438/1‑150 (MQ=255)
                                           gACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCt                                                                                                         >  1:1809817/1‑150 (MQ=255)
                                               aaGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCa                                                                                                     >  2:938253/1‑150 (MQ=255)
                                               aaGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTCAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCa                                                                                                     <  2:800344/150‑1 (MQ=255)
                                                 gTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAAc                                                                                                   <  1:162623/150‑1 (MQ=255)
                                                 gTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAAc                                                                                                   <  1:691495/150‑1 (MQ=255)
                                                 gTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAAc                                                                                                   <  2:1251694/150‑1 (MQ=255)
                                                    cACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCAt                                                                                                >  1:293406/1‑150 (MQ=255)
                                                    cACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCAt                                                                                                >  1:2188912/1‑150 (MQ=255)
                                                        tGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTcgc                                                                                            >  2:607735/1‑150 (MQ=255)
                                                        tGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTcgc                                                                                            <  1:592366/150‑1 (MQ=255)
                                                        tGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTcgc                                                                                            <  1:252250/150‑1 (MQ=255)
                                                          tttATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGACAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCgcgc                                                                                          <  2:311282/150‑1 (MQ=255)
                                                          tgttagCCGGATAACGCGTTCACGCCGTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCgcgc                                                                                          >  1:2424527/6‑150 (MQ=255)
                                                            tATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCAc                                                                                        <  2:1590647/150‑1 (MQ=255)
                                                              agCCGGATAACGCGTTCACGCCGTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGa                                                                                      >  2:1412498/2‑150 (MQ=255)
                                                               gCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAg                                                                                     >  1:519097/1‑150 (MQ=255)
                                                               gCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAg                                                                                     <  1:453715/150‑1 (MQ=255)
                                                               gCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAg                                                                                     >  2:1658937/1‑150 (MQ=255)
                                                                ccGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCGCTGCCCTTTCAACCATTCGCGCACGAgg                                                                                    >  2:2246052/1‑150 (MQ=255)
                                                                ccGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAgg                                                                                    >  1:98308/1‑150 (MQ=255)
                                                                ccGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCAATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCCAACGCGCACGAgg                                                                                    >  2:1201562/1‑150 (MQ=255)
                                                                   gATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAAc                                                                                 <  2:519916/150‑1 (MQ=255)
                                                                     tACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAg                                                                               <  2:186069/150‑1 (MQ=255)
                                                                        ggCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAgcgc                                                                            >  1:1895698/1‑150 (MQ=255)
                                                                        ggCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCACTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAgg                                                                                    <  1:1808957/142‑1 (MQ=255)
                                                                               aCGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACa                                                                     <  1:80272/150‑1 (MQ=255)
                                                                                  cgTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTg                                                                  <  1:1529926/150‑1 (MQ=255)
                                                                                   gTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTgc                                                                 >  1:1429757/1‑150 (MQ=255)
                                                                                   gTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTCACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTgc                                                                 <  1:1223622/150‑1 (MQ=255)
                                                                                     tATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTgcgc                                                               <  1:665149/150‑1 (MQ=255)
                                                                                      aTCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTgcgcg                                                              <  1:2155620/150‑1 (MQ=255)
                                                                                        ccGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCt                                                            <  1:430790/150‑1 (MQ=255)
                                                                                          ggCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTc                                                          <  2:139678/150‑1 (MQ=255)
                                                                                            cTAACAAAATCGTGAAAACTCAATAAATTGCAGAAATCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCAt                                                        <  2:1954659/150‑1 (MQ=255)
                                                                                             tAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCAtt                                                       <  2:2188912/150‑1 (MQ=255)
                                                                                              aaCAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTg                                                      <  1:1992632/150‑1 (MQ=255)
                                                                                               aCAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGa                                                     <  1:2422135/150‑1 (MQ=255)
                                                                                                cAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGaa                                                    >  1:54214/1‑150 (MQ=255)
                                                                                                 aaaaTCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAg                                                   <  1:633063/150‑1 (MQ=255)
                                                                                                     tCGTGAAAACTCAATAAATTGCAGAAACCCTATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAg                                               <  1:2094398/150‑1 (MQ=255)
                                                                                                       gTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAggg                                             >  1:963412/1‑150 (MQ=255)
                                                                                                       gTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCGTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCCggg                                             >  2:2145099/1‑150 (MQ=255)
                                                                                                       gTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAAAATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAggg                                             >  2:50810/1‑150 (MQ=255)
                                                                                                            aaCTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTc                                        <  2:1049626/150‑1 (MQ=255)
                                                                                                              cTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCATGGTCTTCCa                                      >  2:158352/1‑150 (MQ=255)
                                                                                                              cTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCa                                      <  2:848866/150‑1 (MQ=255)
                                                                                                                cAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCTAGc                                    <  2:1809817/150‑1 (MQ=255)
                                                                                                                cAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGc                                    <  1:1656960/150‑1 (MQ=255)
                                                                                                                 aaTAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCa                                   >  1:2399652/1‑150 (MQ=255)
                                                                                                                   tAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCCCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCaaa                                 >  2:2439324/1‑150 (MQ=255)
                                                                                                                     aaTTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATc                               <  1:2453382/150‑1 (MQ=255)
                                                                                                                      aTTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATcc                              >  1:978450/1‑150 (MQ=255)
                                                                                                                       ttGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCca                             <  1:1887030/150‑1 (MQ=255)
                                                                                                                          cAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCcatca                          <  2:1088801/150‑1 (MQ=255)
                                                                                                                          cAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCcatca                          <  1:1602367/150‑1 (MQ=255)
                                                                                                                            gAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCCATCata                        <  1:2081003/150‑1 (MQ=255)
                                                                                                                              aaCCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCCATCATATg                      <  2:2270778/150‑1 (MQ=255)
                                                                                                                                  ccATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCCATCATATGCgcc                  <  2:1779627/150‑1 (MQ=255)
                                                                                                                                  aagtttgCCTGATGCACAATGTTTATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCCATCATATGCgcc                  >  1:1705403/7‑150 (MQ=255)
                                                                                                                                   cATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCCATCATATGCgcca                 <  1:2390502/150‑1 (MQ=255)
                                                                                                                                       aGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCCATCATATGCgccagcgg             <  1:305145/150‑1 (MQ=255)
                                                                                                                                          cTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCcat                            >  2:761468/1‑132 (MQ=255)
                                                                                                                                          cTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCcat                            <  1:761468/132‑1 (MQ=255)
                                                                                                                                             aTAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCCATCATATGCgccagcggccagcg       <  2:666780/150‑1 (MQ=255)
                                                                                                                                             aTAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCCAACATATGCgccagcggccagcg       <  2:1800870/150‑1 (MQ=255)
                                                                                                                                               aaaCATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCCATCATATGCGCCAGCGGCCAgcgcc     >  2:1912760/1‑149 (MQ=255)
                                                                                                                                               aaaCATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCCATCATATGCGCCAGCGGCCAGCGCa     >  2:445308/1‑150 (MQ=255)
                                                                                                                                               aaaCATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCCATCATATGCGCCAGCGGCCAGCGCa     <  2:2424527/150‑1 (MQ=255)
                                                                                                                                                aaCATTGTGCATCAGACAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCCATCATATGCGCCAGCGGCCAGCGCAc    >  2:972428/1‑150 (MQ=255)
                                                                                                                                                 aCATTGTGCATCAGGCAAACTTCAAGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAAAAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCCATCATATGCGCCAGCGGCCAGCGCAcc   >  1:400283/1‑150 (MQ=255)
                                                                                                                                                  cATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCCATCATATGCGCCAGCGGCCAGCGCACCt  <  1:722665/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCCATCATATGCGCCAGCGGCCAGCGCACCT  >  NC_000913/3525415‑3525710

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