Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,632,658 A→G G642G (GGA→GGG slt → soluble lytic murein transglycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,632,6580AG98.5% 244.4 / ‑3.2 67coding (1926/1938 nt)sltsoluble lytic murein transglycosylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (25/41);  minor base . (0/1);  total (25/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.38e-01

TCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGACGTCGTTACTGATCCGCACGTTTATGATATGCTATCGTACTCTTTAGCGAGTACAACCGGGGGAGGCATTTTGCTTCCCCCGCTAACAATGGCGACATATTATGGCCCAACAATCACCCTATTCAGCAGC  >  NC_000913/4632522‑4632788
                                                                                                                                        |                                                                                                                                  
tCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGcgt                                                                                                                                 <  1:42221/140‑1 (MQ=255)
tCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGcgt                                                                                                                                 <  2:535601/140‑1 (MQ=255)
      cAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTAc                                                                                                                           <  2:404833/140‑1 (MQ=255)
         tGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGa                                                                                                                        >  2:621419/1‑140 (MQ=255)
            cATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGCGTCGTTACTGATCCg                                                                                                                    <  1:575431/140‑1 (MQ=255)
             aTTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTATTGATCCg                                                                                                                    <  1:23127/140‑1 (MQ=255)
             aTTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCg                                                                                                                    <  1:967547/140‑1 (MQ=255)
             aTTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCg                                                                                                                    <  1:621419/140‑1 (MQ=255)
             aTTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCg                                                                                                                    <  1:173910/140‑1 (MQ=255)
             aTTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCg                                                                                                                    <  2:1128867/140‑1 (MQ=255)
             aTTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCg                                                                                                                    <  1:1116155/140‑1 (MQ=255)
              tttGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGc                                                                                                                   <  2:29646/140‑1 (MQ=255)
               ttGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCa                                                                                                                  <  1:437291/140‑1 (MQ=255)
               ttGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCa                                                                                                                  <  1:314096/140‑1 (MQ=255)
               ttGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCa                                                                                                                  <  1:403864/140‑1 (MQ=255)
                 gTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCg                                                                                                                    <  2:942677/136‑1 (MQ=255)
                 gTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCg                                                                                                                    >  1:942677/1‑136 (MQ=255)
                     agagTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCACGTTTa                                                                                                            <  1:490223/140‑1 (MQ=255)
                       agTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCACGTTTATg                                                                                                          >  2:427905/1‑140 (MQ=255)
                         tATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCACGTTTATGat                                                                                                        >  2:278319/1‑140 (MQ=255)
                           ttCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATcc                                                                                                                     >  2:667778/1‑125 (MQ=255)
                           ttCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATcc                                                                                                                     <  1:667778/125‑1 (MQ=255)
                           ttCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCACGTTTATGatat                                                                                                      >  2:189195/1‑140 (MQ=255)
                           ttCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCACGTTTATGatat                                                                                                      >  2:732465/1‑140 (MQ=255)
                            tCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCACGTTTATGATATg                                                                                                     <  2:1139129/140‑1 (MQ=255)
                                ttCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCACGTTTATGATATGCTAt                                                                                                 >  2:456774/1‑140 (MQ=255)
                                 tctcCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCACGTTTATGATATGCTATc                                                                                                >  1:120579/1‑140 (MQ=255)
                                   tcCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCTACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCACGTTTATGATATGCTATCGt                                                                                              <  2:1022074/140‑1 (MQ=255)
                                   tcCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCACGtt                                                                                                              <  1:646484/124‑1 (MQ=255)
                                   tcCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCACGtt                                                                                                              >  2:646484/1‑124 (MQ=255)
                                   tcCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCACGTTTATGATATGCTATCGt                                                                                              <  2:869941/140‑1 (MQ=255)
                                           gcgcgGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCACGTTTATGATATGCTATCGTACTCTTTa                                                                                      <  1:667565/140‑1 (MQ=255)
                                            cgcgGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCACGTTTATGATATGCTATCGTACTCTTTAg                                                                                     >  1:844612/1‑140 (MQ=255)
                                                gTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCACGTTTATGATATGCTATCGTACTCTTTAGCGAg                                                                                 <  2:844612/140‑1 (MQ=255)
                                                 ttATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCACGTTTATGATATGCTATCGTACTCTTTAGCGAGt                                                                                <  1:278319/140‑1 (MQ=255)
                                                  tATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCACGTTTATGATATGCTATCGTACTCTTTAGCGAGTa                                                                               <  1:214230/140‑1 (MQ=255)
                                                           aaCGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCACGTTTATGATATGCTATCGTACTCTTTAGCGAGTACAACCgggg                                                                      >  1:192320/1‑140 (MQ=255)
                                                            aCGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCACGTTTATGATATGCTATCGTACTCTTTAGCGAGTACAACCggggg                                                                     <  2:632039/140‑1 (MQ=255)
                                                               tGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCACGTTTATGATATGCTATCGTACTCTTTAGCGAGTACAACCGGGGGAgg                                                                  <  1:702920/140‑1 (MQ=255)
                                                                gCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCACGTTTATGATATGCTATCGTACTCTTTAGCGAGTACAACCGGGGGAGGc                                                                 >  2:1012333/1‑140 (MQ=255)
                                                                gCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCACGTTTATGATATGCTATCGTACTCTTTAGCGAGTACAACCGGGGGAGGc                                                                 >  1:1130353/1‑140 (MQ=255)
                                                                  tGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCACGTTTATGATATGCTATCGTACTCTTTAGCGAGTACAACCGGGGGAGGCAt                                                               >  1:405385/1‑140 (MQ=255)
                                                                      ttATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCACGTTTATGATATGCTATCGTACTCTTTAGCGAGTACAACCGGGGGAGGCATTTTg                                                           <  2:192320/140‑1 (MQ=255)
                                                                              cTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCACGTTTATGATATGCTATCGTACTCTTTAGCGAGTACAACCGGGGGAGGCATTTTGCTTccccc                                                   >  2:467521/1‑140 (MQ=255)
                                                                              cTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCACGTTTATGATATGCTATCGTACTCTTTAGCGAGTACAACCGGGGGAGGCATTTTGCTTccccc                                                   >  2:539567/1‑140 (MQ=255)
                                                                               ttACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCACGTTTATGATATGCTATCGTACTCTTTAGCGAGTACAACCGGGGGAGGCATTTTGCTTCCCCCg                                                  <  1:456774/140‑1 (MQ=255)
                                                                                       gCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCACGTTTATGATATGCTATCGTACTCTTTAGCGAGTACAACCGGGGGAGGCATTTTGCTTCCCCCGCTAACAAt                                          >  1:128520/1‑140 (MQ=255)
                                                                                         tATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCACGTTTATGATATGCTATCGTACTCTTTAGCGAGTACAACCGGGGGAGGCATTTTGCTTCCCCCGCTAACAATgg                                        <  1:1012333/140‑1 (MQ=255)
                                                                                         tATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCACGTTTATGATATGCTATCGTACTCTTTAGCGAGTACAACCGGGGGAGGCATTTTGCTTCCCCCGCTAACAATgg                                        <  1:288059/140‑1 (MQ=255)
                                                                                           tttCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCACGTTTATGATATGCTATCGTACTCTTTAGCGAGTACAACCGGGGGAGGCATTTTGCTTCCCCCGCTAACAATGGCg                                      <  1:539567/140‑1 (MQ=255)
                                                                                             tCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCACGTTTATGATATGCTATCGTACTCTTTAGCGAGTACAACCGGGGGAGGCATTTTGCTTCCCCCGCTAACAATGGCGAc                                    <  1:524336/140‑1 (MQ=255)
                                                                                                   gggATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCACGTTTATGATATGCTATCGTACTCTTTAGCGAGTACAACCGGGGGAGGCATTTTGCTTCCCCCGCTAACAATGGCGACAtatta                              >  2:155462/1‑140 (MQ=255)
                                                                                                   gggATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCACGTTTATGATATGCTATCGTACTCTTTAGCGAGTACAACCGGGGGAGGCATTTTGCTTCCCCCGCTAACAATGGCGACAtatta                              <  1:189195/140‑1 (MQ=255)
                                                                                                     gATAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCACGTTTATGATATGCTATCGTACTCTTTAGCGAGTACAACCGGGGGAGGCATTTTGCTTCCCCCGCTAACAATGGCGACATATTATg                            >  2:278424/1‑140 (MQ=255)
                                                                                                       tAAACCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCACGTTTATGATATGCTATCGTACTCTTTAGCGAGTACAACCGGGGGAGGCATTTTGCTTCCCCCGCTAACAATGGCGACATATTATGGc                          >  1:1124360/1‑140 (MQ=255)
                                                                                                          aCCGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCACGTTTATGATATGCTATCGTACTCTTTAGCGAGTACAACCGGGGGAGGCATTTTGCTTCCCCCGCTAACAATGGCGACATATTATGGCCca                       >  2:9415/1‑140 (MQ=255)
                                                                                                            cGACGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCACGTTTATGATATGCTATCGTACTCTTTAGCGAGTACAACCGGGGGAGGCATTTTGCTTCCCCCGCTAACAATGGCGACATATTATGGCCcaac                     <  2:1124360/140‑1 (MQ=255)
                                                                                                              aCGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCACGTTTATGATATGCTATCGTACTCTTTAGCGAGTACAACCGGGGGAGGCATTTTGCTTCCCCCGCTAACAATGGCGACATATTATGGCCcaacaa                   <  1:732465/140‑1 (MQ=255)
                                                                                                              aCGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCACGTTTATGATATGCTATCGTACTCTTTAGCGAGTACAACCGGGGGAGGCATTTTGCTTCCCCCGCTAACAATGGCGACATATTATGGCCcaacaa                   <  2:120579/140‑1 (MQ=255)
                                                                                                              aCGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCACGTTTATGATATGCTATCGTACTCTTTAGCGAGTACAACCGGGGGAGGCATTTTGCTTCCCCCGCTAACAATGGCGACATATTATGGCCcaacaa                   <  1:427905/140‑1 (MQ=255)
                                                                                                               cGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCACGTTTATGATATGCTATCGTACTCTTTAGCGAGTACAACCGGGGGAGGCATTTTGCTTCCCCCGCTAACAATGGCGACATATTATGGCCcaacaa                   >  2:516718/1‑139 (MQ=255)
                                                                                                               cGTTGATGAGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCACGTTTATGATATGCTATCGTACTCTTTAGCGAGTACAACCGGGGGAGGCATTTTGCTTCCCCCGCTAACAATGGCGACATATTATGGCCcaacaa                   <  1:516718/139‑1 (MQ=255)
                                                                                                                   gatgaGCGCCACGGAATGGGGGCGTCGTTACTGATCCGCACGTTTATGATATGCTATCGTACTCTTTAGCGAGTACAACCGGGGGAGGCATTTTGCTTCCCCCGCTAACAATGGCGACATATTATGGCCCAACAATCAcc              <  1:155462/140‑1 (MQ=255)
                                                                                                                         cgcCACGGAATGGGGGCGTCGTTACTGATCCGCACGTTTATGATATGCTATCGTACTCTTTAGCGAGTACAACCGGGGGAGGCATTTTGCTTCCCCCGCTAACAATGGCGACATATTATGGCCCAACAATCACCCTATTc        >  2:871815/1‑140 (MQ=255)
                                                                                                                          gcCACGGAATGGGGGCGTCGTTACTGATCCGCACGTTTATGATATGCTATCGTACTCTTTAGCGAGTACAACCGGGGGAGGCATTTTGCTTCCCCCGCTAACAATGGCGACATATTATGGCCCAACAATCACCCTATTca       >  2:1068612/1‑140 (MQ=255)
                                                                                                                               ggAATGGGGGCGTCGTTACTGATCCGCACGTTTATGATATGCTATCGTACTCTTTAGCGAGTACAACCGGGGGAGGCATTTTGCTTCCCCCGCTAACAATGGCGACATATTATGGCCCAACAATCACCCTATTCagcagc  <  1:871815/140‑1 (MQ=255)
                                                                                                                               ggAATGGGGGCGTCGTTACTGATCCGCACGTTTATGATATGCTATCGTACTCTTTAGCGAGTACAACCGGGGGAGGCATTTTGCTTCCCCCGCTAACAATGGCGACATATTATGGCCCAACAATCACCCTATTCagcagc  <  1:467521/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                  
TCGACGCAGTGGCATTTGTCGAGAGTATTCCATTCTCCGAGACGCGCGGTTATGTGAAGAACGTGCTGGCTTATGACGCTTACTACCGCTATTTCATGGGGGATAAACCGACGTTGATGAGCGCCACGGAATGGGGACGTCGTTACTGATCCGCACGTTTATGATATGCTATCGTACTCTTTAGCGAGTACAACCGGGGGAGGCATTTTGCTTCCCCCGCTAACAATGGCGACATATTATGGCCCAACAATCACCCTATTCAGCAGC  >  NC_000913/4632522‑4632788

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