Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 837,650 T→C 7.5% intergenic (+214/+15) ybiC → / ← ybiJ putative dehydrogenase/DUF1471 family putative periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913837,6500TC 7.5% 213.5 / 3.4 64intergenic (+214/+15)ybiC/ybiJputative dehydrogenase/DUF1471 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base T (33/26);  new base C (3/2);  total (36/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAGCCAGTTTCGCTTCCAGGGCATCGAGGGTGGATGCGCCA  >  NC_000913/837512‑837797
                                                                                                                                          |                                                                                                                                                   
tCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTa                                                                                                                                           >  3:1765536/1‑149 (MQ=255)
   gCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGa                                                                                                                                        >  1:32014/1‑149 (MQ=255)
            cagGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACaaaatcgtc                                                                                                                                 >  3:553413/1‑141 (MQ=255)
                 aaaTACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAAtt                                                                                                                          >  4:1558680/1‑149 (MQ=255)
                       ttGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCa                                                                                                                    <  3:601894/149‑1 (MQ=255)
                          tgCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAAGGCCTTATCAGTCCTACGAATTACTTATAAATTACCGCa                                                                                                                    >  4:1137713/1‑146 (MQ=255)
                          tgCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCa                                                                                                                    <  3:1137713/146‑1 (MQ=255)
                               gACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACtt                                                                                                                                 <  2:448133/128‑1 (MQ=255)
                               gACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACtt                                                                                                                                 >  1:448133/1‑128 (MQ=255)
                                  gTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTaa                                                                                                         <  1:870805/149‑1 (MQ=255)
                                       gTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTAt                                                                                                    >  3:1008179/1‑149 (MQ=255)
                                       gTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTAt                                                                                                    >  1:1698588/1‑149 (MQ=255)
                                        tAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCa                                                                                                                    >  3:1557844/1‑132 (MQ=255)
                                        tAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCa                                                                                                                    <  4:1557844/132‑1 (MQ=255)
                                        tAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTAtt                                                                                                   <  3:1558680/149‑1 (MQ=255)
                                          aaCGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAAt                                                                                                                                      <  3:713180/112‑1 (MQ=255)
                                          aaCGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAAt                                                                                                                                      >  4:713180/1‑112 (MQ=255)
                                             gCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTAttgttgt                                                                                              <  4:1765536/149‑1 (MQ=255)
                                             gCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTAtt                                                                                                   >  3:1713630/1‑144 (MQ=255)
                                             gCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTAtt                                                                                                   <  4:1713630/144‑1 (MQ=255)
                                               cTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTAttgttgtt                                                                                             <  1:1136955/148‑1 (MQ=255)
                                               cTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTACTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTAttgttgtt                                                                                             >  2:1136955/1‑148 (MQ=255)
                                                  aTCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTtggtg                                                                                         >  1:611797/1‑149 (MQ=255)
                                                    ccGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGc                                                                                       <  2:611797/149‑1 (MQ=255)
                                                         cAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAg                                                                                  <  2:1698588/149‑1 (MQ=255)
                                                           aCGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAg                                                                                  >  4:69533/1‑147 (MQ=255)
                                                           aCGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAg                                                                                  <  3:69533/147‑1 (MQ=255)
                                                            cGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGa                                                                               <  4:1008179/149‑1 (MQ=255)
                                                                   gATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGt                                                                         <  3:317770/148‑1 (MQ=255)
                                                                   gATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGt                                                                         >  4:317770/1‑148 (MQ=255)
                                                                                 gATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCgg                                                          <  2:32014/149‑1 (MQ=255)
                                                                                   tAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGaa                                                                                                                                       <  3:557298/70‑1 (MQ=255)
                                                                                   tAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGaa                                                                                                                                       <  3:829829/70‑1 (MQ=255)
                                                                                   tAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGaa                                                                                                                                       >  4:557298/1‑70 (MQ=255)
                                                                                   tAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGaa                                                                                                                                       >  4:829829/1‑70 (MQ=255)
                                                                                      gACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTCCGAATTa                                                                                                                                    <  2:1018883/70‑1 (MQ=255)
                                                                                      gACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTa                                                                                                                                    >  1:1018883/1‑70 (MQ=255)
                                                                                                          ggCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAGCCAGttt                                 >  1:219188/1‑149 (MQ=255)
                                                                                                            cATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTtggt                                                                                          >  4:346649/1‑90 (MQ=255)
                                                                                                            cATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTtggt                                                                                          <  3:346649/90‑1 (MQ=255)
                                                                                                                    gacGGATAAGGCGTTTACGCCGCATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTtggtg                                                                                         >  4:1808219/3‑83 (MQ=25)
                                                                                                                    gCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAGCCAGTTTCGCTTCCAgg                       >  1:692079/1‑149 (MQ=255)
                                                                                                                          taaggcgTTTACGCCGCATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCt                                                        >  3:1361306/4‑110 (MQ=255)
                                                                                                                          taaggcgTTTACGCCGCATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCt                                                        >  1:1456583/4‑110 (MQ=255)
                                                                                                                          taaggcgTTTACGCCGCATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCt                                                        <  2:1456583/107‑1 (MQ=255)
                                                                                                                          taaggcgTTTACGCCGCATCCGTCCTACGAAGTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCt                                                        <  4:1361306/107‑1 (MQ=255)
                                                                                                                            cggcgTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAGCCAGTTTCGCTTCCAGGGCATCGAg               >  2:1393721/1‑149 (MQ=255)
                                                                                                                                 tAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTAttgttgtt                                                                                             <  2:1705882/66‑1 (MQ=255)
                                                                                                                                 tAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTAttgttgtt                                                                                             <  2:223455/66‑1 (MQ=255)
                                                                                                                                 tAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTAttgttgtt                                                                                             >  1:78033/1‑66 (MQ=255)
                                                                                                                                 tAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTAttgttgtt                                                                                             <  2:78033/66‑1 (MQ=255)
                                                                                                                                 tAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTAttgttgtt                                                                                             >  1:1705882/1‑66 (MQ=255)
                                                                                                                                 tAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTAttgttgtt                                                                                             >  1:223455/1‑66 (MQ=255)
                                                                                                                                  aaaCGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCa                                                                                                                    >  3:550535/1‑42 (MQ=255)
                                                                                                                                  aaaCGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCa                                                                                                                    <  4:550535/42‑1 (MQ=255)
                                                                                                                                   aaCGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAGCCAGTTTCGCTTCCAGGGCATCGa                <  2:1873415/141‑1 (MQ=255)
                                                                                                                                   aaCGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAGCCAGTTTCGCTTCCAGGGCATCGa                >  1:1873415/1‑141 (MQ=255)
                                                                                                                                   aaCGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAGCCAGTTTCGCTTCCAGGGCATCGAGGGTGGAt        >  1:1206289/1‑149 (MQ=255)
                                                                                                                                   aaCGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAGCCAGTTTCGCTTCCAGGGCATCGAGGGTGGAt        >  1:1177931/1‑149 (MQ=255)
                                                                                                                                   aaCGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAGCCAGTTTCGCTTCCAGGGCATCGAGGGTGGAt        >  3:1605329/1‑149 (MQ=255)
                                                                                                                                    aCGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAGCCAGTTTCGCTTCCTGGGCATCGAGGGTGGATg       >  1:1000376/1‑149 (MQ=255)
                                                                                                                                    aCGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAGCCAGTTTCGCTTCCAGGGCATCGAGGGTGGAt        <  1:1079574/148‑1 (MQ=255)
                                                                                                                                    aCGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAGCCAGTTTCGCTTCCAGGGCATCGAGGGTGGAt        >  2:1079574/1‑148 (MQ=255)
                                                                                                                                        cTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAGCCAGTTTCGCTTCCAGGGCATCGAGGGTGGATGCGcc   >  1:513093/1‑149 (MQ=255)
                                                                                                                                         ttATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGAt                                                                              >  4:1083702/1‑73 (MQ=255)
                                                                                                                                         ttATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGAt                                                                              <  3:1083702/73‑1 (MQ=255)
                                                                                                                                         ttATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCa                                         <  4:976527/110‑1 (MQ=255)
                                                                                                                                         ttATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCa                                         >  3:976527/1‑110 (MQ=255)
                                                                                                                                         ttATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAGCCAGTTTCGCTTCCAGGGCATCGAGGGTGGATGCGCCa  >  1:541789/1‑149 (MQ=255)
                                                                                                                                          tATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAGCCAGTTTCGCTTCCAGGGCATCGa                >  1:515350/1‑134 (MQ=255)
                                                                                                                                          tATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAGCCAGTTTCGCTTCCAGGGCATCGa                <  2:515350/134‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                   
TCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAGCCAGTTTCGCTTCCAGGGCATCGAGGGTGGATGCGCCA  >  NC_000913/837512‑837797

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑