Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,522,780 C→G 5.4% intergenic (‑29/‑91) hofM ← / → mrcA DNA catabolic putative pilus assembly protein/penicillin‑binding protein 1a, murein transglycosylase and transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,522,7800CG 5.4% 467.4 / 13.8 129intergenic (‑29/‑91)hofM/mrcADNA catabolic putative pilus assembly protein/penicillin‑binding protein 1a, murein transglycosylase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base C (61/61);  new base G (5/2);  total (66/63)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

Only 100 of 134 total aligned reads displayed.

CTCCAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTCTGCTGG  >  NC_000913/3522632‑3522922
                                                                                                                                                    |                                                                                                                                              
   cAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGa                                                                                                                                             <  3:1060778/149‑1 (MQ=255)
     gCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGAc                                                                                                                                            >  4:12105/1‑148 (MQ=255)
     gCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGAc                                                                                                                                            <  3:12105/148‑1 (MQ=255)
          aaCCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGAc                                                                                                                                            <  1:905740/143‑1 (MQ=255)
          aaCCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGAc                                                                                                                                            >  2:905740/1‑143 (MQ=255)
          aaCCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACgg                                                                                                                                          <  3:535509/145‑1 (MQ=255)
           aCCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGtata                                                                                                                                     <  3:1161142/149‑1 (MQ=255)
            ccGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTatat                                                                                                                                    <  1:1179557/149‑1 (MQ=255)
                   caGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCg                                                                                                                             >  1:1398260/1‑149 (MQ=255)
                            aaaaaGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGc                                                                                                                              <  1:1201422/139‑1 (MQ=255)
                            aaaaaGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGc                                                                                                                              >  2:1201422/1‑139 (MQ=255)
                             aaaaGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACggg                                                                                                                                         <  2:684072/127‑1 (MQ=255)
                             aaaaGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACggg                                                                                                                                         >  1:684072/1‑127 (MQ=255)
                              aaaGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCaa                                                                                                                                 <  2:1118535/134‑1 (MQ=255)
                                aGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGc                                                                                                                              <  1:1330570/135‑1 (MQ=255)
                                aGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGc                                                                                                                              >  2:1330570/1‑135 (MQ=255)
                                aGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCtt                                                                                                                <  1:196760/149‑1 (MQ=255)
                                   aTTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCAGAGTGACGGGtg                                                                                                                                       >  4:298181/1‑122 (MQ=255)
                                   aTTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCAGAGTGACGGGt                                                                                                                                        >  3:991730/1‑122 (MQ=255)
                                    ttCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTtata                                                                                                            >  3:677117/1‑149 (MQ=255)
                                       ttttGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCTTTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCaa                                                                                                                                 <  2:500874/125‑1 (MQ=255)
                                       ttttGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCaa                                                                                                                                 >  1:500874/1‑125 (MQ=255)
                                             gcCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTtat                                                                                                             >  1:646092/1‑139 (MQ=255)
                                             gcCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGCATATCAATGCGTCTGGCTTGCCTTtat                                                                                                             <  2:646092/139‑1 (MQ=255)
                                               cccGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCgt                                                                                                 >  2:1228308/1‑149 (MQ=255)
                                                 cgaaCGATCGCAACCGCTACCGCTTCTTGCTGTTATAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCAGAGTGACggg                                                                                                                                         >  3:971875/4‑107 (MQ=255)
                                                  gCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACgg                                                                                                                                          <  4:245188/105‑1 (MQ=255)
                                                  gCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACgg                                                                                                                                          >  3:245188/1‑105 (MQ=255)
                                                   cACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCgtttgt                                                                                             <  4:677117/149‑1 (MQ=255)
                                                    aCGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCgtttgtt                                                                                            <  2:1398260/149‑1 (MQ=255)
                                                      gATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTtat                                                                                                             <  2:1263613/130‑1 (MQ=255)
                                                       aTCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCAGAGTGACggg                                                                                                                                         >  2:645139/1‑101 (MQ=255)
                                                       aTCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGc                                                                                                                              <  1:1273450/112‑1 (MQ=255)
                                                       aTCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTtat                                                                                         <  4:1453333/149‑1 (MQ=255)
                                                            aaCCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATaaa                                                                                      <  3:493856/147‑1 (MQ=255)
                                                            aaCCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATaaa                                                                                      >  4:493856/1‑147 (MQ=255)
                                                               cgctaccgctTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGcc                                                                                 >  1:1000268/1‑149 (MQ=255)
                                                                  taccgctTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGcc                                                                                 <  3:244437/146‑1 (MQ=255)
                                                                     cgctTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATg                                                                           >  4:921138/1‑149 (MQ=255)
                                                                      gctTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTAtactac                                                                                                        >  3:726952/1‑119 (MQ=255)
                                                                      gctTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTAtactac                                                                                                        <  4:726952/119‑1 (MQ=255)
                                                                        ttcttGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGaaa                                                                        <  1:1228308/149‑1 (MQ=255)
                                                                        ttcttGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGaaa                                                                        <  3:118512/149‑1 (MQ=255)
                                                                          cttGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACt                                                                      >  4:554429/1‑149 (MQ=255)
                                                                             gCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTaaa                                                                   >  3:489692/1‑149 (MQ=255)
                                                                             gCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTaaa                                                                   <  3:921138/149‑1 (MQ=255)
                                                                                 ttgtAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGc                                                                                                                   <  3:1210916/97‑1 (MQ=255)
                                                                                 ttgtAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTACAAggg                                                               >  4:944033/1‑149 (MQ=255)
                                                                                 ttgtAAATGCAAACCAATTTGCCAGATCTTAAATGCCAGTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGc                                                                                                                   >  4:1210916/1‑97 (MQ=255)
                                                                                 ttgtAAATGCAAACCAATTTGCCAGATCGTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTaaa                                                                   <  4:667465/145‑1 (MQ=255)
                                                                                     aaaTGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGc                                                                                                                              >  3:1385900/1‑82 (MQ=255)
                                                                                     aaaTGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCACCTTATCACCCGTCACTCTGACGGGTATATCAATGc                                                                                                                              <  4:1385900/82‑1 (MQ=255)
                                                                                       aTGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAAttt                                                         >  1:342928/1‑149 (MQ=255)
                                                                                        tGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTc                                                        <  4:489692/149‑1 (MQ=255)
                                                                                            aaCCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGt                                                    >  1:1072433/1‑149 (MQ=255)
                                                                                              ccAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCgtttgtt                                                                                            >  1:1105411/1‑107 (MQ=255)
                                                                                              ccAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCgtttgtt                                                                                            <  2:1105411/107‑1 (MQ=255)
                                                                                               cAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGAc                                                                                                                                            <  1:1120840/58‑1 (MQ=255)
                                                                                               cAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGaa                                                 <  2:1000268/149‑1 (MQ=255)
                                                                                                aaTTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGtat                                                                                                                                      <  1:1036594/63‑1 (MQ=255)
                                                                                                aaTTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGtat                                                                                                                                      >  2:1036594/1‑63 (MQ=255)
                                                                                                aaTTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAg                                                >  2:204575/1‑149 (MQ=255)
                                                                                                     gCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCGGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGt                                           >  1:1437445/1‑149 (MQ=255)
                                                                                                      ccAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCgtttgt                                                                                             <  4:355191/98‑1 (MQ=255)
                                                                                                      ccAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCgtttgt                                                                                             >  3:355191/1‑98 (MQ=255)
                                                                                                      ccAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAAc                                                                                     >  2:896785/1‑106 (MQ=255)
                                                                                                      ccAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAAc                                                                                     <  1:896785/106‑1 (MQ=255)
                                                                                                        aGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTaaa                                        <  1:1387990/149‑1 (MQ=255)
                                                                                                         gATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCt                                                                                                                 <  3:1295514/75‑1 (MQ=255)
                                                                                                          aTCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAAt                                                                  >  4:1241270/1‑121 (MQ=255)
                                                                                                           tCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTtata                                                                                                            <  4:543968/78‑1 (MQ=255)
                                                                                                           tCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTtata                                                                                                            >  3:543968/1‑78 (MQ=255)
                                                                                                           tCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGaaa                                                                        <  2:1385781/114‑1 (MQ=255)
                                                                                                            cTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTAt                                    >  1:250472/1‑149 (MQ=255)
                                                                                                              tAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTAttt                                  <  1:204575/149‑1 (MQ=255)
                                                                                                               aaaTGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTAtttt                                 <  3:554429/149‑1 (MQ=255)
                                                                                                                aaTGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTAtttt                                 >  2:1149927/1‑148 (MQ=255)
                                                                                                                aaTGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTAtttt                                 <  1:1149927/148‑1 (MQ=255)
                                                                                                                  tGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTAtacta                                                                                                         <  1:858154/74‑1 (MQ=255)
                                                                                                                  tGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTAtacta                                                                                                         >  2:858154/1‑74 (MQ=255)
                                                                                                                     cATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCgtttgtt                                                                                            >  3:617513/1‑84 (MQ=255)
                                                                                                                          tgatgaTCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGc                      >  3:1377290/1‑149 (MQ=255)
                                                                                                                               atCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCt                 <  4:1377290/149‑1 (MQ=255)
                                                                                                                                  tcCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTcc                                                       >  4:11149/1‑108 (MQ=255)
                                                                                                                                  tcCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTcc                                                       <  3:11149/108‑1 (MQ=255)
                                                                                                                                  tcCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGtt              >  2:1095227/1‑149 (MQ=255)
                                                                                                                                   ccTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTAttttt                                <  4:536138/130‑1 (MQ=255)
                                                                                                                                    cTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTc                                                        >  2:1132893/1‑105 (MQ=255)
                                                                                                                                    cTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTc                                                        <  1:1132893/105‑1 (MQ=255)
                                                                                                                                     ttATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCgtttgt                                                                                             >  1:103312/1‑67 (MQ=255)
                                                                                                                                     ttATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCgtttgt                                                                                             <  2:103312/67‑1 (MQ=255)
                                                                                                                                       aTCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCAtt         >  2:1137029/1‑149 (MQ=255)
                                                                                                                                        tCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTc        >  1:736417/1‑149 (MQ=255)
                                                                                                                                          aCCCGTCAGAGTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGAt                            <  2:180456/127‑1 (MQ=255)
                                                                                                                                             cGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGa                                                              >  4:454215/1‑90 (MQ=255)
                                                                                                                                              gTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTCTGCTgg  <  1:1095227/149‑1 (MQ=255)
                                                                                                                                                   tctGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCtt                        <  1:1323467/122‑1 (MQ=255)
                                                                                                                                                   tctGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCtt                        >  2:1323467/1‑122 (MQ=255)
                                                                                                                                                   tctGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTctgc     >  1:423228/1‑141 (MQ=255)
                                                                                                                                                   tctGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTctgc     <  2:423228/141‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
CTCCAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTCTGCTGG  >  NC_000913/3522632‑3522922

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