Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,525,558 A→G 5.7% intergenic (+135/+31) 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,5580AG 5.7% 481.7 / 13.1 157intergenic (+135/+31)mrcA/nudEpeptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcA/ADP‑sugar diphosphatase NudE
Reads supporting (aligned to +/- strand):  ref base A (64/84);  new base G (3/6);  total (67/90)
Fisher's exact test for biased strand distribution p-value = 7.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

Only 100 of 166 total aligned reads displayed.

CAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCCATCATATGCGCCAGCGGCCAGCGCA  >  NC_000913/3525409‑3525707
                                                                                                                                                     |                                                                                                                                                     
caGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATa                                                                                                                                                       <  2:154462/150‑1 (MQ=255)
 aGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATaa                                                                                                                                                      <  2:967192/150‑1 (MQ=255)
  ggAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATaaa                                                                                                                                                     >  2:722665/1‑150 (MQ=255)
    aaTTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACa                                                                                                                                                   <  1:2325897/150‑1 (MQ=255)
      ttGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACAtt                                                                                                                                                 <  2:2264956/150‑1 (MQ=255)
        gTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATtgt                                                                                                                                               >  2:1656960/1‑150 (MQ=255)
         ttCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATtgtg                                                                                                                                              <  1:1838063/150‑1 (MQ=255)
          tCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCTTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGc                                                                                                                                             >  1:2074978/1‑150 (MQ=255)
           cTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCa                                                                                                                                            >  1:1084098/1‑150 (MQ=255)
           cTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCa                                                                                                                                            <  1:1218463/150‑1 (MQ=255)
           cTGATTAAAAAGGCGCTTAGGCGCCTTTTCAGTTTGCTGACAAAGTCAACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCa                                                                                                                                            >  2:1788247/1‑150 (MQ=255)
            tGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCAt                                                                                                                                           >  1:1088801/1‑150 (MQ=255)
             gATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATc                                                                                                                                          >  1:530074/1‑150 (MQ=255)
             aactgaaaaGGCGCCGAAGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAGCCCCATAAGCCTGATAAACATTGTGCATc                                                                                                                                          >  1:984583/6‑150 (MQ=255)
              aTTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCa                                                                                                                                         <  1:1612677/150‑1 (MQ=255)
               ttAAAAAGGCGCTTCGGCGCCTTTTCAGTTTACTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAg                                                                                                                                        >  2:764301/1‑150 (MQ=255)
                tgaaaaGGCGCCGAAGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAgg                                                                                                                                       >  2:168189/3‑150 (MQ=255)
                tAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATGCACAATGTTTATCAgg                                                                                                                                       <  2:453012/150‑1 (MQ=255)
                 aaaaaGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATGCACAATGTTTATCAGGc                                                                                                                                      <  1:810224/150‑1 (MQ=255)
                 aaaaaGGAGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCTGATACTGCGTGAACGCTTTATCGGGCTAACAAAATCTTGAATACTCTATTAATTGCAGAAACCCCATAAGCCTGATGCACAATGTTATCAGGct                                                                                                                                     <  2:1239360/150‑2 (MQ=255)
                  aaaaGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATGCACAATGTTTATCAGGct                                                                                                                                     <  2:451559/150‑2 (MQ=255)
                  aaaaGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCa                                                                                                                                     >  2:1691320/1‑150 (MQ=255)
                      ggCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACt                                                                                                                                 >  1:437265/1‑150 (MQ=255)
                      ggCGCTTAGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACt                                                                                                                                 >  2:1641147/1‑150 (MQ=255)
                         gcTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCTATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCa                                                                                                                              >  2:2094398/1‑150 (MQ=255)
                          cTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCAc                                                                                                                             <  1:301417/150‑1 (MQ=255)
                            tCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGc                                                                                                                           <  2:1450958/150‑1 (MQ=255)
                               gcgcCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCAtt                                                                                                                        <  2:1525436/150‑1 (MQ=255)
                                cgcCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCAttt                                                                                                                       >  2:1602367/1‑150 (MQ=255)
                                  ccTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTac                                                                                                                     >  1:1355469/1‑150 (MQ=255)
                                     tttCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACt                                                                                                                  >  1:640801/1‑150 (MQ=255)
                                        cAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGc                                                                                                               >  1:965069/1‑150 (MQ=255)
                                          gTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGccc                                                                                                             <  1:45116/150‑1 (MQ=255)
                                            ttGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCt                                                                                                           >  1:2400070/1‑150 (MQ=255)
                                              gCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGc                                                                                                         >  1:1463227/1‑150 (MQ=255)
                                                 gACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCt                                                                                                      >  1:2194438/1‑150 (MQ=255)
                                                 gACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCt                                                                                                      >  1:1809817/1‑150 (MQ=255)
                                                    aaaGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTc                                                                                                   >  2:305145/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                                                                                                <  1:749099/150‑1 (MQ=255)
                                                          cACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCAt                                                                                             >  1:2188912/1‑150 (MQ=255)
                                                              tGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTcgc                                                                                         <  1:592366/150‑1 (MQ=255)
                                                              tGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTcgc                                                                                         >  2:607735/1‑150 (MQ=255)
                                                              tGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTcgc                                                                                         <  1:252250/150‑1 (MQ=255)
                                                                tttATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCgcgc                                                                                       <  2:1519363/150‑1 (MQ=255)
                                                                tttATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGACAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCgcgc                                                                                       <  2:311282/150‑1 (MQ=255)
                                                                tgttagCCGGATAACGCGTTCACGCCGTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCgcgc                                                                                       >  1:2424527/6‑150 (MQ=255)
                                                                    agCCGGATAACGCGTTCACGCCGTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGa                                                                                   >  2:1412498/2‑150 (MQ=255)
                                                                     gCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAg                                                                                  >  1:519097/1‑150 (MQ=255)
                                                                      ccGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAgg                                                                                 >  1:98308/1‑150 (MQ=255)
                                                                      ccGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCAATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCCAACGCGCACGAgg                                                                                 >  2:1201562/1‑150 (MQ=255)
                                                                        ggATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGaa                                                                               >  2:2422135/1‑150 (MQ=255)
                                                                              ggCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAgcgc                                                                         >  1:1895698/1‑150 (MQ=255)
                                                                              ggCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCACTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAgg                                                                                 <  1:1808957/142‑1 (MQ=255)
                                                                                        cgTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGACCAGCGCACTGACAATg                                                               >  2:1880690/1‑150 (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)
                                                                                            aTCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTgcgcg                                                           <  1:764301/150‑1 (MQ=255)
                                                                                              ccGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCt                                                         <  1:430790/150‑1 (MQ=255)
                                                                                                  cTAACAAAATCGTGAAAACTCAATAAATTGCAGAAATCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCAt                                                     <  2:1954659/150‑1 (MQ=255)
                                                                                                    aaCAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTg                                                   <  1:1992632/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)
                                                                                                            cGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAgg                                           <  2:1084098/150‑1 (MQ=255)
                                                                                                             gTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCAATGAAGTCAggg                                          >  2:611389/1‑150 (MQ=255)
                                                                                                             gTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCGTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCCggg                                          >  2:2145099/1‑150 (MQ=255)
                                                                                                              tGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGt                                         <  1:938253/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                                <  2:508641/150‑1 (MQ=255)
                                                                                                                       aaTAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCa                                <  2:2317656/150‑1 (MQ=255)
                                                                                                                             ttGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCca                          <  1:1887030/150‑1 (MQ=255)
                                                                                                                                cAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCcatca                       <  1:1602367/150‑1 (MQ=255)
                                                                                                                                cAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCcatca                       <  1:1658937/150‑1 (MQ=255)
                                                                                                                                 aGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCcatcat                      <  2:250849/150‑1 (MQ=255)
                                                                                                                                 aGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCcatcat                      <  2:1326481/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)
                                                                                                                                           tttgCCTGATGCACAATGTTTATCAGGCAACCTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCGTCATTGAAGTCAGGGTCTTCCAGCAAATCCATCATATGCgccagc            >  2:2190838/4‑150 (MQ=255)
                                                                                                                                             aGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCCATCATATGCgccagcgg          <  1:305145/150‑1 (MQ=255)
                                                                                                                                               ccTGATGCACAATGTTTATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCCATCATATGCgccagcggcc        >  1:2435042/1‑150 (MQ=255)
                                                                                                                                               ccTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCCATCATATGCgccagcggcc        <  2:437265/150‑1 (MQ=255)
                                                                                                                                                cTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCcat                         >  2:761468/1‑132 (MQ=255)
                                                                                                                                                   aTAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCCAACATATGCgccagcggccagcg    <  2:1800870/150‑1 (MQ=255)
                                                                                                                                                    tAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCCATCATATGCGCCAGCGGCCAgcgc   >  1:891810/1‑150 (MQ=255)
                                                                                                                                                    caaaCATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCCATCATATGCGCCAGCGGCCAgcgc   >  1:966278/2‑150 (MQ=255)
                                                                                                                                                     aaaCATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCCATCATATGCGCCAGCGGCCAgcgcc  >  2:1912760/1‑149 (MQ=255)
                                                                                                                                                     aaaCATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCCATCATATGCGCCAGCGGCCAGCGCa  <  2:2424527/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCCATCATATGCGCCAGCGGCCAGCGCA  >  NC_000913/3525409‑3525707

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: