Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,832,194 A→T 12.4% intergenic (+28/+25) 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,1940AT12.4% 604.5 / 60.8 193intergenic (+28/+25)yicH/yicIputative inner membrane‑anchored periplasmic AsmA family protein/putative alpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base A (86/83);  new base T (12/12);  total (98/95)
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.95e-01

Only 100 of 211 total aligned reads displayed.

TACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCA  >  NC_000913/3832052‑3832338
                                                                                                                                              |                                                                                                                                                
tACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGt                                                                                                                                               <  1:429316/146‑1 (MQ=255)
  cGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTAccc                                                                                                                                           <  2:1122749/148‑1 (MQ=255)
               aTTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTAcc                                                                                                                                            >  2:558624/1‑134 (MQ=255)
               aTTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGa                                                                                                                                                  <  4:851750/128‑1 (MQ=255)
                     tGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTAcc                                                                                                                                            <  4:507220/128‑1 (MQ=255)
                     tGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTAcc                                                                                                                                            >  3:507220/1‑128 (MQ=255)
                     tGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGa                                                                                                                                                  >  3:992132/1‑122 (MQ=255)
                     tGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAcc                                                                                                                                            <  2:98369/128‑1 (MQ=255)
                      gCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTAcc                                                                                                                                            <  1:612542/127‑1 (MQ=255)
                      gCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTAcc                                                                                                                                            >  2:612542/1‑127 (MQ=255)
                                  gcAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGt                                                                                                           >  3:115257/1‑148 (MQ=255)
                                   cAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCg                                                                                                                                          >  1:1076819/1‑116 (MQ=255)
                                   cAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCg                                                                                                                                          <  2:1076819/116‑1 (MQ=255)
                                          gcaAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTa                                                                                                                                              <  3:440077/105‑1 (MQ=255)
                                          gcaAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTa                                                                                                                                              >  4:440077/1‑105 (MQ=255)
                                            aaGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTAGGACTCAGa                                                                                                                                                  <  1:966782/99‑1 (MQ=255)
                                                    cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACa                                                                                                                     <  1:1049251/120‑1 (MQ=255)
                                                    cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACa                                                                                                                     >  2:1049251/1‑120 (MQ=255)
                                                    cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTcc                                                                                                >  3:207591/1‑141 (MQ=255)
                                                    cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTcc                                                                                                <  4:207591/141‑1 (MQ=255)
                                                      aaCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCt                                                                                       >  4:316621/1‑148 (MQ=255)
                                                        ccATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGa                                                                                                                                                  >  2:845420/1‑87 (MQ=255)
                                                        ccATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCAGGTACGACTCAGa                                                                                                                                                  <  1:845420/87‑1 (MQ=255)
                                                         cATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTa                                                                                                                                              >  4:992933/1‑90 (MQ=255)
                                                         cATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTAccc                                                                                                                                           <  3:958925/93‑1 (MQ=255)
                                                           tGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTAcc                                                                                                                                            >  2:691746/1‑90 (MQ=255)
                                                           tGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTAcc                                                                                                                                            <  1:691746/90‑1 (MQ=255)
                                                           tGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAg                                                                                                                           <  2:285705/107‑1 (MQ=255)
                                                               tGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGa                                                                                                                                                  >  3:541336/1‑80 (MQ=255)
                                                               tGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGa                                                                                                                                                  <  4:541336/80‑1 (MQ=255)
                                                               tGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTcc                                                                                                <  3:812017/130‑1 (MQ=255)
                                                                gCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTg                                                                                            >  3:133161/1‑133 (MQ=255)
                                                                gCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTg                                                                                            <  4:133161/133‑1 (MQ=255)
                                                                     cATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAcc                                                                                                                                            >  1:14792/1‑80 (MQ=255)
                                                                     cATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGccc                                                                          >  4:777846/1‑146 (MQ=255)
                                                                       tCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCa                                                                                                                                  <  1:131384/88‑1 (MQ=255)
                                                                       tCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCa                                                                                                                                  >  2:131384/1‑88 (MQ=255)
                                                                         gTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTc                                                                    <  3:316621/148‑1 (MQ=255)
                                                                               aGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCg                                                                                                                                          >  1:11578/1‑72 (MQ=255)
                                                                               aGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTc                                                              <  1:204947/148‑1 (MQ=255)
                                                                                     gAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTg                                                        >  1:4014/1‑148 (MQ=255)
                                                                                      aaGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAAt                                                                                                              <  2:116837/93‑1 (MQ=255)
                                                                                      aaGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAAt                                                                                                              >  1:116837/1‑93 (MQ=255)
                                                                                      aaGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGAAGTACCCGGCATTCCATCGATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCTCTTTCAGCCTGc                                                       <  4:987248/148‑1 (MQ=255)
                                                                                         ttACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTg                                                                                            >  3:897713/1‑108 (MQ=255)
                                                                                              tCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCa                                                                                                    >  3:767050/1‑95 (MQ=255)
                                                                                                gACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCtt                                             >  2:152521/1‑148 (MQ=255)
                                                                                                 aCGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTg                                                                                                            <  2:450306/84‑1 (MQ=255)
                                                                                                 aCGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTg                                                                                                            >  1:387592/1‑84 (MQ=255)
                                                                                                 aCGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTg                                                                                                            <  2:387592/84‑1 (MQ=255)
                                                                                                        cGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTcaccac                                                                                <  2:501218/105‑1 (MQ=255)
                                                                                                         ggTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCAt                                    <  4:870300/148‑1 (MQ=255)
                                                                                                         ggTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCAt                                    <  2:4014/148‑1 (MQ=255)
                                                                                                                  aaCATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTca                                                                                    >  4:122079/1‑91 (MQ=255)
                                                                                                                     aTCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCAttt                                  >  4:367726/1‑138 (MQ=255)
                                                                                                                     aTCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCAttt                                  <  3:367726/138‑1 (MQ=255)
                                                                                                                      tCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTa                                 <  1:876359/138‑1 (MQ=255)
                                                                                                                      tCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTa                                 >  2:876359/1‑138 (MQ=255)
                                                                                                                       cAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGAcc                                      <  3:835130/132‑1 (MQ=255)
                                                                                                                         aTCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCCTCAATAGa                                                                                                                          >  4:829237/1‑46 (MQ=255)
                                                                                                                          tctcGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCt                                                                                       >  3:1034827/1‑80 (MQ=255)
                                                                                                                          tctcGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCt                                                                                       <  4:1034827/80‑1 (MQ=255)
                                                                                                                             cGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTg                                                                                            <  3:785449/72‑1 (MQ=255)
                                                                                                                               ccGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAAc                                                   >  2:262042/1‑111 (MQ=255)
                                                                                                                                 gggTACGTCTGAGTCGTACCCGGCATTCCATCAATAGTTTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGa                                        <  4:766584/120‑1 (MQ=255)
                                                                                                                                  ggTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTg                                                                                            >  3:64969/1‑67 (MQ=37)
                                                                                                                                  ggTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTg                                                                                            <  4:64969/67‑1 (MQ=37)
                                                                                                                                  ggTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGt                                                >  4:944025/1‑111 (MQ=255)
                                                                                                                                    tACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTc                                                                                                          >  3:1064701/1‑51 (MQ=255)
                                                                                                                                    tACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTc                                                                                                          <  4:1064701/51‑1 (MQ=255)
                                                                                                                                     aCGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAg        <  4:680890/148‑1 (MQ=255)
                                                                                                                                     aCGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAg        >  3:680890/1‑148 (MQ=255)
                                                                                                                                      cGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGc                                           >  1:484554/1‑112 (MQ=255)
                                                                                                                                        aCTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCt                                              >  1:1090620/1‑107 (MQ=255)
                                                                                                                                        aCTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCt                                              <  2:1090620/107‑1 (MQ=255)
                                                                                                                                        aCTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAgtgt     >  4:794217/1‑148 (MQ=255)
                                                                                                                                         cTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTc                                                              >  1:519973/1‑90 (MQ=255)
                                                                                                                                          tCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAAc                                                   <  4:597579/100‑1 (MQ=255)
                                                                                                                                          tCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAAc                                                   >  3:597579/1‑100 (MQ=255)
                                                                                                                                          tCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACt                               <  3:276271/120‑1 (MQ=255)
                                                                                                                                          tCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACt                               >  4:276271/1‑120 (MQ=255)
                                                                                                                                          tCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGcac          >  4:395478/1‑141 (MQ=255)
                                                                                                                                          tCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGcac          <  3:395478/141‑1 (MQ=255)
                                                                                                                                          tCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAgtgt     >  1:826475/1‑146 (MQ=255)
                                                                                                                                          tCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAgtgt     <  2:826475/146‑1 (MQ=255)
                                                                                                                                           cAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTc                                                                                                 >  3:994209/1‑53 (MQ=255)
                                                                                                                                           cAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTcaccac                                                                                <  4:64170/70‑1 (MQ=255)
                                                                                                                                           cAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTcaccac                                                                                >  3:64170/1‑70 (MQ=255)
                                                                                                                                           cAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGAcc                                      <  4:659947/112‑1 (MQ=255)
                                                                                                                                           cAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGAcc                                      >  3:659947/1‑112 (MQ=255)
                                                                                                                                           cAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTa                                 >  3:852228/1‑117 (MQ=255)
                                                                                                                                           cAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACtt                              <  2:803806/120‑1 (MQ=255)
                                                                                                                                           cAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACtt                              >  1:803806/1‑120 (MQ=255)
                                                                                                                                           cAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTc                            <  1:365096/122‑1 (MQ=255)
                                                                                                                                           cAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTc                            >  2:365096/1‑122 (MQ=255)
                                                                                                                                           cAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAg        >  1:266401/1‑142 (MQ=255)
                                                                                                                                           cAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCa  >  4:1128936/1‑148 (MQ=255)
                                                                                                                                           cAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCa  >  1:560369/1‑148 (MQ=255)
                                                                                                                                             gACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTgcgc                >  1:158229/1‑132 (MQ=255)
                                                                                                                                              aCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGccc                                                                          <  2:164911/73‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                
TACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCA  >  NC_000913/3832052‑3832338

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: