Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 475,971 C→G 16.3% intergenic (+20/+11) ybaY → / ← ybaZ outer membrane lipoprotein/excision repair protein, alkyltransferase‑like protein ATL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913475,9710CG16.3% 615.9 / 98.5 226intergenic (+20/+11)ybaY/ybaZouter membrane lipoprotein/excision repair protein, alkyltransferase‑like protein ATL
Reads supporting (aligned to +/- strand):  ref base C (91/98);  new base G (21/16);  total (112/114)
Fisher's exact test for biased strand distribution p-value = 3.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 234 total aligned reads displayed.

GCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCT  >  NC_000913/475872‑476070
                                                                                                   |                                                                                                   
gcgcAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGc                                                                                                     <  1:473921/100‑1 (MQ=255)
 cgcAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGcc                                                                                                    >  2:524012/1‑100 (MQ=255)
 cgcAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGcc                                                                                                    >  2:624132/1‑100 (MQ=255)
 cgcAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGc                                                                                                     >  1:2550038/1‑99 (MQ=255)
  gcAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGccc                                                                                                   >  1:49829/1‑100 (MQ=255)
    aaCGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCgg                                                                                                 <  1:1226682/100‑1 (MQ=255)
     aCGACTTCCGTACTTTCGGCTTCACCATCTCAGTTGAATCCGTCTTCCGCAGTTCCAGTTCCTCCGGAATGTTTCGCGCCGTTAGCCCTCTCCGCCCgg                                                                                                 <  1:403793/99‑1 (MQ=255)
         cTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGGGCGGAGAgg                                                                                            <  2:4263567/100‑1 (MQ=255)
         cTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGGGCGGAGAgg                                                                                            <  2:1485748/100‑1 (MQ=255)
         cTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAgg                                                                                            <  1:1332520/100‑1 (MQ=255)
         cTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAgg                                                                                            <  2:746515/100‑1 (MQ=255)
          tACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGGGCGGAGAggg                                                                                           <  2:3337766/100‑1 (MQ=255)
          tACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAggg                                                                                           >  2:2950414/1‑100 (MQ=255)
          tACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAggg                                                                                           <  2:1118305/100‑1 (MQ=255)
          tACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAggg                                                                                           <  2:1329188/100‑1 (MQ=255)
           aCCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCGTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGGGCggg                                                                                           >  2:1840054/1‑99 (MQ=255)
           aCCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGGGCGGAGAGGGt                                                                                          <  2:3350556/100‑1 (MQ=255)
           aCCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGt                                                                                          >  2:2258358/1‑100 (MQ=255)
            ccGTACCTCCGACTTCCCCAACTCAGGTGAATCCGTCTTCCCCAGTTCCAGCTCCTGCGCAATATTAAGCGACGTTAACCCTCTCCGCCCGCAGAggg                                                                                           <  1:3987241/98‑1 (MQ=255)
             cGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGGGCGGAGAGGGTta                                                                                        <  2:1858207/100‑1 (MQ=255)
               tACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTtagt                                                                                      <  2:4204523/100‑1 (MQ=255)
               tACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTtagt                                                                                      >  1:3568656/1‑100 (MQ=255)
               tACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTtagt                                                                                      <  2:3713320/100‑1 (MQ=255)
                 ccTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTtagtag                                                                                    >  1:1963597/1‑100 (MQ=255)
                 ccTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCAAGCCAATAAGCAATATTAAGCGACGTTAACACTCTCAGCCCGGAGAGGGTtagtag                                                                                    >  2:3077701/1‑100 (MQ=255)
                     cGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGTCGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTcc                                                                                <  1:2341822/100‑1 (MQ=255)
                     cGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTcc                                                                                <  2:2723401/100‑1 (MQ=255)
                     cGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTcc                                                                                <  1:2244347/100‑1 (MQ=255)
                      gACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCa                                                                               <  2:3884831/100‑1 (MQ=255)
                        cTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGc                                                                             <  2:1371330/100‑1 (MQ=255)
                        cTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGc                                                                             >  2:2504295/1‑100 (MQ=255)
                         ttCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCg                                                                            <  2:4003643/100‑1 (MQ=255)
                         ctCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCg                                                                            <  2:2676278/99‑1 (MQ=255)
                          tCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGa                                                                           >  2:262266/1‑100 (MQ=255)
                           cACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTACGCGACGCTTCCCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGAt                                                                          >  1:4171361/1‑100 (MQ=255)
                           cACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGag                                                                          >  1:1084679/1‑99 (MQ=255)
                           cACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGag                                                                          >  1:2676265/1‑99 (MQ=255)
                           cACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGAt                                                                          <  1:786224/100‑1 (MQ=255)
                                aCTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGc                                                                     >  2:2268199/1‑100 (MQ=255)
                                  tCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTg                                                                   <  1:1723625/100‑1 (MQ=255)
                                   cAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGc                                                                  >  2:48904/1‑100 (MQ=255)
                                    aGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCa                                                                 <  2:3211152/100‑1 (MQ=255)
                                     ggTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCaa                                                                >  2:216995/1‑100 (MQ=255)
                                      gTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGTTATCGCTGGaac                                                               <  2:2110704/100‑2 (MQ=255)
                                      gTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAg                                                               >  2:3954033/1‑100 (MQ=255)
                                          aTCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGa                                                           >  2:1052038/1‑100 (MQ=255)
                                           tCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGAt                                                          >  2:6549/1‑100 (MQ=255)
                                            ccGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGAtt                                                         <  1:4339879/100‑1 (MQ=255)
                                            ccGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGAtt                                                         >  2:1423661/1‑100 (MQ=255)
                                              gTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTg                                                       <  2:3434801/100‑1 (MQ=255)
                                                cTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGcc                                                     >  1:2972457/1‑100 (MQ=255)
                                                  tCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCg                                                   >  2:1629418/1‑100 (MQ=255)
                                                  tCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCg                                                   >  1:893202/1‑100 (MQ=255)
                                                   cGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGGGCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGc                                                  <  2:4157991/100‑1 (MQ=255)
                                                      cAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTc                                               >  2:1586403/1‑100 (MQ=255)
                                                       aGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTcc                                              <  1:2135/100‑1 (MQ=255)
                                                       aGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTcc                                              >  1:3410081/1‑100 (MQ=255)
                                                        gTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTccc                                             <  1:2701491/100‑1 (MQ=255)
                                                           ccAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGAt                                          <  1:309207/100‑1 (MQ=255)
                                                             aGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATAc                                        >  2:2038384/1‑100 (MQ=255)
                                                              gCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATAcc                                       >  2:822687/1‑100 (MQ=255)
                                                              gCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATAcc                                       >  1:2125112/1‑100 (MQ=255)
                                                                tCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCAt                                     <  2:2387126/100‑1 (MQ=255)
                                                                 ccTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATc                                    >  1:114350/1‑100 (MQ=255)
                                                                  cTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGGGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATca                                   >  1:4279879/1‑100 (MQ=255)
                                                                   tACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATcan                                  <  2:3038090/100‑2 (MQ=255)
                                                                   tACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATcac                                  <  1:668861/100‑1 (MQ=255)
                                                                    aCGCAATATTAAGCGACGTTAACCCGCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATcaca                                 <  1:3039563/100‑1 (MQ=255)
                                                                     cGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCacac                                >  1:3181380/1‑100 (MQ=255)
                                                                     cGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCacac                                >  1:3817968/1‑100 (MQ=255)
                                                                      gCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACAcc                               >  2:4232314/1‑100 (MQ=255)
                                                                      gCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACAcc                               >  2:4050766/1‑100 (MQ=255)
                                                                            ttAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGc                         <  2:3476305/100‑1 (MQ=255)
                                                                             tAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGcc                        <  1:19964/100‑1 (MQ=255)
                                                                             tAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGcc                        <  2:1811817/100‑1 (MQ=255)
                                                                               aGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAg                      >  1:3953165/1‑100 (MQ=255)
                                                                               aGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAg                      >  2:3494610/1‑100 (MQ=255)
                                                                                gCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGt                     <  1:1559456/100‑1 (MQ=255)
                                                                                  gACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCGGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTaa                   <  1:2846236/100‑1 (MQ=255)
                                                                                    cGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATg                 <  2:3678532/100‑1 (MQ=255)
                                                                                     gTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGc                >  2:2125179/1‑100 (MQ=255)
                                                                                      ttAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGcc               >  1:4217978/1‑100 (MQ=255)
                                                                                      ttAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGcc               <  1:1648969/100‑1 (MQ=255)
                                                                                       tAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCt              <  1:3299028/100‑1 (MQ=255)
                                                                                        aaCCCTCTCCGGGCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTg             >  2:2152218/1‑100 (MQ=255)
                                                                                        aaCCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTg             >  2:2593216/1‑100 (MQ=255)
                                                                                         aCCCTCTCCGGGCGGAGAGGGTTTGTAGTTCCAGCGATAACGCTGCAAGTGGATTTGCCCGCTTCGCGATACCATCACACCTTCTGACAGTAATGCCTGt            >  1:3523069/1‑100 (MQ=255)
                                                                                         aCCCTCTCCGGGCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGt            >  2:4065005/1‑100 (MQ=255)
                                                                                         aCCCTCTCCGGGCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGt            >  2:434596/1‑100 (MQ=255)
                                                                                         aCCCTCTCCGGGCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGt            >  2:2458806/1‑100 (MQ=255)
                                                                                          cccTCTCCGGGCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTc           >  2:3931569/1‑100 (MQ=255)
                                                                                          ccatctcCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTc           >  2:1687105/4‑100 (MQ=255)
                                                                                           ccTCTCCGGGCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCg          >  1:1450686/1‑100 (MQ=255)
                                                                                            ctctCCGGGCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGc         >  1:727853/1‑100 (MQ=255)
                                                                                            ctctCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGc         >  1:566226/1‑100 (MQ=255)
                                                                                             tctcCGCCCGGAGAGGGTTAGTAGTTACAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCt        >  2:2682356/1‑100 (MQ=255)
                                                                                              ctcCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTg       >  2:1719226/1‑100 (MQ=255)
                                                                                               tcCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTg       >  2:3609844/1‑99 (MQ=255)
                                                                                               tcCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGa      >  1:4311710/1‑100 (MQ=255)
                                                                                                   cccGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCt  >  1:245231/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCT  >  NC_000913/475872‑476070

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