Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,342,081 A→G Y37H (TAT→CAT)  yfaL ← putative ATP‑binding component of a transport system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,342,0810AG90.2% 592.6 / 45.6 204Y37H (TAT→CAT) yfaLputative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand):  ref base A (10/10);  new base G (74/110);  total (84/120)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.69e-01

Only 100 of 208 total aligned reads displayed.

AAACCAGCCATTGCCCATCAGCGATACTCCAGTCCTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATATCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGAAAAGAGATGCAATCATTGACGGGAGTAAAGATAAATACTCCTTGC  >  NC_000913/2341991‑2342170
                                                                                          |                                                                                         
aaaCCAGCCATTGCCCATCAGCGATACTCCAGTCCTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTc                                                                                         >  1:3194318/1‑93 (MQ=255)
aaaCCAGCCATTGCCCATCAGCGATACTCCAGTCCTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATATc                                                                                         >  1:4197172/1‑93 (MQ=255)
  aCCAGCCATTGCCCATCAGCGATACTCCAGTCCTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTccc                                                                                       <  2:2722843/93‑1 (MQ=255)
  aCCAGCCATTGCCCATCAGCGATACTCCAGTCCTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTccc                                                                                       >  1:998609/1‑93 (MQ=255)
  aCCAGCCATTGCCCATCAGCGATACTCCAGTCCTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATATccc                                                                                       <  1:122379/93‑1 (MQ=255)
   ccAGCCATTGCCCATCAGCGATACTCCAGTCCTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCt                                                                                      >  1:3231472/1‑93 (MQ=255)
   ccAGCCATTGCCCATCAGCGATACTCCAGTCCTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCt                                                                                      <  1:1934657/93‑1 (MQ=255)
    cAGCCATTGCCCATCAGCGATACTCCAGTCCTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATATCCCTg                                                                                     >  2:6402724/1‑93 (MQ=255)
          ttGCCCATCAGCGATACTCCAGTCCTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGa                                                                               <  1:4498283/93‑1 (MQ=255)
             cccATCAGCGATACTCCAGTCCTGCGTAATGCATGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATc                                                                            >  1:6319558/1‑93 (MQ=255)
              cgaTCAGCGATGCGCCAGTCCTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCa                                                                           <  2:5800170/91‑1 (MQ=255)
              ccATCAGCGATACTCCAGTCCTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCa                                                                           >  2:4889843/1‑93 (MQ=255)
               cATCAGCGATACTCCAGTCCTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAg                                                                          >  2:4900273/1‑93 (MQ=255)
               cATCAGCGATACTCCAGTCCTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAg                                                                          <  1:2931026/93‑1 (MQ=255)
                aTCAGCGATACTCCAGTCCTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGt                                                                         <  1:1832941/93‑1 (MQ=255)
                aTCAGCGATACTCCAGTCCTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGt                                                                         <  1:3341992/93‑1 (MQ=255)
                aTCAGCGATACTCCAGTCCTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATATCCCTGGCATGAATCAGt                                                                         <  1:3497097/93‑1 (MQ=255)
                 tCAGCGATACTCCAGTCCTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTg                                                                        >  2:318885/1‑93 (MQ=255)
                 tCAGCGATACTCCAGTCCTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATATCCCTGGCATGAATCAGTg                                                                        >  1:454335/1‑93 (MQ=255)
                  cAGCGATACTCCAGTCCTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGa                                                                       >  2:143887/1‑93 (MQ=255)
                  cAGCGATACTCCAGTCCTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGa                                                                       <  2:2246041/93‑1 (MQ=255)
                    gCGATACTCCAGTCCTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATATCCCTGGCATGAATCAGTGAcc                                                                     <  1:1474378/93‑1 (MQ=255)
                    gCGATACTCCAGTCCTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATATCCCTGGCATGAATCAGTGAcc                                                                     >  2:3935473/1‑93 (MQ=255)
                      gATACTCCAGTCCTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGAccgc                                                                   >  1:6238697/1‑93 (MQ=255)
                       aTACTCCAGTCCTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGAccgcc                                                                  >  1:5489598/1‑93 (MQ=255)
                       aTACTCCAGTCCTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATATCCCTGGCATGAATCAGTGAccgcc                                                                  >  1:1169337/1‑93 (MQ=255)
                        tACTCCAGTCCTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGAccgccg                                                                 <  2:6002012/93‑1 (MQ=255)
                          cTCCAGTCCTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCcgcg                                                               >  2:557439/1‑93 (MQ=255)
                          cTCCAGTCCTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCcgcg                                                               >  1:1600323/1‑93 (MQ=255)
                           tCCAGTCCTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGa                                                              <  2:2024017/93‑1 (MQ=255)
                           tCCAGTCCTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGa                                                              <  1:2255724/93‑1 (MQ=255)
                           tCCAGTCCTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGa                                                              <  2:3518814/93‑1 (MQ=255)
                            ccAGTCCTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACCACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGac                                                             >  1:3167816/1‑93 (MQ=255)
                             cAGTCCTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGaca                                                            >  1:3671808/1‑93 (MQ=255)
                             cAGTCCTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGaca                                                            >  2:1831312/1‑93 (MQ=255)
                              aGTCCTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGacac                                                           <  1:4817642/93‑1 (MQ=255)
                              aGTCCTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGacac                                                           <  2:4081377/93‑1 (MQ=255)
                               gTCCTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCAAGTCCCTGGCATGAATCAGTGACCGCCGCGACAcc                                                          <  2:2642197/93‑1 (MQ=255)
                                tCCTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCg                                                         >  1:2824355/1‑93 (MQ=255)
                                 ccTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGt                                                        <  2:3642352/93‑1 (MQ=255)
                                 ccTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGt                                                        <  1:1554866/93‑1 (MQ=255)
                                 ccTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATATCCCTGGCATGAATCAGTGACCGCCGCGACACCGt                                                        >  1:1656144/1‑93 (MQ=255)
                                  cTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCggt                                                       >  2:6317709/1‑91 (MQ=255)
                                  cTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGtt                                                       >  1:6221688/1‑93 (MQ=255)
                                  cTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGtt                                                       <  2:5702268/93‑1 (MQ=255)
                                   tGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTa                                                      <  1:5899422/93‑1 (MQ=255)
                                   tGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTa                                                      <  1:3526561/93‑1 (MQ=255)
                                   tGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTa                                                      <  1:1415096/93‑1 (MQ=255)
                                      gTAATGCCTGAAAGGCTTTGCCTGCTGGCCTTACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCa                                                   <  1:6477398/93‑1 (MQ=255)
                                      gTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCa                                                   >  1:5601569/1‑93 (MQ=255)
                                      gTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCa                                                   <  2:2971288/93‑1 (MQ=255)
                                        aaTGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCaga                                                 <  2:239361/93‑1 (MQ=255)
                                        aaTGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTCACCGCCGCGACACCGTTCGCagg                                                 >  2:1302621/1‑92 (MQ=255)
                                         aTGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCagag                                                <  1:1393223/93‑1 (MQ=255)
                                         aTGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCagag                                                >  2:1122567/1‑93 (MQ=255)
                                          tGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAgaga                                               >  1:1698167/1‑93 (MQ=255)
                                           gCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCTGCGACACCGTTAGCAGAGaa                                              <  1:3774877/93‑1 (MQ=255)
                                            ccTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGaaa                                             >  2:1505016/1‑93 (MQ=255)
                                              tGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGAAAag                                           >  1:2985215/1‑93 (MQ=255)
                                              tGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGAAAag                                           <  1:1484595/93‑1 (MQ=255)
                                              tGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATATCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGAAAag                                           >  2:1237783/1‑93 (MQ=255)
                                               gAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATATCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGAAAaga                                          <  1:874681/93‑1 (MQ=255)
                                                 aaGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGAAAAgaga                                        <  2:3047269/93‑1 (MQ=255)
                                                      tttGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGAAAAGAGATGCaa                                   <  1:2849079/93‑1 (MQ=255)
                                                      tttGCCTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGAAAAGAGATGCaa                                   <  2:6057445/93‑1 (MQ=255)
                                                          ccTGCTGGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGAAAAGAGATGCAATCAt                               <  2:1089288/93‑1 (MQ=255)
                                                           ctgctgGCCTGACACCTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGAAAAGAGATGCAATCAtt                              <  1:5012215/93‑1 (MQ=255)
                                                             gctgGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGAAAAGAGATGCAATCATTGa                            >  1:816301/1‑93 (MQ=255)
                                                             gctgGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGAAAAGAGATGCAATCATTGa                            <  2:4688634/93‑1 (MQ=255)
                                                             gctgGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGAAAAGAGATGCAATCATTGa                            <  1:5576911/93‑1 (MQ=255)
                                                              ctgGCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGAAAAGAGATGCAATCATTGAc                           <  1:2646157/93‑1 (MQ=255)
                                                                ggCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGAAAAGAGATGCAATCATTGACgg                         <  1:5546172/93‑1 (MQ=255)
                                                                ggCCTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGAAAAGAGATGCAATCATTGACgg                         <  1:5824332/93‑1 (MQ=255)
                                                                 gcgtGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGAAAAGAGATGCAATCATTGACggg                        <  2:1054575/90‑1 (MQ=255)
                                                                  ccTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGAAAAGAGATGCAATCATTGACGGGa                       <  2:1271856/93‑1 (MQ=255)
                                                                  ccTGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGAAAAGAGATGCAATCATTGACGGGa                       <  2:4648118/93‑1 (MQ=255)
                                                                    tGACAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGAAAAGAGATGCAATCATTGACGGGAGt                     <  1:3023657/93‑1 (MQ=255)
                                                                     gACAACTCGCTTTGACATCATATCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGAAAAGAGATGCAATCATTGACGGGAGTa                    <  2:1625244/93‑1 (MQ=255)
                                                                       cAACTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGAAAAGAGATGCAATCATTGACGGGAGTaaa                  >  1:5229472/1‑93 (MQ=255)
                                                                         aCTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGAAAAGAGATGCAATCATTGACGGGAGTAAAGa                >  2:1211579/1‑93 (MQ=255)
                                                                         aCTCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGAAAAGAGATGCAATCATTGACGGGAGTAAAGa                >  2:6144011/1‑93 (MQ=255)
                                                                           tCGCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGAAAAGAGATGCAATCATTGACGGGAGTAAAGATa              >  1:4241082/1‑93 (MQ=255)
                                                                             gCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGAAAAGAGATGCAATCATTGACGGGAGTAAAGATaaa            <  2:2282525/93‑1 (MQ=255)
                                                                             gCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGAAAAGAGATGCAATCATTGACGGGAGTAAAGATaaa            >  2:5366946/1‑93 (MQ=255)
                                                                             gCTTTGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGAAAAGAGATGCAATCATTGACGGGAGTAAAGATaaa            <  2:4271834/93‑1 (MQ=255)
                                                                               tttGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGAAAAGAGATGCAATCATTGACGGGAGTAAAGATAAATa          >  1:4052256/1‑93 (MQ=255)
                                                                                ttGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGAAAAGAGATGCAATCATTGACGGGAGTAAAGATAAATAc         <  2:5296079/93‑1 (MQ=255)
                                                                                ttGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGAAAAGAGATGCAATCATTGACGGGAGTAAAGATAAATAc         >  2:5919697/1‑93 (MQ=255)
                                                                                ttGACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGAAAAGAGATGCAATCATTGACGGGAGTAAAGATAAATAc         <  2:595708/93‑1 (MQ=255)
                                                                                  gACATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGAAAAGAGATGCAATCATTGACGGGAGTAAAGATAAATACTc       <  2:3507083/93‑1 (MQ=255)
                                                                                   aCATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGAAAAGAGATGCAATCATTGACGGGAGTAAAGATAAATACTcc      <  1:3704608/93‑1 (MQ=255)
                                                                                   aCATCATGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGAAAAGAGATGCAATCATTGACGGGAGTAAAGATAAATACTcc      >  1:6591352/1‑93 (MQ=255)
                                                                                    catcatGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGAAAAGAGATGCAATCATTGACGGGAGTAAAGATAAATACTCCt     <  2:2827676/93‑1 (MQ=255)
                                                                                    catcatGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGAAAAGAGATGCAATCATTGACGGGAGTAAAGATAAATACTCCt     <  1:3942521/93‑1 (MQ=255)
                                                                                    catcatGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGAAAAGAGATGCAATCATTGACGGGAGTAAAGATAAATACTCCt     >  1:5281620/1‑93 (MQ=255)
                                                                                      tcatGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGAAAAGAGATGCAATCATTGACGGGAGTAAAGATAAATACTCCTTg   <  2:395921/93‑1 (MQ=255)
                                                                                      tcatGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGAAAAGAGATGCAATCATTGACGGGAGTAAAGATAAATACTCCTTg   <  1:544052/93‑1 (MQ=255)
                                                                                      tcatGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGAAAAGAGATGCAATCATTGACGGGAGCAAAGATAAATACTCCTTg   >  1:4339340/1‑93 (MQ=255)
                                                                                       catGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGAAAAGAGATGCAATCATTGACGGGAGTAAAGATAAATACTCCTTGc  >  2:3378838/1‑93 (MQ=255)
                                                                                       catGTCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGAAAAGAGATGCAATCATTGACGGGAGTAAAGATAAATACTCCTTGc  <  2:3311154/93‑1 (MQ=255)
                                                                                          |                                                                                         
AAACCAGCCATTGCCCATCAGCGATACTCCAGTCCTGCGTAATGCCTGAAAGGCTTTGCCTGCTGGCCTGACAACTCGCTTTGACATCATATCCCTGGCATGAATCAGTGACCGCCGCGACACCGTTAGCAGAGAAAAGAGATGCAATCATTGACGGGAGTAAAGATAAATACTCCTTGC  >  NC_000913/2341991‑2342170

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