Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,302,652 A→G Y513H (TAT→CAT)  yjcS ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,302,6520AG99.7% 1210.1 / NA 383Y513H (TAT→CAT) yjcShypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (178/204);  total (178/205)
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.99e-01

Only 100 of 390 total aligned reads displayed.

TGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTCAGCAGTTCTGCCG  >  NC_000913/4302562‑4302741
                                                                                          |                                                                                         
tggtgCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGc                                                                                         <  2:6609419/93‑1 (MQ=255)
  gtgCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGccc                                                                                       <  2:3906569/93‑1 (MQ=255)
   tgCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCa                                                                                      <  2:3065667/93‑1 (MQ=255)
     ccGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGc                                                                                    >  2:3527562/1‑93 (MQ=255)
     ccGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGc                                                                                    <  1:2552788/93‑1 (MQ=255)
     ccGTCGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGc                                                                                    <  1:3840675/93‑1 (MQ=255)
      cGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCt                                                                                   >  1:1647065/1‑93 (MQ=255)
       gTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTg                                                                                  >  1:7275652/1‑93 (MQ=255)
       gTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTg                                                                                  >  2:5447364/1‑93 (MQ=255)
        tGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGt                                                                                 <  2:294981/93‑1 (MQ=255)
        tGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGt                                                                                 >  1:2838958/1‑93 (MQ=255)
        tGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGt                                                                                 >  2:6059795/1‑93 (MQ=255)
          gCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGATGGTc                                                                               >  1:92986/1‑93 (MQ=255)
           cTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCa                                                                              >  1:7840995/1‑93 (MQ=255)
            tGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCaa                                                                             <  1:2692795/93‑1 (MQ=255)
             gAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGTCCCAGCTGTTCaaa                                                                            <  2:7167419/93‑1 (MQ=255)
             gAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCaaa                                                                            >  2:1209136/1‑93 (MQ=255)
             gAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCaaa                                                                            >  1:5491567/1‑93 (MQ=255)
               aCTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTGGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGt                                                                          <  1:2145973/93‑1 (MQ=255)
               aCTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGt                                                                          <  1:5758832/93‑1 (MQ=255)
               aCTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGAGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGt                                                                          >  2:1691945/1‑93 (MQ=255)
                 ttATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTa                                                                        <  1:7137740/93‑1 (MQ=255)
                 ttATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTa                                                                        <  2:1486127/93‑1 (MQ=255)
                 ttATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTa                                                                        <  2:2736298/93‑1 (MQ=255)
                  tATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTAt                                                                       <  1:3528592/93‑1 (MQ=255)
                    tGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTc                                                                     <  2:4784294/93‑1 (MQ=255)
                    tGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTc                                                                     >  2:4181343/1‑93 (MQ=255)
                       aCCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCACAGTTCTTCGCt                                                                  >  2:3209418/1‑93 (MQ=255)
                        ccccTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCtt                                                                 <  1:4002157/93‑1 (MQ=255)
                         cccTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTg                                                                >  1:1758832/1‑93 (MQ=255)
                           cTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCa                                                              <  1:3282297/93‑1 (MQ=255)
                           cTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCa                                                              >  1:2359868/1‑93 (MQ=255)
                           cTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCa                                                              <  1:7634264/93‑1 (MQ=255)
                            ttCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAg                                                             >  1:314288/1‑93 (MQ=255)
                             tCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGa                                                            <  1:830863/93‑1 (MQ=255)
                              cgcgCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGAt                                                           >  1:3023413/1‑93 (MQ=255)
                                cgcAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTTATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGAttc                                                         >  1:1284071/1‑93 (MQ=255)
                                cgcAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGAttc                                                         <  2:2638466/93‑1 (MQ=255)
                                cgcAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGAttc                                                         <  2:7199944/93‑1 (MQ=255)
                                cgcAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGAttc                                                         >  2:1171840/1‑93 (MQ=255)
                                 gcAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGAttct                                                        <  1:3474015/93‑1 (MQ=255)
                                 gcAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGAttct                                                        >  1:6553674/1‑93 (MQ=255)
                                  cAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGAttctt                                                       <  2:7520483/93‑1 (MQ=255)
                                  cAGCTCTTTCGCGCCGGTCAGGTAGAAACCCCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGAttctt                                                       <  1:6721001/93‑1 (MQ=255)
                                  cAGCGCTTTCGCGCCGGTCAGGTAGAACCCGCGCTATGTGGCGCATTCCGGCTGATTGCCTAGCTGTTCAGAGTTATTCGCTTGCAGAttctt                                                       <  2:3656821/93‑1 (MQ=255)
                                   aGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGAttcttc                                                      <  2:789852/93‑1 (MQ=255)
                                   aGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGAttcttc                                                      >  2:19004/1‑93 (MQ=255)
                                   aGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGAttcttc                                                      <  2:6949577/93‑1 (MQ=255)
                                    gttctTTCGCGCCGGTCAGTTAGAAACCGTGCTATGTGGCGGACTCGTCCTGATGGCCCAGCTGTTCAAAGTTATTCTCGTGCAGATTCTTcg                                                     >  1:7018431/3‑93 (MQ=255)
                                    gCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTcg                                                     >  2:2828047/1‑93 (MQ=255)
                                    gCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTcg                                                     >  1:7483242/1‑93 (MQ=255)
                                    gCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTcg                                                     <  1:4811018/93‑1 (MQ=255)
                                    gCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTcg                                                     <  2:5537796/93‑1 (MQ=255)
                                     ctctTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGTCCTGATGTCTCAGCTGTTCAAAGGTATTCGCTTGCAGATTCTTcgc                                                    <  2:530513/93‑1 (MQ=255)
                                     ctctTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTcgc                                                    >  1:1162797/1‑93 (MQ=255)
                                       ctTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGa                                                  >  2:4991676/1‑93 (MQ=255)
                                        tttCGCGCCGGTCAGGTATAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGAc                                                 >  2:2808636/1‑93 (MQ=255)
                                         ttCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGAcc                                                <  1:6267113/93‑1 (MQ=255)
                                          tCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCt                                               <  2:1846161/93‑1 (MQ=255)
                                           cgcgCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTg                                              <  2:5169904/93‑1 (MQ=255)
                                           cgcgCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTg                                              <  2:5566960/93‑1 (MQ=255)
                                             cgcCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGt                                            >  2:5702646/1‑93 (MQ=255)
                                               ccGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCa                                          >  2:1399626/1‑93 (MQ=255)
                                                cGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCAc                                         >  1:3859918/1‑93 (MQ=255)
                                                 ggTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTcccc                                        <  2:6835904/93‑1 (MQ=255)
                                                 ggTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCAcc                                        >  2:7762413/1‑93 (MQ=255)
                                                  gTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCAccc                                       >  1:2063780/1‑93 (MQ=255)
                                                   tCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCg                                      >  1:5297084/1‑93 (MQ=255)
                                                   tCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCg                                      >  1:6141716/1‑93 (MQ=255)
                                                    cAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCgg                                     >  2:3465811/1‑93 (MQ=255)
                                                      ggTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTGAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGt                                   >  1:4113531/1‑93 (MQ=255)
                                                       gTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGtt                                  <  2:6732932/93‑1 (MQ=255)
                                                        tAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTg                                 >  1:1049414/1‑93 (MQ=255)
                                                        tAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTg                                 >  1:1852475/1‑93 (MQ=255)
                                                        tAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTg                                 <  1:4976344/93‑1 (MQ=255)
                                                         aGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTgg                                >  2:2117712/1‑93 (MQ=255)
                                                         aGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTgg                                <  2:4263810/93‑1 (MQ=255)
                                                          gAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTggc                               >  2:7310859/1‑93 (MQ=255)
                                                          gAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTggc                               <  2:5912681/93‑1 (MQ=255)
                                                              ccGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCa                           >  2:1118216/1‑93 (MQ=255)
                                                               cgcgCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCaa                          <  2:3399433/93‑1 (MQ=255)
                                                                gcgcCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAAt                         >  1:3646443/1‑93 (MQ=255)
                                                                gcgcCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAAt                         <  1:4394464/93‑1 (MQ=255)
                                                                 cgcCATGTGGCGGACTCGTCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATc                        <  2:2378673/93‑1 (MQ=255)
                                                                  gcCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCa                       >  1:137381/1‑93 (MQ=255)
                                                                   ccATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCAc                      <  2:4634582/93‑1 (MQ=255)
                                                                    cATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCAcc                     <  2:7432169/93‑1 (MQ=255)
                                                                     aTGTGGCGGACTAGGCCTGATGGCCCAGCTGTTCAAAGTTATTAGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCt                    <  2:7016128/93‑1 (MQ=255)
                                                                      tgtgGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGCGTTGGCGGCAATCACCTg                   >  2:4388402/1‑93 (MQ=255)
                                                                       gtgGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGt                  <  1:4908258/93‑1 (MQ=255)
                                                                           cGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTca              <  2:7709019/93‑1 (MQ=255)
                                                                             gACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTcagc            >  1:4703440/1‑93 (MQ=255)
                                                                              aCTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTcagca           >  2:2541938/1‑93 (MQ=255)
                                                                              aCTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTcagca           >  1:2959935/1‑93 (MQ=255)
                                                                              aCTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTcagca           <  1:5881008/93‑1 (MQ=255)
                                                                              aCTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTcagca           <  2:4223679/93‑1 (MQ=255)
                                                                                tCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTCAGCAGt         >  2:6335283/1‑93 (MQ=255)
                                                                                 cGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTCAGCAGtt        <  1:4606076/93‑1 (MQ=255)
                                                                                 cGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTCAGCAGtt        <  2:5848533/93‑1 (MQ=255)
                                                                                  ggCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTCAGCAGTTc       <  1:5034586/93‑1 (MQ=255)
                                                                                          |                                                                                         
TGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTCAGCAGTTCTGCCG  >  NC_000913/4302562‑4302741

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