Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,934,276 C→T H125Y (CAC→TAC)  rbsB → D‑ribose periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,934,2760CT99.2% 1400.8 / ‑4.4 393H125Y (CAC→TAC) rbsBD‑ribose periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (192/197);  total (193/200)
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.97e-01

Only 100 of 411 total aligned reads displayed.

GCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTG  >  NC_000913/3934186‑3934365
                                                                                          |                                                                                         
gCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCtac                                                                                         >  2:7104099/1‑93 (MQ=255)
 cAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCtaca                                                                                        >  1:6126494/1‑93 (MQ=255)
  aGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACAt                                                                                       <  2:2409331/93‑1 (MQ=255)
  aGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACAt                                                                                       >  2:3403391/1‑93 (MQ=255)
  aGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGAGAGCTACAt                                                                                       <  2:5988301/93‑1 (MQ=255)
   gTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACAtt                                                                                      >  2:1236747/1‑93 (MQ=255)
   gTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACAtt                                                                                      <  2:6677041/93‑1 (MQ=255)
   gTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACAtt                                                                                      <  2:2587692/93‑1 (MQ=255)
   gTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACAtg                                                                                      >  1:4429595/1‑92 (MQ=255)
     gggTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGc                                                                                    <  1:6287854/93‑1 (MQ=255)
        tAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTc                                                                                 <  1:777276/93‑1 (MQ=255)
        tAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTc                                                                                 <  1:4905045/93‑1 (MQ=255)
          aTGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTg                                                                               >  2:7165839/1‑93 (MQ=255)
          aTGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTg                                                                               <  2:5045471/93‑1 (MQ=255)
           tGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTACCACTCTTAACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCGACATTGCTTCTGa                                                                              <  2:7645758/93‑1 (MQ=255)
             cTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATa                                                                            >  2:1723984/1‑93 (MQ=255)
              tgtgAAGATGGCTAACCAGGCGAACATCCCGGTTATCCCTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATaa                                                                           >  1:6862348/1‑93 (MQ=255)
              tgtgAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATaa                                                                           <  1:4149100/93‑1 (MQ=255)
               gtgAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAAc                                                                          >  2:4251535/1‑93 (MQ=255)
                 gAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCCGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGt                                                                        <  1:3695705/93‑1 (MQ=255)
                  aaGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTa                                                                       <  1:2290747/93‑1 (MQ=255)
                  aaGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTa                                                                       <  2:3916681/93‑1 (MQ=255)
                   aGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTAc                                                                      <  2:3093923/93‑1 (MQ=255)
                   aGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTAc                                                                      <  1:4303010/93‑1 (MQ=255)
                   aGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTAc                                                                      <  2:5546474/93‑1 (MQ=255)
                   aGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTAc                                                                      >  1:970911/1‑93 (MQ=255)
                    gATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACt                                                                     <  2:7521644/93‑1 (MQ=255)
                      tGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTgg                                                                   >  2:1017503/1‑93 (MQ=255)
                       ggCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTggg                                                                  >  1:10526/1‑93 (MQ=255)
                        gCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGggc                                                                 <  2:812659/93‑1 (MQ=255)
                        gCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGggc                                                                 <  1:6009327/93‑1 (MQ=255)
                          tAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGggcgg                                                               >  2:6724984/1‑93 (MQ=255)
                           aaCCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGggcggc                                                              >  2:2645424/1‑93 (MQ=255)
                            aCCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCa                                                             >  2:6201517/1‑93 (MQ=255)
                              cAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCaaa                                                           >  1:3088278/1‑93 (MQ=255)
                                ggCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAAt                                                         >  1:2255081/1‑93 (MQ=255)
                                  cGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCg                                                       <  1:7213998/93‑1 (MQ=255)
                                   gAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGc                                                      >  2:6419966/1‑93 (MQ=255)
                                    aaCATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGcg                                                     >  1:1998993/1‑92 (MQ=255)
                                    aaCATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCt                                                     <  1:3766182/93‑1 (MQ=255)
                                    aaCATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCt                                                     >  1:7844434/1‑93 (MQ=255)
                                     aCATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTg                                                    >  1:4794439/1‑93 (MQ=255)
                                      cATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTgg                                                   >  2:7427053/1‑93 (MQ=255)
                                       aTCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTggg                                                  >  1:1416816/1‑92 (MQ=255)
                                       aTCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGt                                                  <  1:3524616/93‑1 (MQ=255)
                                       aTCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGt                                                  >  2:1860807/1‑93 (MQ=255)
                                        tCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTg                                                 >  1:5134252/1‑93 (MQ=255)
                                        tCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTg                                                 >  2:136276/1‑93 (MQ=255)
                                          ccGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGAt                                               >  2:2221956/1‑93 (MQ=255)
                                           cGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGAtt                                              >  2:2758942/1‑93 (MQ=255)
                                            ggTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTa                                             <  2:6064503/93‑1 (MQ=255)
                                             gTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTAc                                            >  2:3692614/1‑93 (MQ=255)
                                             gTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTAc                                            <  1:5947908/93‑1 (MQ=255)
                                              ttATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGATCTACAGTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACa                                           <  1:845143/93‑1 (MQ=255)
                                              ttATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACa                                           <  2:5412118/93‑1 (MQ=255)
                                              ttATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACa                                           >  2:4545123/1‑93 (MQ=255)
                                              ttATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACa                                           <  2:1403817/93‑1 (MQ=255)
                                                   aCTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATcgcg                                      >  1:1526418/1‑93 (MQ=255)
                                                   aCTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATcgcg                                      >  1:7346854/1‑93 (MQ=255)
                                                    ctctTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCga                                     >  2:7275040/1‑93 (MQ=255)
                                                     tctTGACCGCCAGGCAACGAAAGGTGAAGTTGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCgaa                                    <  1:6522328/93‑1 (MQ=255)
                                                     tctTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCgaa                                    >  1:4041112/1‑93 (MQ=255)
                                                         gACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGaaa                                >  2:1159425/1‑93 (MQ=255)
                                                         gACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGaaa                                <  1:4639887/93‑1 (MQ=255)
                                                           ccGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGc                              <  2:1122646/93‑1 (MQ=255)
                                                           ccGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGc                              <  2:3176245/93‑1 (MQ=255)
                                                           ccGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGc                              <  1:2084640/93‑1 (MQ=255)
                                                            cGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCg                             >  1:703060/1‑93 (MQ=255)
                                                              ccAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCggg                           >  2:5644936/1‑93 (MQ=255)
                                                              ccAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCggg                           <  2:3723388/93‑1 (MQ=255)
                                                              ccAGGCAACCAAAGGTGAAGTCGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCggg                           <  2:2507760/93‑1 (MQ=255)
                                                                aGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTg                         >  2:3865055/1‑93 (MQ=255)
                                                                aGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTg                         >  1:3656534/1‑93 (MQ=255)
                                                                aGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTg                         >  1:4540531/1‑93 (MQ=255)
                                                                 ggCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGa                        <  2:6163534/93‑1 (MQ=255)
                                                                  gCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGaa                       <  2:2912013/93‑1 (MQ=255)
                                                                     aCGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGt                    >  1:3033128/1‑93 (MQ=255)
                                                                     aCGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGt                    >  2:4901341/1‑93 (MQ=255)
                                                                     aCGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGt                    <  2:590060/93‑1 (MQ=255)
                                                                      cGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTg                   <  1:5831290/93‑1 (MQ=255)
                                                                       gAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGc                  >  2:5686979/1‑93 (MQ=255)
                                                                        aaaGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGcc                 >  1:1321281/1‑93 (MQ=255)
                                                                         aaGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCa                >  1:1752106/1‑93 (MQ=255)
                                                                           ggTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCaaa              >  2:3236174/1‑93 (MQ=255)
                                                                            gTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAg             <  1:3247602/93‑1 (MQ=255)
                                                                            gTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAg             <  2:6282598/93‑1 (MQ=255)
                                                                              gAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGtt           >  1:2630898/1‑93 (MQ=255)
                                                                              gAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGtt           >  1:5698148/1‑93 (MQ=255)
                                                                               aaGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTa          <  2:6809120/93‑1 (MQ=255)
                                                                               aaGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTa          <  2:213006/93‑1 (MQ=255)
                                                                               aaGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTa          >  1:5589019/1‑93 (MQ=255)
                                                                                 gtggtgAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATc        >  1:7428113/1‑93 (MQ=255)
                                                                                 atggtgAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATc        <  1:2576578/92‑1 (MQ=255)
                                                                                  tggtgAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCg       <  1:5421854/93‑1 (MQ=255)
                                                                                   ggtgAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGa      <  2:4120010/93‑1 (MQ=255)
                                                                                   ggtgAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGa      <  2:4408702/93‑1 (MQ=255)
                                                                                   ggtgAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGa      <  2:7815535/93‑1 (MQ=255)
                                                                                     tgAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGc    <  1:2941323/93‑1 (MQ=255)
                                                                                     tgAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGc    <  1:3439471/93‑1 (MQ=255)
                                                                                       aGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTg  <  2:4666722/93‑1 (MQ=255)
                                                                                          |                                                                                         
GCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTG  >  NC_000913/3934186‑3934365

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