Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,522,780 C→G 5.2% 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.2% 410.2 / 9.8 115intergenic (‑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 (52/57);  new base G (3/3);  total (55/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 119 total aligned reads displayed.

TCCAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTCTGCTGGGAGCA  >  NC_000913/3522633‑3522927
                                                                                                                                                   |                                                                                                                                                   
tcCAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCActct                                                                                                                                                    >  1:685325/1‑149 (MQ=255)
tcCAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCActct                                                                                                                                                    >  3:758073/1‑149 (MQ=255)
 ccAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTg                                                                                                                                                   >  2:379481/1‑149 (MQ=255)
  cAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGa                                                                                                                                                  <  1:403863/149‑1 (MQ=255)
  cAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGa                                                                                                                                                  <  1:455694/149‑1 (MQ=255)
       gCAACCGCCAGCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGGTGCCGCTTCTTGCAGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGt                                                                                                                                             <  1:555756/149‑1 (MQ=255)
       gCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGt                                                                                                                                             >  2:400801/1‑149 (MQ=255)
        cagccGGCAGCAGCGTTGGAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGGGTCTTGCTGTTGTAAATGGCAACCAATTTGCCAGATCTTAAATGCCATGGTGATGATCTCCTTATCACCCGTCACAGCGACGGGtg                                                                                                                                            <  2:594557/146‑2 (MQ=255)
        cAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGta                                                                                                                                            <  2:558490/149‑1 (MQ=255)
          aCCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGtata                                                                                                                                          <  1:19830/149‑1 (MQ=255)
            cGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTatat                                                                                                                                         <  3:190919/148‑1 (MQ=255)
            cGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTatat                                                                                                                                         >  4:190919/1‑148 (MQ=255)
               caccaGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAAt                                                                                                                                     <  2:202433/149‑1 (MQ=255)
                accaGAGTTGCAAAAAGCATTCTTTTACGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCAGTTGTAAATGCAAACCAAATTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATg                                                                                                                                    >  4:482737/1‑149 (MQ=255)
                   aGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGAc                                                                                                                                                 >  2:470277/1‑133 (MQ=255)
                   aGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGAc                                                                                                                                                 <  1:470277/133‑1 (MQ=255)
                      gTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATTCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTg                                                                                                                              >  4:347721/1‑149 (MQ=255)
                      gTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTg                                                                                                                              >  4:118653/1‑149 (MQ=255)
                                 cATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGTAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCAGAGTGACGGGtgataa                                                                                                                                        >  1:248262/1‑123 (MQ=255)
                                  aTTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTtat                                                                                                                  <  1:761688/149‑1 (MQ=255)
                                      ttttGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTAtacta                                                                                                              >  2:401901/1‑149 (MQ=255)
                                      ttttGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTAtacta                                                                                                              >  2:633420/1‑149 (MQ=255)
                                              cccGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCgt                                                                                                      >  2:279530/1‑149 (MQ=255)
                                                  cACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCgtttgt                                                                                                  <  4:711534/149‑1 (MQ=255)
                                                   aCGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCgtttgtt                                                                                                 >  1:355886/1‑149 (MQ=255)
                                                   aCGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCgtttgtt                                                                                                 <  3:320291/149‑1 (MQ=255)
                                                     gATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTta                                                                                               <  2:597811/149‑1 (MQ=255)
                                                         gCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATaaa                                                                                           <  2:565345/149‑1 (MQ=255)
                                                         gCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATaaa                                                                                           <  3:159029/149‑1 (MQ=255)
                                                             ccgctaccgctTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGc                                                                                       >  1:732439/1‑149 (MQ=255)
                                                                   ccgctTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAAt                                                                                 <  1:33390/149‑1 (MQ=255)
                                                                   ccgctTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAAt                                                                                 <  4:559422/149‑1 (MQ=255)
                                                                    cgctTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATg                                                                                >  1:132931/1‑149 (MQ=255)
                                                                     gctTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCAGAGTGACgg                                                                                                                                               <  2:713736/85‑1 (MQ=255)
                                                                     gctTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGa                                                                               >  3:80959/1‑149 (MQ=255)
                                                                     gctTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGa                                                                               <  1:435598/149‑1 (MQ=255)
                                                                       ttcttGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATTATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGaaa                                                                             <  4:485449/149‑1 (MQ=255)
                                                                       ttcttGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGaaa                                                                             <  3:728024/149‑1 (MQ=255)
                                                                       ttcttGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGaaa                                                                             >  3:142920/1‑149 (MQ=255)
                                                                            gCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATAGCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTaaa                                                                        <  4:437335/149‑1 (MQ=255)
                                                                              tgttgtAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATg                                                                      >  1:424756/1‑149 (MQ=255)
                                                                                 tgtAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGa                                                                   <  4:142920/149‑1 (MQ=255)
                                                                                        gCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTcc                                                            >  3:536657/1‑149 (MQ=255)
                                                                                         cAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCa                                                           <  2:732439/149‑1 (MQ=255)
                                                                                         cAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCa                                                           <  3:745356/149‑1 (MQ=255)
                                                                                         cAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCa                                                           >  1:528533/1‑149 (MQ=255)
                                                                                          aaaCCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAg                                                          >  2:150495/1‑149 (MQ=255)
                                                                                          aaaCCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACCGGTATATCAATGCGTTTGGCTTGCCTTTATACTACCTCGCGTTTGTTTATATACTGACCATATGATACTATATGGGAAAttttaag                                                          >  2:322097/1‑145 (MQ=255)
                                                                                           aaCCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGt                                                         >  1:360417/1‑149 (MQ=255)
                                                                                            aCCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTg                                                        >  1:105953/1‑149 (MQ=255)
                                                                                              cAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGaa                                                      <  1:755335/149‑1 (MQ=255)
                                                                                              cAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGaa                                                      <  4:86781/149‑1 (MQ=255)
                                                                                              cAATTTGCCAGATCTTAAATGCCATTGTGATGATATCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAACCTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGaa                                                      <  1:633420/149‑1 (MQ=255)
                                                                                                aTTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGt                                                    >  3:657265/1‑149 (MQ=255)
                                                                                                 tttGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGtt                                                   <  1:160287/149‑1 (MQ=255)
                                                                                                   tGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCg                                                 <  2:528533/149‑1 (MQ=255)
                                                                                                       aGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTaaa                                             <  2:98199/149‑1 (MQ=255)
                                                                                                       aGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTaaa                                             >  4:683153/1‑149 (MQ=255)
                                                                                                        gATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAg                                            >  4:74879/1‑149 (MQ=255)
                                                                                                        gATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAg                                            >  4:731982/1‑149 (MQ=255)
                                                                                                            ttAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTAtt                                        <  1:150495/149‑1 (MQ=255)
                                                                                                             tAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCTAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTAttt                                       <  3:432453/149‑1 (MQ=255)
                                                                                                             tAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTAttt                                       <  4:407418/149‑1 (MQ=255)
                                                                                                             tAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAAATAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTAttt                                       <  1:309735/149‑1 (MQ=255)
                                                                                                                aTGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTAtttttt                                    <  4:80959/149‑1 (MQ=255)
                                                                                                                   ccATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGAt                                 <  2:132931/149‑1 (MQ=255)
                                                                                                                    cATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATc                                >  1:660326/1‑149 (MQ=255)
                                                                                                                     aTTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATcc                               <  2:87811/149‑1 (MQ=255)
                                                                                                                     aTTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATcc                               <  3:382116/149‑1 (MQ=255)
                                                                                                                     aTTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATcc                               >  4:542110/1‑149 (MQ=255)
                                                                                                                      ttGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTc                                                             >  3:276676/1‑118 (MQ=255)
                                                                                                                      ttGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTc                                                             <  4:276676/118‑1 (MQ=255)
                                                                                                                         tgatgaTCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGc                           <  3:482737/149‑1 (MQ=255)
                                                                                                                         tgatgaTCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGc                           <  4:164204/149‑1 (MQ=255)
                                                                                                                          gatgatCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCa                          >  3:259437/1‑149 (MQ=255)
                                                                                                                           atgatCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAg                         >  1:677201/1‑149 (MQ=255)
                                                                                                                             gatCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTc                       >  2:464333/1‑149 (MQ=255)
                                                                                                                               tctcCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTg                     >  3:633719/1‑149 (MQ=255)
                                                                                                                                ctcCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGt                    >  2:520802/1‑149 (MQ=255)
                                                                                                                                 tcCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCt                      <  3:571300/146‑1 (MQ=255)
                                                                                                                                 tcCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCt                      >  4:571300/1‑146 (MQ=255)
                                                                                                                                 tcCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGtt                   <  2:211388/149‑1 (MQ=255)
                                                                                                                                 tcCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGtt                   >  4:68119/1‑149 (MQ=255)
                                                                                                                                 tcCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGtt                   <  1:299401/149‑1 (MQ=255)
                                                                                                                                  ccTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCCTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTg                  >  3:23759/1‑149 (MQ=255)
                                                                                                                                   cTTATCTCCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGc                 <  1:741742/149‑1 (MQ=255)
                                                                                                                                   cTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGc                 >  3:688634/1‑149 (MQ=255)
                                                                                                                                    ttATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCa                >  1:765910/1‑149 (MQ=255)
                                                                                                                                       tCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTc             >  4:730987/1‑149 (MQ=255)
                                                                                                                                        cACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTct            >  2:118623/1‑149 (MQ=255)
                                                                                                                                         aCCCGTCAGAGTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTctg           >  3:244771/1‑149 (MQ=255)
                                                                                                                                         aCCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTctg           <  3:23647/149‑1 (MQ=255)
                                                                                                                                           ccGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGaa                                                      >  4:552884/1‑104 (MQ=255)
                                                                                                                                           ccGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGaa                                                      <  3:552884/104‑1 (MQ=255)
                                                                                                                                             gTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTCTGCTgg       <  3:628498/149‑1 (MQ=255)
                                                                                                                                             gTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTCTGCTgg       <  4:756392/149‑1 (MQ=255)
                                                                                                                                              tCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAAc                                                                                          >  3:747175/1‑65 (MQ=255)
                                                                                                                                              tCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAAc                                                                                          <  4:747175/65‑1 (MQ=255)
                                                                                                                                                 ctctGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTCTGCTGGGAGc   >  4:622028/1‑149 (MQ=255)
                                                                                                                                                  tctGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTCTGCTGGGAGCa  >  3:564250/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TCCAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTCTGCTGGGAGCA  >  NC_000913/3522633‑3522927

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