Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,832,191 C→G 17.7% intergenic (+25/+28) 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,1910CG17.7% 432.2 / 85.0 170intergenic (+25/+28)yicH/yicIputative inner membrane‑anchored periplasmic AsmA family protein/putative alpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base C (71/69);  new base G (17/13);  total (88/82)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-02

Only 100 of 172 total aligned reads displayed.

AGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACC  >  NC_000913/3832093‑3832290
                                                                                                  |                                                                                                   
agcaAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCa                                                                                                    <  2:2863391/100‑1 (MQ=255)
 gcaAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAg                                                                                                   <  1:509183/100‑1 (MQ=255)
    aGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACg                                                                                                <  2:1864513/100‑1 (MQ=255)
      tGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTa                                                                                              <  1:2534031/100‑1 (MQ=255)
         cTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAccc                                                                                           >  1:3286308/1‑100 (MQ=255)
         cTCAAACCATGAATGATGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAccc                                                                                           >  1:687890/1‑100 (MQ=255)
          tCAAACCATGAATGCTGGCATCGTTTGCAGCGGTGAAGTTACAGCGACGGAGCCGGTGCAGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCg                                                                                          <  2:1090815/100‑1 (MQ=255)
          tCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCg                                                                                          <  1:2632277/100‑1 (MQ=255)
           cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCgg                                                                                         <  1:2981732/100‑1 (MQ=39)
           cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCgg                                                                                         <  2:2959563/100‑1 (MQ=39)
           cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCgg                                                                                         >  1:2265075/1‑100 (MQ=255)
           cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCgg                                                                                         >  2:389330/1‑100 (MQ=255)
           cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCgg                                                                                         >  2:1204280/1‑100 (MQ=255)
           cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCgg                                                                                         >  1:1869367/1‑100 (MQ=255)
           cAAACCATGAATGCTGGCATCGTTTCCAGCAGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCgg                                                                                         <  2:80997/100‑1 (MQ=39)
            aaaCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCGGc                                                                                        <  1:1101802/100‑1 (MQ=39)
            aaaCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCGGc                                                                                        <  2:1702148/100‑1 (MQ=39)
            aaaCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGc                                                                                        <  2:2327137/100‑1 (MQ=255)
             aaCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCTATCTCGCCGGGTACGTCTGAGTCGTACCCGGcg                                                                                       <  2:2405329/100‑2 (MQ=38)
             aaCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCa                                                                                       >  1:182864/1‑100 (MQ=255)
              aCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCGGcga                                                                                      <  1:2328663/100‑3 (MQ=38)
              aCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCAt                                                                                      >  2:2577969/1‑100 (MQ=255)
               ccATGAATGCTGGCATCGTTTCCAGCTGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTTGCCGGGTACGTCTGAGTCGTACCCGGcgag                                                                                     <  1:1533721/100‑4 (MQ=37)
               ccATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCAtt                                                                                     <  1:3210484/100‑1 (MQ=255)
               ccATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCAtt                                                                                     <  1:1460035/100‑1 (MQ=255)
                catgcatgCTGGCATAGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGCTAAACATCAATCTCGCCGGGTACGACTCAGACGTAGCCGGCATTc                                                                                    <  2:3407092/100‑1 (MQ=255)
                cATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTc                                                                                    >  2:521019/1‑100 (MQ=255)
                   gAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCAt                                                                                 <  2:1685668/100‑1 (MQ=255)
                   gAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCAt                                                                                 <  2:728586/100‑1 (MQ=255)
                     aTGCTGGCATCGTTTCCAGCGGTGAAGTTACGTCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGCCTCAGACGTACCCGGCATTCCATc                                                                                <  2:2783179/99‑1 (MQ=255)
                     aTGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCa                                                                               <  1:2758374/100‑1 (MQ=255)
                      tGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCaa                                                                              <  2:3553931/100‑1 (MQ=255)
                      tGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCa                                                                               <  1:2971252/99‑1 (MQ=255)
                        cTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGGAAACATCAATCTCGCCGGGTACGACCCAGACGTACCCCGAATTCCATCaa                                                                              >  2:1861809/1‑98 (MQ=255)
                          ggCATCGTTTCCAGCGGTGAAGTTACATCGACGTAGCCCGTGCGGTAAACATCAATCGCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGa                                                                          <  2:2379567/100‑1 (MQ=255)
                            cATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGAtt                                                                        >  2:2489809/1‑100 (MQ=255)
                              tCGTTTCCAGCGGCGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTAc                                                                      >  1:2477273/1‑100 (MQ=255)
                                 tttCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAAc                                                                   <  2:261967/100‑1 (MQ=255)
                                  ttCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACg                                                                  <  1:2857407/100‑1 (MQ=255)
                                    ccAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTa                                                                >  2:1764712/1‑100 (MQ=255)
                                    ccAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTa                                                                >  1:1290216/1‑100 (MQ=255)
                                     cAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTaa                                                               >  2:3568716/1‑100 (MQ=255)
                                      aGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTTAt                                                              >  2:3639304/1‑100 (MQ=255)
                                      aGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAAt                                                              <  1:3142542/100‑1 (MQ=255)
                                      aGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAAt                                                              >  1:3263771/1‑100 (MQ=255)
                                        cGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTg                                                            <  1:1223384/100‑1 (MQ=255)
                                        cGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTg                                                            <  1:293214/100‑1 (MQ=255)
                                          gcgAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTc                                                          <  1:2357144/98‑1 (MQ=255)
                                          gTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTc                                                          >  2:1570813/1‑100 (MQ=255)
                                              aGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAgcg                                                      <  1:3527186/100‑1 (MQ=255)
                                              aGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAgcg                                                      <  1:150227/100‑1 (MQ=255)
                                               gTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTTCATCAATAGATTACAACGTAATTGTCAgcgc                                                     >  2:2724158/1‑100 (MQ=255)
                                               gTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAgcgc                                                     >  2:266269/1‑100 (MQ=255)
                                                 tACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCAt                                                   >  1:894964/1‑100 (MQ=255)
                                                   cATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCCTCAATAGATTACAACGTACTTGTCAGCGCATTc                                                 >  2:714977/1‑100 (MQ=255)
                                                   cATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCCTTc                                                 >  2:3021991/1‑100 (MQ=255)
                                                   cATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCAGTCCATCAATAGATTACAACGTAATTGTCAGCGCATTc                                                 >  2:2451274/1‑100 (MQ=255)
                                                     tCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTccc                                               <  1:2882067/100‑1 (MQ=255)
                                                      cGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCt                                              >  1:2416424/1‑100 (MQ=255)
                                                          ggAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTg                                            <  2:1131201/98‑1 (MQ=255)
                                                          ggAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAg                                          <  1:2095413/100‑1 (MQ=255)
                                                          ggAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCCTGAg                                          <  1:1597124/100‑1 (MQ=255)
                                                           gAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAgg                                         <  2:1058071/100‑1 (MQ=255)
                                                           gAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAgg                                         <  1:2095910/100‑1 (MQ=255)
                                                            aGCCGGTGCCGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCCATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGc                                        <  2:879210/100‑1 (MQ=255)
                                                              ccGGTGGGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCtt                                      <  2:2540552/100‑1 (MQ=255)
                                                               nggTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTCTCAGCGCATTCCCTTGAGGCTTc                                     >  2:1935164/2‑100 (MQ=255)
                                                                ggTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTca                                    >  1:93966/1‑100 (MQ=255)
                                                                  tGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTcacc                                  <  2:2583196/100‑1 (MQ=255)
                                                                     ggTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCCTTCCCTTGAGGCTTcaccac                                >  1:3275100/1‑99 (MQ=255)
                                                                      gTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCAccacca                              >  2:304964/1‑100 (MQ=255)
                                                                        aaaCATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGc                            >  1:1860380/1‑100 (MQ=255)
                                                                         aaCATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGcc                           >  1:148363/1‑100 (MQ=255)
                                                                         aaCATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGcc                           <  2:586707/100‑1 (MQ=255)
                                                                           cATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGcccc                         >  1:325427/1‑100 (MQ=255)
                                                                             gatggaatgCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTg                       >  1:363782/9‑100 (MQ=37)
                                                                              cAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGc                      <  1:2233970/100‑1 (MQ=255)
                                                                               aaTCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCt                     <  1:1754845/100‑1 (MQ=255)
                                                                                ggaatgCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCCGCGCATTCCCTTCAGGCTTCACCACCAGCCCCTGCTc                    >  2:221603/6‑100 (MQ=37)
                                                                                aTCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTc                    <  1:1657993/100‑1 (MQ=255)
                                                                                 tctcGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCa                   >  1:1151322/1‑100 (MQ=255)
                                                                                 gaatgCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCa                   >  1:1249765/5‑100 (MQ=38)
                                                                                   tcGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACt                 <  1:712027/100‑1 (MQ=255)
                                                                                    tgCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACtt                >  1:811777/2‑100 (MQ=39)
                                                                                     gcagggTACGACTCAGACGGACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACttt               <  1:2720067/97‑1 (MQ=255)
                                                                                     accGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACttt               >  2:1399773/2‑100 (MQ=39)
                                                                                      ccGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTCATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTc              >  2:1434802/1‑100 (MQ=39)
                                                                                      ccGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTc              <  1:1627821/100‑1 (MQ=255)
                                                                                       cGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCa             >  2:2039723/1‑100 (MQ=255)
                                                                                       cGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCa             >  2:2506026/1‑100 (MQ=255)
                                                                                        gggTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAg            >  1:2608818/1‑100 (MQ=255)
                                                                                          gTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGcc          <  2:687887/100‑1 (MQ=255)
                                                                                           taagACTCAGACGTACCCGCCATTCCAGCAATAGATTACAACCTAATTGTCAGCGCATTCACTTGAGGCTGCTCCACCAGCCCCTGCTCACTTCCAg            <  2:214024/94‑1 (MQ=255)
                                                                                             cGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGc       <  1:75622/100‑1 (MQ=255)
                                                                                              gACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGc       >  2:1825055/1‑99 (MQ=255)
                                                                                              gACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCg      >  1:3561800/1‑100 (MQ=255)
                                                                                              gACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCg      >  1:3561845/1‑100 (MQ=255)
                                                                                               aCTCAGACGTACCCGGCATTCCATCAATAGATTCCAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCCCACTTTCAGCCTGc       >  2:13250/1‑98 (MQ=255)
                                                                                                  cAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAAcc  <  1:1916623/100‑1 (MQ=255)
                                                                                                  cAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAAcc  >  2:3386747/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
AGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACC  >  NC_000913/3832093‑3832290

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