Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,259,919 C→T R156H (CGC→CAC)  fruK ← fructose‑1‑phosphate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,259,9190CT99.2% 444.8 / ‑5.0 133R156H (CGC→CAC) fruKfructose‑1‑phosphate kinase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (55/76);  total (56/77)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01

Only 100 of 136 total aligned reads displayed.

AGGTTTCACCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGCGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGCTC  >  NC_000913/2259828‑2260008
                                                                                           |                                                                                         
aGGTTTCACCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGCg                                                                                          <  2:3134386/93‑1 (MQ=255)
 ggTTTCACCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGa                                                                                         <  2:3457607/93‑1 (MQ=255)
  gTTTCACCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAg                                                                                        >  2:2609099/1‑93 (MQ=255)
  gTTTCACCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAg                                                                                        <  2:4154708/93‑1 (MQ=255)
   tttCACCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGt                                                                                       >  2:1479530/1‑93 (MQ=255)
   tttCACCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGt                                                                                       <  2:4315806/93‑1 (MQ=255)
      cACCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCac                                                                                    >  1:264120/1‑92 (MQ=255)
      cACCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCAt                                                                                    >  1:3239304/1‑93 (MQ=255)
      cACCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCAt                                                                                    >  2:2592016/1‑93 (MQ=255)
       aCCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATc                                                                                   >  2:1483687/1‑93 (MQ=255)
       aCCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATc                                                                                   <  1:260165/93‑1 (MQ=255)
        ccAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATcc                                                                                  <  2:3874746/93‑1 (MQ=255)
         cAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCa                                                                                 >  2:4146699/1‑93 (MQ=255)
          aGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAg                                                                                >  1:465505/1‑93 (MQ=255)
           gCCACGGTGCCGCTTTCAAGCCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGt                                                                               <  1:1506103/93‑1 (MQ=255)
           gCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGt                                                                               >  1:3754222/1‑93 (MQ=255)
             cACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCg                                                                             <  2:2317394/93‑1 (MQ=255)
              aCGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCgg                                                                            >  1:4568890/1‑93 (MQ=255)
              aCGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCgg                                                                            >  2:1353327/1‑93 (MQ=255)
              aCGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCgg                                                                            >  1:465162/1‑93 (MQ=255)
              aCGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCgg                                                                            <  1:3054267/93‑1 (MQ=255)
                  taccGCTGTCAAACCTGCTACTAACGATTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACACAGGTGAGTCATCCAGTCGGTGaa                                                                        <  1:1576462/91‑1 (MQ=255)
                   gCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAAc                                                                       <  2:4457904/93‑1 (MQ=255)
                      tctatCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCt                                                                    <  1:2359158/89‑1 (MQ=255)
                        tttCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTc                                                                  <  1:3462302/93‑1 (MQ=255)
                         ttCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTcc                                                                 <  2:1442348/93‑1 (MQ=255)
                           cAACCCTGCTACTAACGCTTCACGTCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACCCTTCCgg                                                               <  2:3839854/93‑1 (MQ=255)
                              aCCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCt                                                            <  1:2988381/93‑1 (MQ=255)
                              aCCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCt                                                            <  2:4733022/93‑1 (MQ=255)
                                cTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGa                                                          <  2:2456109/93‑1 (MQ=255)
                                 tGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGAc                                                         >  2:771536/1‑93 (MQ=255)
                                 tGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGAc                                                         <  1:3687182/93‑1 (MQ=255)
                                   ctactaACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGc                                                       <  1:819665/93‑1 (MQ=255)
                                    tactacCGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGcc                                                      <  2:3328252/93‑1 (MQ=255)
                                    tactaACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGcc                                                      <  2:2564581/93‑1 (MQ=255)
                                    tactaACGCTTCACGGCTAATATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGATGAACGCTTCCGGGCTGACGcc                                                      <  1:3013085/93‑1 (MQ=255)
                                     actaACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCt                                                     <  2:3472755/93‑1 (MQ=255)
                                     actaACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCt                                                     >  1:1456106/1‑93 (MQ=255)
                                      ctaACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCta                                                    >  2:51589/1‑92 (MQ=255)
                                      ctaACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTg                                                    >  1:2953354/1‑93 (MQ=255)
                                       taACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGa                                                   <  1:3359711/93‑1 (MQ=255)
                                       taACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGa                                                   <  2:2000403/93‑1 (MQ=255)
                                       taACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGa                                                   <  1:4628251/93‑1 (MQ=255)
                                        aaCGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGAc                                                  <  1:3220255/93‑1 (MQ=255)
                                        aaCGCTTCACGGCTAATATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGAc                                                  <  1:1554139/93‑1 (MQ=255)
                                          cGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACgg                                                <  1:548349/93‑1 (MQ=255)
                                          cGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACgg                                                <  2:85558/93‑1 (MQ=255)
                                           gCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGt                                               >  1:1260320/1‑93 (MQ=255)
                                            cTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTa                                              <  1:4561102/93‑1 (MQ=255)
                                             ttCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTaa                                             <  2:4321250/93‑1 (MQ=255)
                                               cACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGc                                           <  2:1520222/93‑1 (MQ=255)
                                                aCGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCt                                          >  2:3707080/1‑93 (MQ=255)
                                                aCGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCt                                          >  2:171375/1‑93 (MQ=255)
                                                  ggCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTc                                        <  1:653989/93‑1 (MQ=255)
                                                  ggCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTc                                        <  1:4485136/93‑1 (MQ=255)
                                                  ggCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTc                                        >  1:3931274/1‑93 (MQ=255)
                                                  ggCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTc                                        >  1:400650/1‑93 (MQ=255)
                                                  ggCTACTAGCAAAGATAATGAAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTc                                        <  2:866719/93‑1 (MQ=255)
                                                   gCTACTATCAAAGATAATGCAAGGACACTGGCGACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTcc                                       <  1:1493808/93‑1 (MQ=255)
                                                   gCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTcc                                       <  1:86045/93‑1 (MQ=255)
                                                     tactaTCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGc                                     <  2:2809193/93‑1 (MQ=255)
                                                     tactaTCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGc                                     <  1:1957532/93‑1 (MQ=255)
                                                     tactaTCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGc                                     <  2:4671610/93‑1 (MQ=255)
                                                      actaTCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCt                                    >  1:1990285/1‑93 (MQ=255)
                                                      actaTCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCt                                    >  1:249184/1‑93 (MQ=255)
                                                        taTCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGa                                  >  2:2523999/1‑93 (MQ=255)
                                                        taTCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGa                                  >  1:3678773/1‑93 (MQ=255)
                                                        taTCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGa                                  <  2:2952452/93‑1 (MQ=255)
                                                         agcAAAGAAAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGGCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGac                                 <  2:1918426/91‑1 (MQ=255)
                                                         aTCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGac                                 <  2:3614054/93‑1 (MQ=255)
                                                         aTCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGac                                 >  2:920411/1‑93 (MQ=255)
                                                         aTCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGac                                 >  1:3463105/1‑93 (MQ=255)
                                                          tCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGaca                                >  1:4056893/1‑93 (MQ=255)
                                                          tCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGaca                                <  2:1969298/93‑1 (MQ=255)
                                                          tCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGaca                                >  1:4704907/1‑93 (MQ=255)
                                                          tCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGaca                                >  1:389437/1‑93 (MQ=255)
                                                          tCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGaca                                >  2:2796862/1‑93 (MQ=255)
                                                           cAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGacac                               >  1:1736811/1‑93 (MQ=255)
                                                               gATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGAc                           <  2:1620396/93‑1 (MQ=255)
                                                                aTAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGAcc                          <  2:231824/93‑1 (MQ=255)
                                                                 tAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCa                         >  1:3450908/1‑93 (MQ=255)
                                                                  aaTGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCat                        <  1:162650/93‑1 (MQ=255)
                                                                  aaTGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCat                        >  2:3806324/1‑93 (MQ=255)
                                                                   aTGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCata                       <  1:1803713/93‑1 (MQ=255)
                                                                   aTGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCata                       <  1:1803711/93‑1 (MQ=255)
                                                                        tggACAATGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGaa                  >  1:2396127/2‑93 (MQ=255)
                                                                         ggACACTGACTGCGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAAc                 <  1:3916742/93‑1 (MQ=255)
                                                                         ggACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAAc                 >  1:3029991/1‑93 (MQ=255)
                                                                         ggACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAAc                 <  2:3889318/93‑1 (MQ=255)
                                                                            cacTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGa              >  2:916880/1‑93 (MQ=255)
                                                                             acTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGAc             <  1:61430/93‑1 (MQ=255)
                                                                               tGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACGGACCATATCGAACTGACCg           >  2:3883294/1‑93 (MQ=255)
                                                                                   tACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGcc       <  2:2427559/93‑1 (MQ=255)
                                                                                   tACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGcc       <  1:462689/93‑1 (MQ=255)
                                                                                   tACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGcc       >  1:2765149/1‑93 (MQ=255)
                                                                                      gCAGGAGAGTCCTCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAgt    >  2:3238668/1‑92 (MQ=255)
                                                                                       cAGGTGAGTCATCCAGTGGGTGAACGCTTCCGGGCTGACGCCTGAGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGCTc  <  1:114887/93‑1 (MQ=255)
                                                                                        aGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGCTc  >  2:3595969/1‑93 (MQ=255)
                                                                                        aGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGCTc  >  1:1413019/1‑93 (MQ=255)
                                                                                        aGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGCTc  >  2:3595968/1‑93 (MQ=255)
                                                                                           |                                                                                         
AGGTTTCACCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGCGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGCTC  >  NC_000913/2259828‑2260008

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