Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,831,764 C→G 6.7% A436A (GCC→GCG yicH → putative inner membrane‑anchored periplasmic AsmA family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,831,7640CG 6.7% 978.1 / 44.9 270A436A (GCC→GCGyicHputative inner membrane‑anchored periplasmic AsmA family protein
Reads supporting (aligned to +/- strand):  ref base C (123/129);  new base G (9/9);  total (132/138)
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.27e-01

Only 100 of 299 total aligned reads displayed.

CGCCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAG  >  NC_000913/3831620‑3831911
                                                                                                                                                |                                                                                                                                                   
        gCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCaa                                                                                                                                          >  4:453815/1‑148 (MQ=255)
        gCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCaa                                                                                                                                          <  3:424810/148‑1 (MQ=255)
        gCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTCAATGCCGCCGCCGCGACATTCaa                                                                                                                                          <  1:955650/148‑1 (MQ=255)
         cAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAAt                                                                                                                                         <  4:318861/148‑1 (MQ=255)
               gATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACa                                                                                                                                               <  3:993205/136‑1 (MQ=255)
               gATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAAt                                                                                                                                         <  1:222832/142‑1 (MQ=255)
                  cATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGAt                                                                                                                                >  1:158139/1‑148 (MQ=255)
                   aTTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATg                                                                                                                               <  2:246490/148‑1 (MQ=255)
                     tGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGccgccgccg                                                                                                                                                    >  2:656959/1‑125 (MQ=255)
                     tGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGccgccgccg                                                                                                                                                    <  1:656959/125‑1 (MQ=255)
                     tGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGtt                                                                                                                             <  2:852321/148‑1 (MQ=255)
                     tGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGtt                                                                                                                             <  4:852251/148‑1 (MQ=255)
                                      tCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCaa                                                                                                                                          <  3:1097608/118‑1 (MQ=255)
                                      tCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGgc                                                                                                            >  4:359587/1‑148 (MQ=255)
                                         cGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCgtg                                                                                                                                     <  1:93849/120‑1 (MQ=255)
                                         cGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCt                                                                                                         >  1:906926/1‑148 (MQ=255)
                                          gTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTtcgtcg                                                                                                                        >  3:763310/1‑132 (MQ=255)
                                            ggCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGccgccgccg                                                                                                                                                    <  2:147972/102‑1 (MQ=255)
                                                 cTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCa                                                                                                 <  1:320466/148‑1 (MQ=255)
                                                     ccGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCgtg                                                                                                                                     >  2:163218/1‑108 (MQ=255)
                                                     ccGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAAcag                                                                                             >  3:874873/1‑148 (MQ=255)
                                                        cgcTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATg                                                                                                                               >  3:1068640/1‑111 (MQ=255)
                                                            cGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCg                                                                                                    >  2:907557/1‑134 (MQ=255)
                                                            cGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACg                                                                                      >  3:853305/1‑148 (MQ=255)
                                                             gATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACgg                                                                                     <  3:112545/148‑1 (MQ=255)
                                                               tGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGc                                                                                                   <  3:60588/132‑1 (MQ=255)
                                                               tGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGc                                                                                                   >  4:60588/1‑132 (MQ=255)
                                                                      ggTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGa                                                                                                       <  2:702482/121‑1 (MQ=255)
                                                                       gTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAAc                                                                                               >  2:489496/1‑128 (MQ=255)
                                                                           cAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGATGAACGCCAACAGCAGCACgg                                                                                     <  4:302956/134‑1 (MQ=255)
                                                                            aaCCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACgg                                                                                     >  4:96459/1‑133 (MQ=255)
                                                                             aCCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTc                                                                                                                  >  3:955929/1‑103 (MQ=255)
                                                                               cTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGa                                                                                  <  4:894242/133‑1 (MQ=255)
                                                                                tGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTg                                                                    <  4:66034/146‑1 (MQ=255)
                                                                                  aCGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGtt                                                                                                                             <  3:273698/87‑1 (MQ=255)
                                                                                    gCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCa                                                              <  2:313177/148‑1 (MQ=255)
                                                                                        gTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTa                                                          <  2:939269/148‑1 (MQ=255)
                                                                                        gTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTa                                                          <  3:515194/148‑1 (MQ=255)
                                                                                         ttACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCa                                                              >  3:375137/1‑143 (MQ=255)
                                                                                         ttACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTAc                                                         <  4:745466/148‑1 (MQ=255)
                                                                                          tACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCgtgt                                                                                                                                    >  2:971779/1‑72 (MQ=255)
                                                                                          tACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCgt                                                                                                                                      >  2:362855/1‑70 (MQ=255)
                                                                                          tACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCgt                                                                                                                                      <  1:362855/70‑1 (MQ=255)
                                                                                           aCCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGt                                                                                                                   >  1:985758/1‑88 (MQ=255)
                                                                                            ccGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGa                                                      >  1:167952/1‑148 (MQ=255)
                                                                                              gATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCa                                                                            <  1:120642/124‑1 (MQ=255)
                                                                                              gATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGaa                                                                       <  1:286992/129‑1 (MQ=255)
                                                                                               aTCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGaa                                                                                 <  2:202508/118‑1 (MQ=255)
                                                                                                tCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCa                                                              >  1:300206/1‑136 (MQ=255)
                                                                                                tCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGaaaa                                                   <  2:219548/147‑1 (MQ=255)
                                                                                                tCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGaaaa                                                   >  2:190842/1‑147 (MQ=255)
                                                                                                   tAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGc                                                                                                                <  2:586889/83‑1 (MQ=255)
                                                                                                   tAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAatat                                                                              <  2:1112194/117‑1 (MQ=255)
                                                                                                   tAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAat                                                                                >  3:227959/1‑115 (MQ=255)
                                                                                                    aaaTGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTa                                                          >  1:448976/1‑136 (MQ=255)
                                                                                                    aaaTGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTa                                                          <  2:448976/136‑1 (MQ=255)
                                                                                                     aaTGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCAt                                             <  3:148485/148‑1 (MQ=255)
                                                                                                       tGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCa                                                              >  3:808732/1‑129 (MQ=255)
                                                                                                            cGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCgtgt                                                                                                                                    <  4:1073898/54‑1 (MQ=255)
                                                                                                             gTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTg                                                                                                                                  <  4:1104991/55‑1 (MQ=255)
                                                                                                               cTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCgtgt                                                                                                                                    >  4:394568/1‑51 (MQ=255)
                                                                                                                  gAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAatat                                                                              <  3:664478/102‑1 (MQ=255)
                                                                                                                  gAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTg                                                                    >  1:591534/1‑112 (MQ=255)
                                                                                                                  gAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCAttt                                                           <  2:776179/121‑1 (MQ=255)
                                                                                                                   aGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCAc                                                                                       >  2:108167/1‑92 (MQ=255)
                                                                                                                     tGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAAcag                                                                                             <  1:212704/84‑1 (MQ=255)
                                                                                                                     tGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAAcag                                                                                             <  3:972492/84‑1 (MQ=255)
                                                                                                                     tGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGaaaa                                                   >  3:358579/1‑126 (MQ=255)
                                                                                                                        cTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCt                                                                                                               <  4:783019/63‑1 (MQ=255)
                                                                                                                        cTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATc                                                                             >  1:693551/1‑97 (MQ=255)
                                                                                                                          cGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACagcagc                                                                                         >  3:316424/1‑83 (MQ=255)
                                                                                                                          cGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCAttt                                                           <  4:342475/113‑1 (MQ=255)
                                                                                                                          cGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGa                                       >  4:171805/1‑133 (MQ=255)
                                                                                                                            gCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAAcagc                                                                                            <  1:782310/78‑1 (MQ=255)
                                                                                                                             cGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGa                                                                   <  2:508198/102‑1 (MQ=255)
                                                                                                                             cGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGa                                       <  2:298009/130‑1 (MQ=255)
                                                                                                                                 tCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAAc                                                                                               <  4:145705/70‑1 (MQ=255)
                                                                                                                                 tCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAg                                                                  >  4:99090/1‑99 (MQ=255)
                                                                                                                                  cTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTg                                                                                   >  1:5223/1‑81 (MQ=255)
                                                                                                                                  cTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAAc                                                                      >  1:491953/1‑94 (MQ=255)
                                                                                                                                  cTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTg                                                                    <  1:826842/96‑1 (MQ=255)
                                                                                                                                  cTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGcc                             <  4:687729/135‑1 (MQ=255)
                                                                                                                                  cTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACaaa                >  4:495792/1‑148 (MQ=255)
                                                                                                                                  cTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACaaa                >  4:564514/1‑148 (MQ=255)
                                                                                                                                   tGAATGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAata                                                                               >  3:703117/1‑84 (MQ=38)
                                                                                                                                   tGAATGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTa                                                          <  3:338875/105‑1 (MQ=255)
                                                                                                                                    gAATGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTgg                                        >  1:192133/1‑122 (MQ=255)
                                                                                                                                     aaTGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCg                               <  2:816130/130‑1 (MQ=255)
                                                                                                                                      aTGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACgg                                                                                     <  1:400747/75‑1 (MQ=38)
                                                                                                                                      aTGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACgg                                                                                     >  2:400747/1‑75 (MQ=38)
                                                                                                                                      aTGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGcc                             >  1:378180/1‑131 (MQ=255)
                                                                                                                                      aTGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGcc                             <  2:378180/131‑1 (MQ=255)
                                                                                                                                        gccgccgccGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAgtg                         <  1:1089347/133‑1 (MQ=255)
                                                                                                                                         tcgcggcggcGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACt                                                        >  1:721663/2‑101 (MQ=255)
                                                                                                                                           gccgccgcGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAgtgt                        <  4:870434/131‑1 (MQ=255)
                                                                                                                                            ccgccgcGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTg                                                                    <  3:470229/86‑1 (MQ=255)
                                                                                                                                            ccgccgcGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGaa                                      >  1:1062589/1‑116 (MQ=255)
                                                                                                                                               ccgcGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAAc                                                                                               <  4:259179/56‑1 (MQ=255)
                                                                                                                                               ccgcGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTc    <  1:43392/147‑1 (MQ=255)
                                                                                                                                               ccgcGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTc    >  2:43392/1‑147 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
CGCCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAG  >  NC_000913/3831620‑3831911

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