Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,832,188 A→T 7.1% intergenic (+22/+31) yicH → / ← yicI putative inner membrane‑anchored periplasmic AsmA family protein/putative alpha‑glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,832,1880AT 7.1% 582.3 / 21.6 165intergenic (+22/+31)yicH/yicIputative inner membrane‑anchored periplasmic AsmA family protein/putative alpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base A (75/78);  new base T (6/6);  total (81/84)
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.96e-01

Only 100 of 168 total aligned reads displayed.

CCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGT  >  NC_000913/3832040‑3832335
                                                                                                                                                    |                                                                                                                                                   
cccGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGa                                                                                                                                                     <  3:307291/149‑1 (MQ=255)
  cGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACt                                                                                                                                                   <  2:214844/149‑1 (MQ=255)
      gAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGa                                                                                                                                               <  2:389669/149‑1 (MQ=255)
          ccTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTa                                                                                                                                           >  3:27718/1‑149 (MQ=255)
            tACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAcc                                                                                                                                         >  2:243161/1‑149 (MQ=255)
             aCGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAccc                                                                                                                                        <  3:152571/149‑1 (MQ=255)
              cGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCg                                                                                                                                       <  3:88765/149‑1 (MQ=255)
               ggTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCgg                                                                                                                                      >  1:65327/1‑149 (MQ=255)
                 ttAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCa                                                                                                                                    <  4:186559/149‑1 (MQ=255)
                    gCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTc                                                                                                                                 >  1:461243/1‑149 (MQ=255)
                      gggCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCGGcga                                                                                                                                   >  1:68279/1‑143 (MQ=255)
                      gggCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCa                                                                                                                               >  1:176673/1‑149 (MQ=255)
                       ggCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCAt                                                                                                                              <  1:431327/149‑1 (MQ=255)
                       ggCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCAt                                                                                                                              >  4:170494/1‑149 (MQ=255)
                       ggCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCAt                                                                                                                              <  3:343444/149‑1 (MQ=255)
                        gCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATc                                                                                                                             <  3:55115/149‑1 (MQ=255)
                         cAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCa                                                                                                                            >  2:424633/1‑149 (MQ=255)
                          aaTTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCaa                                                                                                                           <  2:444537/149‑1 (MQ=255)
                           aTTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGa                                                                                                                                               >  4:223602/1‑128 (MQ=255)
                           aTTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGa                                                                                                                                               <  3:223602/128‑1 (MQ=255)
                            ttGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATa                                                                                                                         <  1:114530/149‑1 (MQ=255)
                                    cGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCAtt                                                                                                                                  >  4:276739/1‑132 (MQ=255)
                                     gTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAAc                                                                                                                >  3:418/1‑149 (MQ=255)
                                         aCGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTaa                                                                                                            <  3:287627/149‑1 (MQ=255)
                                         aCGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTaa                                                                                                            <  2:139316/149‑1 (MQ=255)
                                          cgccgcAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAAt                                                                                                           <  4:427466/149‑1 (MQ=255)
                                            ccgcAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTg                                                                                                         >  3:109302/1‑149 (MQ=255)
                                              gcAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTc                                                                                                       >  4:10581/1‑149 (MQ=255)
                                               cAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCa                                                                                                      <  3:371177/149‑1 (MQ=255)
                                                    cagctAGGGACTCAAACCTAGAGTGCAGGCACCGGTTCCGGCGGTGAAACTGCAGCTACGGAGCCCGTGCCTTACTCACAAATCTCGCCGCGTACGTCACATACGCACCCGGCATTCCATCAATATATTACAACGTAATTGTCatctcg                                                                                                 <  4:179959/149‑6 (MQ=255)
                                                     agcaAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGGGCAt                                                                                                >  1:437035/1‑149 (MQ=255)
                                                         aGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTccc                                                                                            <  2:461243/149‑1 (MQ=255)
                                                          gTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCt                                                                                           >  2:197586/1‑149 (MQ=255)
                                                             aCTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGa                                                                                        <  3:465650/149‑1 (MQ=255)
                                                             aCTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGa                                                                                        <  3:169037/149‑1 (MQ=255)
                                                              cTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAccc                                                                                                                                        <  2:123783/100‑1 (MQ=255)
                                                              cTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAccc                                                                                                                                        >  1:123783/1‑100 (MQ=255)
                                                              cTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAg                                                                                       >  1:355310/1‑149 (MQ=255)
                                                              cTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCCCATTCCCTTGAg                                                                                       >  1:29130/1‑149 (MQ=255)
                                                               tCAACCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGa                                                                                                                                                     <  2:33792/86‑1 (MQ=255)
                                                               tCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGAc                                                                                                                                              >  3:643/1‑93 (MQ=255)
                                                               tCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGAc                                                                                                                                              <  4:643/93‑1 (MQ=255)
                                                               tCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAgg                                                                                      >  4:259791/1‑149 (MQ=255)
                                                               tCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAgg                                                                                      >  2:229579/1‑149 (MQ=255)
                                                               tCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAg                                                                                       <  2:398239/148‑1 (MQ=255)
                                                                  aaCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCtt                                                                                   >  3:27525/1‑149 (MQ=255)
                                                                     cATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTcac                                                                                >  4:181715/1‑149 (MQ=255)
                                                                              tGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAccc                                                                                                                                        <  4:208559/84‑1 (MQ=255)
                                                                              tGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAccc                                                                                                                                        >  3:208559/1‑84 (MQ=255)
                                                                               ggCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGcccc                                                                      <  4:393483/149‑1 (MQ=255)
                                                                                  aTCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATa                                                                                                                         <  2:9652/95‑1 (MQ=255)
                                                                                  aTCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGa                                                                                        <  1:306722/128‑1 (MQ=255)
                                                                                  aTCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGa                                                                                        >  2:306722/1‑128 (MQ=255)
                                                                                  aTCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGc                                                                   >  4:419914/1‑149 (MQ=255)
                                                                                  aTCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGc                                                                   >  1:202173/1‑149 (MQ=255)
                                                                                        tCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACtt                                                             <  3:163194/149‑1 (MQ=255)
                                                                                           aGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCa                                                          <  1:32314/149‑1 (MQ=255)
                                                                                               gTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCt                                                      <  1:229579/149‑1 (MQ=255)
                                                                                               gTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCt                                                      >  3:92169/1‑149 (MQ=255)
                                                                                                 gAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGc                                                    >  2:421127/1‑149 (MQ=255)
                                                                                                     ttACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGa                                                                                        <  2:141111/109‑1 (MQ=255)
                                                                                                     ttACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGa                                                                                        >  1:141111/1‑109 (MQ=255)
                                                                                                     ttACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAAc                                                <  4:287396/149‑1 (MQ=255)
                                                                                                         aTCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTc                                            >  3:45936/1‑149 (MQ=255)
                                                                                                         aTCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTc                                            <  2:289334/149‑1 (MQ=255)
                                                                                                          tCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCt                                           >  4:219923/1‑149 (MQ=255)
                                                                                                          tCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCt                                           >  3:324526/1‑149 (MQ=255)
                                                                                                             aCGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGc                                        <  4:88196/149‑1 (MQ=255)
                                                                                                              cGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCa                                       <  1:384868/149‑1 (MQ=255)
                                                                                                               ggAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAg                                      >  2:288777/1‑149 (MQ=255)
                                                                                                                    cGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCAt                                 >  3:346471/1‑149 (MQ=255)
                                                                                                                       tGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCCGCGAACCGTCTTGCAGACCCTTTa                              <  1:405073/149‑1 (MQ=255)
                                                                                                                         cGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACt                            <  2:324023/149‑1 (MQ=255)
                                                                                                                          agTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACtt                           <  1:421127/148‑1 (MQ=255)
                                                                                                                           gTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACttt                          >  4:122801/1‑149 (MQ=255)
                                                                                                                             aaaCATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACtt                           <  2:415971/146‑1 (MQ=255)
                                                                                                                             aaaCATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTca                        <  3:170494/149‑1 (MQ=255)
                                                                                                                                 aTCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTc                         >  1:161460/1‑144 (MQ=255)
                                                                                                                                 aTCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTc                         <  2:161460/144‑1 (MQ=255)
                                                                                                                                 aTCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCacaac                    >  2:67527/1‑149 (MQ=255)
                                                                                                                                 aTCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCacaac                    >  3:132075/1‑149 (MQ=255)
                                                                                                                                  tCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCacaac                    >  3:271686/1‑148 (MQ=255)
                                                                                                                                    aaTCTCGCCGGGTACGACTTAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACAtt                 <  4:324526/149‑1 (MQ=255)
                                                                                                                                    aaTCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACAtt                 <  3:419914/149‑1 (MQ=255)
                                                                                                                                     aTCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTca                                                                                 <  2:266272/84‑1 (MQ=255)
                                                                                                                                     aTCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTca                                                                                 >  1:266272/1‑84 (MQ=255)
                                                                                                                                     aTCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTca                        >  2:96047/1‑141 (MQ=255)
                                                                                                                                     aTCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTca                        <  1:96047/141‑1 (MQ=255)
                                                                                                                                         cGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCa            <  2:37972/149‑1 (MQ=255)
                                                                                                                                          gCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCa                                                          >  1:341567/1‑102 (MQ=255)
                                                                                                                                          gCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAg           <  4:27718/149‑1 (MQ=255)
                                                                                                                                           ccGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAgg          <  2:355310/149‑1 (MQ=255)
                                                                                                                                            cGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGa                                                                                        >  2:177113/1‑70 (MQ=38)
                                                                                                                                             gggTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGca        >  1:405759/1‑149 (MQ=255)
                                                                                                                                             gggTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGca        >  3:193854/1‑149 (MQ=255)
                                                                                                                                             gggTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGca        >  3:226546/1‑149 (MQ=255)
                                                                                                                                               gTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGcaca      <  4:346471/149‑1 (MQ=255)
                                                                                                                                                tACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAg     <  1:197586/149‑1 (MQ=255)
                                                                                                                                                  cGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAgtg   >  4:222155/1‑149 (MQ=255)
                                                                                                                                                   gACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAgtgt  >  2:72052/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGT  >  NC_000913/3832040‑3832335

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: