Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,832,188 A→T 17.6% 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,1880AT17.6% 384.0 / 78.0 172intergenic (+22/+31)yicH/yicIputative inner membrane‑anchored periplasmic AsmA family protein/putative alpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base A (72/69);  new base T (17/13);  total (89/83)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.43e-02

Only 100 of 173 total aligned reads displayed.

AGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCG  >  NC_000913/3832090‑3832286
                                                                                                  |                                                                                                  
 gcagcaAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACt                                                                                                  >  1:3598577/1‑100 (MQ=255)
 gcagcaAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACt                                                                                                  >  1:1693459/1‑100 (MQ=255)
   agcaAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCa                                                                                                <  2:2863391/100‑1 (MQ=255)
    gcaAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAg                                                                                               >  1:586709/1‑100 (MQ=255)
    gcaAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAg                                                                                               >  2:2581380/1‑100 (MQ=255)
     caAGTGACTCAAACCATGAATGCTGGCATTGTTTCTAGCGGTGAAGTTACAGCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGa                                                                                              <  2:1081716/100‑1 (MQ=255)
     caAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGa                                                                                              <  2:1628882/100‑1 (MQ=255)
       aGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACg                                                                                            <  2:1864513/100‑1 (MQ=255)
         tGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTa                                                                                          <  1:2534031/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)
              cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCgg                                                                                     <  1:1875209/100‑1 (MQ=39)
              cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCgg                                                                                     >  2:389330/1‑100 (MQ=255)
              cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCgg                                                                                     >  1:2265075/1‑100 (MQ=255)
              cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCgg                                                                                     >  2:1204280/1‑100 (MQ=255)
              cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGGACGTCTGAGGCGTACCCgg                                                                                     >  1:1702156/1‑100 (MQ=255)
              cAAACCATGAATGCTGGCATCGTTTCCAGCAGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCgg                                                                                     <  2:80997/100‑1 (MQ=39)
               aaaCCATTAATGCTGGCATTGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGGGTCGTACCCGGc                                                                                    <  1:2390229/100‑1 (MQ=255)
               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)
                 aCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCGGcga                                                                                  <  1:2328663/100‑3 (MQ=38)
                 aCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCAt                                                                                  >  2:2577969/1‑100 (MQ=255)
                  ccATGAATGCTGGCATCGTTTCCAGCTGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTTGCCGGGTACGTCTGAGTCGTACCCGGcgag                                                                                 <  1:1533721/100‑4 (MQ=37)
                  ccATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCAtt                                                                                 <  1:1460035/100‑1 (MQ=255)
                   catgcatgCTGGCATAGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGCTAAACATCAATCTCGCCGGGTACGACTCAGACGTAGCCGGCATTc                                                                                <  2:3407092/100‑1 (MQ=255)
                   cATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTc                                                                                >  2:521019/1‑100 (MQ=255)
                   cATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTc                                                                                >  1:481422/1‑100 (MQ=255)
                     tGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCa                                                                              <  1:3219137/100‑1 (MQ=255)
                      gAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCAt                                                                             <  2:728586/100‑1 (MQ=255)
                         tGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGCGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCa                                                                           >  2:2371040/1‑99 (MQ=255)
                          gCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAAt                                                                         <  1:1849941/100‑1 (MQ=255)
                               cATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGAtt                                                                    <  2:1050592/100‑1 (MQ=255)
                               cATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGAtt                                                                    >  2:2489809/1‑100 (MQ=255)
                                 tCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTAc                                                                  <  1:1790434/100‑1 (MQ=255)
                                 tCGTTTCCAGCGGCGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTAc                                                                  >  1:2477273/1‑100 (MQ=255)
                                    tttCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAAc                                                               <  2:261967/100‑1 (MQ=255)
                                     ttCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACg                                                              <  1:2857407/100‑1 (MQ=255)
                                      tCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGt                                                             <  2:3422131/100‑1 (MQ=255)
                                      tCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGt                                                             >  1:871137/1‑100 (MQ=255)
                                       ccAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTa                                                            >  2:1764712/1‑100 (MQ=255)
                                       ccAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTa                                                            >  1:1290216/1‑100 (MQ=255)
                                        cAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAATCATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTaa                                                           >  2:1936794/1‑100 (MQ=255)
                                        cAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTaa                                                           >  2:3568716/1‑100 (MQ=255)
                                         aGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTTAt                                                          >  2:3639304/1‑100 (MQ=255)
                                         aGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAAt                                                          >  1:3263771/1‑100 (MQ=255)
                                         aGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAAt                                                          <  1:3142542/100‑1 (MQ=255)
                                         aGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCCCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAAt                                                          <  1:79525/100‑1 (MQ=255)
                                           cGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTg                                                        <  1:293214/100‑1 (MQ=255)
                                           cGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTg                                                        <  1:1223384/100‑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)
                                                      cATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCCTCAATAGATTACAACGTACTTGTCAGCGCATTc                                             >  2:714977/1‑100 (MQ=255)
                                                      cATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCCTTc                                             >  2:3021991/1‑100 (MQ=255)
                                                      cATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCAGTCCATCAATAGATTACAACGTAATTGTCAGCGCATTc                                             >  2:2451274/1‑100 (MQ=255)
                                                       aTCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTcc                                            <  2:2724929/100‑1 (MQ=255)
                                                        tCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTccc                                           <  1:2882067/100‑1 (MQ=255)
                                                        tCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTccc                                           <  1:214707/100‑1 (MQ=255)
                                                         cGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCt                                          >  1:2416424/1‑100 (MQ=255)
                                                         cGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCt                                          >  1:2803938/1‑100 (MQ=255)
                                                          gACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCtt                                         >  2:1335997/1‑100 (MQ=255)
                                                          gACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCtt                                         >  1:3254531/1‑100 (MQ=255)
                                                          gACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCtt                                         >  1:2979169/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)
                                                               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)
                                                                        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)
                                                                                tCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGcccc                     >  2:1268678/1‑98 (MQ=255)
                                                                                gatggaatgCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTg                   >  1:363782/9‑100 (MQ=37)
                                                                                 cAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGc                  <  1:2233970/100‑1 (MQ=255)
                                                                                 cAATCTCGCCGGCTACGACTCAGACGTACCCGGCATTCCATCCATAGATTACAACGTAAATGCCAGCACACTCCCTTGAGGCTTCACCACCAGCCCCtcc                  >  2:1779290/1‑98 (MQ=255)
                                                                                   ggaatgCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCCGCGCATTCCCTTCAGGCTTCACCACCAGCCCCTGCTc                >  2:221603/6‑100 (MQ=37)
                                                                                    tctcGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCa               >  1:1151322/1‑100 (MQ=255)
                                                                                    tctcGCCGGCTACGACTCAGACGTACCCGGCATTCCCTCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTc                <  2:2493938/99‑1 (MQ=255)
                                                                                    gaatgCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCa               >  1:1249765/5‑100 (MQ=38)
                                                                                      tcGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACt             <  1:712027/100‑1 (MQ=255)
                                                                                      atgCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACt             >  1:3424325/3‑100 (MQ=38)
                                                                                      atgCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACt             >  2:1831511/3‑100 (MQ=38)
                                                                                        gcagggTACGACTCAGACGGACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACttt           <  1:2720067/97‑1 (MQ=255)
                                                                                         ccGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTCATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTc          >  2:1434802/1‑100 (MQ=39)
                                                                                         ccGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTc          <  1:1627821/100‑1 (MQ=255)
                                                                                           gggTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAg        >  1:2608818/1‑100 (MQ=255)
                                                                                             gTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGcc      <  2:3376262/100‑1 (MQ=255)
                                                                                             gTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGcc      <  2:687887/100‑1 (MQ=255)
                                                                                              taagACTCAGACGTACCCGCCATTCCAGCAATAGATTACAACCTAATTGTCAGCGCATTCACTTGAGGCTGCTCCACCAGCCCCTGCTCACTTCCAg        <  2:214024/94‑1 (MQ=255)
                                                                                                cGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGc   <  1:75622/100‑1 (MQ=255)
                                                                                                 gACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCg  >  1:3561800/1‑100 (MQ=255)
                                                                                                  |                                                                                                  
AGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCG  >  NC_000913/3832090‑3832286

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: