Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,632,533 G→A 5.1% A601T (GCA→ACA)  slt → lytic murein transglycosylase, soluble

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,632,5330GA 5.1% 279.1 / 3.8 78A601T (GCA→ACA) sltlytic murein transglycosylase, soluble
Reads supporting (aligned to +/- strand):  ref base G (24/50);  new base A (2/2);  total (26/52)
Fisher's exact test for biased strand distribution p-value = 5.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

CAACATTGGCACCAGTTACCTGCAATATGTTTATCAGCAGTTTGGCAATAATCGTATTTTCTCCTCAGCAGCTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGACGTCGTTACTGATC  >  NC_000913/4632394‑4632672
                                                                                                                                           |                                                                                                                                           
cAACATTGGCACCAGTTACCTGCAATATGTTTATCAGCAGTTTGGCAATAATCGTATTTTCTCCTCAGCAGCTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTgg                                                                                                                                             <  1:571104/140‑1 (MQ=255)
cAACATTGGCACCAGTTACCTGCAATATGTTTATCAGCAGTTTGGCAATAATCGTATTTTCTCCTCAGCAGCTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTgg                                                                                                                                             <  1:571296/140‑1 (MQ=255)
   cATTGGCCCCAGTTACCTGCAATATGTTTATCAGCAGTTTGGCAATAATCGTATTTTCTCCTCAGCAGCTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCAt                                                                                                                                          <  1:96765/140‑1 (MQ=255)
   cATTGGCACCAGTTACCTGCTATATGTTTATCAGCAGTTTGGCAATAATCGTATTTTCTCCTCAGCAGCTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCAt                                                                                                                                          <  2:479744/140‑1 (MQ=255)
   cATTGGCACCAGTTACCTGCAATATGTTTATCAGCAGTTTGGCAATAATCGTATTTTCTCCTCAGCAGCTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCAt                                                                                                                                          <  2:353894/140‑1 (MQ=255)
   cATTGGCACCAGTTACCTGCAATATGTTTATCAGCAGTTTGGCAATAATCGTATTTTCTCCTCAGCAGCTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCAt                                                                                                                                          <  2:1020362/140‑1 (MQ=255)
   cATTGGCACCAGTTACATTCAATATGTTTATCAGCATTTTGTCAATAATCGTATTTTCTCCTCAGCAGCTTATAACGCCGGAACAGGGAGGGTGCGAACCTGGCTTGGCAACTGCGCCGGGCGTATCGACGCAGTGGCAt                                                                                                                                          <  1:677878/140‑1 (MQ=255)
     ttGGCACCTGTTACCTGGAATATGTGTATCAGCAGTTTGGCATTAATCGTATTTTCTGCGCAGCAGCTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGGTTGGCAACAGCGCCGGGCGTGTCGACGCAGTGGCAgtt                                                                                                                                        >  2:390099/1‑140 (MQ=255)
       ggCACCAGTTACCTGCAATATGTTTATCAGCAGTTTGGCAATAATCGTATTTTCTCCTCAGCAGCTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGt                                                                                                                                      <  1:279174/140‑1 (MQ=255)
          aCCAGTTACCTGCAATATGTTTATCAGCAGTTTGGCAATAATCGTATTTTCTCCTCAGCAGCTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCga                                                                                                                                   <  2:1099836/140‑1 (MQ=255)
             aGTTACCTGCAATATGTTTATCAGCAGTTTGGCAATAATCGTATTTTCTCCTCAGCAGCTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGagag                                                                                                                                >  2:192558/1‑140 (MQ=255)
                   cTGCCATATGTTTATCAGCAGTTTGGCAATAATCGTATTTTCTCCTCAGCAGCTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTcc                                                                                                                          <  2:351133/140‑1 (MQ=255)
                      cAATATGTTTATCAGCAGTTTGGCAATAATCGTATTTTCTCCTCAGCAGCTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCAtt                                                                                                                       >  1:1283231/1‑140 (MQ=255)
                       aaTATGTTTATCAGCAGTTTGGCAATAATCGTATTTTCTCCTCAGCAGCTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATtc                                                                                                                      >  2:219550/1‑140 (MQ=255)
                       aaTATGTTTATCAGCAGTTTGGCAATAATCGTATTTTCTCCTCAGCAGCTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATtc                                                                                                                      >  2:218358/1‑140 (MQ=255)
                             tttATCAGCAGTTTGGCAATAATCGTATTTTCTCCTCAGCAGCTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCgag                                                                                                                >  1:544981/1‑140 (MQ=255)
                               tATCAGCAGTTTGGCAATAATCGTATTTTCTCCTCAGCAGCTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGAc                                                                                                              <  1:263766/140‑1 (MQ=255)
                                 tCAGCAGTTTGGCAATAATCGTATTTTCTCCTCAGCAGCTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGAcgc                                                                                                            >  2:297829/1‑140 (MQ=255)
                                  cagcagTTTGGCAATACTCGTATTTTCTCCTCAGCAGCTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGAcgcg                                                                                                           <  2:858441/140‑1 (MQ=255)
                                  cagcagTTTGGCAATAATCGTATTTTCTCCTCAGCAGCTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGAcgcg                                                                                                           <  2:1260230/140‑1 (MQ=255)
                                   agcagTTTGGCAATAATCGTATTTTCTCCTCAGCAGCTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGAcgcgc                                                                                                          >  2:1515716/1‑140 (MQ=255)
                                     cagTTTGGCAATAATCGTATTTTCTCCTCAGCAGCTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCgg                                                                                                        >  2:13055/1‑140 (MQ=255)
                                     cagTTTGGCAATAATCGTATTTTCTCCTCAGCAGCTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCgg                                                                                                        >  2:320867/1‑140 (MQ=255)
                                      agTTTGGCAATAATCGTATTTTCTCCTCAGCAGCTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGt                                                                                                       <  1:23455/140‑1 (MQ=255)
                                              aataatCGTATTTTCTCCTCAGCAGCTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTgaa                                                                                               <  2:998048/140‑1 (MQ=255)
                                                taatCGTATTTTCTCCTCAGCAGCTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTgaaga                                                                                             <  2:1468562/140‑1 (MQ=255)
                                                taatCGTATTTTCTCCTCAGCAGCTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTgaaga                                                                                             <  2:1546360/140‑1 (MQ=255)
                                                  atCGTATTTTCTCCTCAGCAGCTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAAc                                                                                           >  2:1348619/1‑140 (MQ=255)
                                                  atCGTATTTTCTCCTCAGCAGCTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAAc                                                                                           >  2:1346621/1‑140 (MQ=255)
                                                     gTATTTTCTCCTCAGCAGCTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTg                                                                                        <  1:574947/140‑1 (MQ=255)
                                                     gTATTTTCTCCTCAGCAGCTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTg                                                                                        <  1:574955/140‑1 (MQ=255)
                                                     gTATTTTCTCCTCAGCAGCTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTg                                                                                        <  2:672084/140‑1 (MQ=255)
                                                     gTATTTTCTCCTCAGCAGCTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTg                                                                                        <  1:302043/140‑1 (MQ=255)
                                                     gTATTTTCTCCTCAGCAGCTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTg                                                                                        <  1:1573552/140‑1 (MQ=255)
                                                     gTATTTTCTCCTCAGCAGCTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTg                                                                                        <  1:915313/140‑1 (MQ=255)
                                                     gTATTTTCTCCTCAGCAGCTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTg                                                                                        <  1:429132/140‑1 (MQ=255)
                                                       aTTTTCTCCTCAGCAGCTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCt                                                                                      <  1:964223/140‑1 (MQ=255)
                                                        ttttCTCCTCAGCAGCTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTg                                                                                     <  1:1403022/140‑1 (MQ=255)
                                                        ttttCTCCTCAGCAGCTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTg                                                                                     <  1:1402987/140‑1 (MQ=255)
                                                          ttCTCCTCAGCAGCTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGc                                                                                   <  1:947602/140‑1 (MQ=255)
                                                               ctcAGCAGCTTATAACGCCGGACCAGGGCTGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGAAGTGGCATTTGTCGAGAGTATTCCATTCGCCGAGAAGCGCGGTTAGGTGAAGAACGTGCTGGCATATg                                                                              >  2:1588938/1‑140 (MQ=255)
                                                                tcAGCAGCTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGa                                                                             <  1:47134/140‑1 (MQ=255)
                                                                  agcagcTTATAACGCCGGTCCAGGGCGGGTGCGAACCTGGGTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACg                                                                           <  1:874968/140‑1 (MQ=255)
                                                                  agcagcTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACg                                                                           <  2:182559/140‑1 (MQ=255)
                                                                  agcagcTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACg                                                                           <  2:167166/140‑1 (MQ=255)
                                                                       cTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTtac                                                                      <  2:961119/140‑1 (MQ=255)
                                                                       cTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTtac                                                                      <  1:457899/140‑1 (MQ=255)
                                                                        ttATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTtact                                                                     <  2:264033/140‑1 (MQ=255)
                                                                            aaCGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCg                                                                 >  2:1538488/1‑140 (MQ=255)
                                                                            aaCGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCg                                                                 <  1:315587/140‑1 (MQ=255)
                                                                            aaCGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCg                                                                 <  1:315373/140‑1 (MQ=255)
                                                                            aaCGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGGGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCg                                                                 <  2:905907/140‑1 (MQ=255)
                                                                                 cGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTAtt                                                            >  2:1432435/1‑140 (MQ=255)
                                                                                 cGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTAtt                                                            >  2:237415/1‑140 (MQ=255)
                                                                                   gACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTc                                                          >  1:545850/1‑140 (MQ=255)
                                                                                   gACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGACATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGtt                                                                                                      <  1:202424/96‑1 (MQ=255)
                                                                                   gACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGACATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGtt                                                                                                      <  1:202293/96‑1 (MQ=255)
                                                                                   gACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGACATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGtt                                                                                                      >  2:202424/1‑96 (MQ=255)
                                                                                   gACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGACATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGtt                                                                                                      >  2:202293/1‑96 (MQ=255)
                                                                                      cAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATg                                                       <  1:1012708/140‑1 (MQ=255)
                                                                                           cGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGa                                                  <  2:1470521/140‑1 (MQ=255)
                                                                                           cGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGa                                                  <  2:114658/140‑1 (MQ=255)
                                                                                              gTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATaa                                               >  1:1624426/1‑140 (MQ=255)
                                                                                              gTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGAGCTGGCTAATGACGCTTACTACCGCTATTTCATGGGGGATaa                                               >  1:1649494/1‑140 (MQ=255)
                                                                                                gCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAAc                                             <  1:1169095/140‑1 (MQ=255)
                                                                                                gCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAAc                                             <  2:1250905/140‑1 (MQ=255)
                                                                                                gCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAAc                                             <  1:1524265/140‑1 (MQ=255)
                                                                                                    aCCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCCTGTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGGGAAGAACGTGCTGGCTTCTGACGCTTACGACCGCTATTTCATGGGGGATAAACCGAc                                         >  2:389790/1‑140 (MQ=255)
                                                                                                     ccTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACg                                        <  2:303722/140‑1 (MQ=255)
                                                                                                        ggCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTtg                                     >  1:1680431/1‑140 (MQ=255)
                                                                                                                  cAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCa                           <  1:732053/140‑1 (MQ=255)
                                                                                                                   aGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCAc                          >  1:437666/1‑140 (MQ=255)
                                                                                                                          ggCGTATCGACGCAGTGGCAGTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGGTCTGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCTCGGAAggg                   >  2:502142/1‑140 (MQ=255)
                                                                                                                             gTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGa                <  1:753702/140‑1 (MQ=255)
                                                                                                                                tCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGAcgt             <  2:263129/140‑1 (MQ=255)
                                                                                                                                tCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGAcgt             <  2:262944/140‑1 (MQ=255)
                                                                                                                                     gCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGACGTCGTTa        <  1:117517/140‑1 (MQ=255)
                                                                                                                                     gCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGACGTCGTTa        <  2:1291392/140‑1 (MQ=255)
                                                                                                                                      cAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGAcgcg                                                                                                           >  2:1290784/1‑40 (MQ=255)
                                                                                                                                      cAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGAcgcg                                                                                                           <  1:1290784/40‑1 (MQ=255)
                                                                                                                                           gCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGACGTCGTTACTGATc  >  1:734724/1‑140 (MQ=255)
                                                                                                                                           |                                                                                                                                           
CAACATTGGCACCAGTTACCTGCAATATGTTTATCAGCAGTTTGGCAATAATCGTATTTTCTCCTCAGCAGCTTATAACGCCGGACCAGGGCGGGTGCGAACCTGGCTTGGCAACAGCGCCGGGCGTATCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGACGTCGTTACTGATC  >  NC_000913/4632394‑4632672

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