Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,621,644 A→G W15R (TGG→CGG)  lplA ← lipoate‑protein ligase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,621,6440AG99.5% 628.9 / ‑5.5 202W15R (TGG→CGG) lplAlipoate‑protein ligase A
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (87/114);  total (87/115)
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.92e-01

Only 100 of 207 total aligned reads displayed.

CTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTT  >  NC_000913/4621554‑4621733
                                                                                          |                                                                                         
ctacCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACcgcg                                                                                         <  2:2220145/93‑1 (MQ=255)
 tacCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCgg                                                                                        >  2:3514508/1‑93 (MQ=255)
  acCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCggg                                                                                       >  2:2913289/1‑93 (MQ=255)
       gtCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTa                                                                                  >  2:2312860/1‑93 (MQ=255)
         cGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAag                                                                                >  1:1389927/1‑93 (MQ=255)
             aTTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAagagca                                                                            >  2:2444164/1‑91 (MQ=255)
             aTTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTc                                                                            >  1:1034715/1‑93 (MQ=255)
             aTTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTc                                                                            <  2:2696915/93‑1 (MQ=255)
              ttGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCa                                                                           <  2:1180323/93‑1 (MQ=255)
              ttGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCa                                                                           >  2:4550216/1‑93 (MQ=255)
              ttGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCa                                                                           >  1:252839/1‑93 (MQ=255)
              ttGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGCCa                                                                           >  2:3681154/1‑93 (MQ=255)
                 cgcCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCagag                                                                        <  1:1371714/93‑1 (MQ=255)
                  gcCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAgaga                                                                       <  1:4341539/93‑1 (MQ=255)
                  gcCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGAGAAAAATACCCTCTTCAACCGCCAGGTTAACCCGCGGGTCGTAAGAGTCAgaga                                                                       <  2:4084632/93‑1 (MQ=255)
                   ccAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGAt                                                                      >  2:530345/1‑93 (MQ=255)
                   ccAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGAt                                                                      >  1:2777614/1‑93 (MQ=255)
                    cAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATg                                                                     <  1:1830103/93‑1 (MQ=255)
                    cAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATg                                                                     >  1:2665130/1‑93 (MQ=255)
                     agagAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGa                                                                    >  1:1300055/1‑93 (MQ=255)
                     agagAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGa                                                                    >  2:234620/1‑93 (MQ=255)
                     agagAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGa                                                                    <  2:4013513/93‑1 (MQ=255)
                     agagAAACAGAACGCGCGGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGa                                                                    <  1:3489401/93‑1 (MQ=255)
                      gagaAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAg                                                                   <  1:2396382/93‑1 (MQ=255)
                        gaAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCa                                                                 >  1:3888383/1‑93 (MQ=255)
                        gaAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCa                                                                 <  1:4164839/93‑1 (MQ=255)
                         aaaCAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAg                                                                >  2:219879/1‑93 (MQ=255)
                         aaaCAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAg                                                                <  2:3371303/93‑1 (MQ=255)
                          aaCAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAgg                                                               <  2:3287850/93‑1 (MQ=255)
                           aCAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGc                                                              <  1:4111327/93‑1 (MQ=255)
                            cAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCg                                                             <  1:120610/93‑1 (MQ=255)
                             aGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGt                                                            <  1:3096098/93‑1 (MQ=255)
                              gAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTa                                                           >  2:1947542/1‑93 (MQ=255)
                               aaCGCGCTGCGTGGCGGGCATTTGGCTAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTaa                                                          >  1:4003522/1‑93 (MQ=255)
                               aaCGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTaa                                                          >  1:3155136/1‑93 (MQ=255)
                                aCGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAACCCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAAt                                                         <  2:2519600/93‑1 (MQ=255)
                                aCGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAAt                                                         >  2:4582418/1‑93 (MQ=255)
                                 cgcgCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAAtg                                                        <  1:3032292/93‑1 (MQ=255)
                                  gcgcTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAAtgt                                                       >  2:2399486/1‑93 (MQ=255)
                                  gcgcTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAAtgt                                                       <  2:1701005/93‑1 (MQ=255)
                                  gcgcTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAAtgt                                                       <  1:2880461/93‑1 (MQ=255)
                                  gcgcTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAAtgt                                                       >  2:3993772/1‑93 (MQ=255)
                                   cgcTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAAtgtg                                                      >  2:267902/1‑93 (MQ=255)
                                   cgcTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAAtgtg                                                      <  2:360873/93‑1 (MQ=255)
                                   cgcTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAAtgtg                                                      <  1:1527380/93‑1 (MQ=255)
                                    gcTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTgg                                                     <  1:3719575/93‑1 (MQ=255)
                                     cTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGa                                                    <  2:669636/93‑1 (MQ=255)
                                     cTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGa                                                    >  1:3679720/1‑93 (MQ=255)
                                     cTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGa                                                    >  2:3413130/1‑93 (MQ=255)
                                      tGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGAc                                                   <  2:4621500/93‑1 (MQ=255)
                                         gTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGGCGTAAGAGTCAGAGTTGAGCGGGCGTAATGTGGACATa                                                >  1:1141473/1‑93 (MQ=255)
                                          tGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATaa                                               <  1:541078/93‑1 (MQ=255)
                                             cGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGa                                            <  2:4182705/93‑1 (MQ=255)
                                              gggCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGAt                                           >  1:57143/1‑93 (MQ=255)
                                               ggCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGAtt                                          >  2:3266358/1‑93 (MQ=255)
                                                gCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGAttt                                         <  2:2588259/93‑1 (MQ=255)
                                                 cATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTc                                        >  1:848499/1‑93 (MQ=255)
                                                 cATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTc                                        <  1:742961/93‑1 (MQ=255)
                                                  aTTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTcc                                       <  2:1250364/93‑1 (MQ=255)
                                                  aTTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTcc                                       <  1:2029432/93‑1 (MQ=255)
                                                   tttGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCt                                      >  1:3930467/1‑93 (MQ=255)
                                                    ttGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCtt                                     >  1:4573318/1‑93 (MQ=255)
                                                    ttGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCtt                                     <  1:99496/93‑1 (MQ=255)
                                                      ggCGAAAACTACACTCTTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTc                                   <  1:1471/93‑1 (MQ=255)
                                                      ggCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTc                                   <  2:567542/93‑1 (MQ=255)
                                                      ggCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTc                                   <  1:1346448/93‑1 (MQ=255)
                                                       gCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCt                                  <  1:917586/93‑1 (MQ=255)
                                                        cGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCtt                                 >  1:2157146/1‑93 (MQ=255)
                                                        cGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCtt                                 >  2:4476858/1‑93 (MQ=255)
                                                        cGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCtt                                 <  1:221351/93‑1 (MQ=255)
                                                        cGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCtt                                 <  1:1259177/93‑1 (MQ=255)
                                                           aaaaTACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGCGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTgta                              >  1:1411012/1‑93 (MQ=255)
                                                              aTACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTgtagta                           >  2:1497484/1‑93 (MQ=255)
                                                               tACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTaa                          >  1:2446641/1‑93 (MQ=255)
                                                               tACACTCTTCCACCGCCAGGCTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTaa                          >  2:2839988/1‑93 (MQ=255)
                                                                acacTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAAt                         <  2:1398004/93‑1 (MQ=255)
                                                                 cacTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATg                        <  2:404311/93‑1 (MQ=255)
                                                                  acTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATgg                       >  2:4167044/1‑93 (MQ=255)
                                                                   ctctTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATggg                      >  2:787524/1‑93 (MQ=255)
                                                                       tCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCAc                  <  1:1737523/93‑1 (MQ=255)
                                                                       tCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCAc                  >  2:2961199/1‑93 (MQ=255)
                                                                         cACGGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCActc                <  1:490160/93‑1 (MQ=255)
                                                                          aCCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCActct               <  2:3773993/93‑1 (MQ=255)
                                                                             gCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCtttt            >  2:2268281/1‑93 (MQ=255)
                                                                             gCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCtttt            <  2:3391837/93‑1 (MQ=255)
                                                                             gCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCtttt            <  2:1114508/93‑1 (MQ=255)
                                                                               cAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCtttttt          >  1:3599260/1‑93 (MQ=255)
                                                                               cAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCtttttt          >  2:3337015/1‑93 (MQ=255)
                                                                                 ggTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTtct        >  1:4532956/1‑93 (MQ=255)
                                                                                  gTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTtctc       <  1:2415986/93‑1 (MQ=255)
                                                                                  gTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTtctc       <  2:2753274/93‑1 (MQ=255)
                                                                                   ttAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCg      >  1:3841985/1‑93 (MQ=255)
                                                                                     aaaCCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCt    <  1:4295262/93‑1 (MQ=255)
                                                                                     aaaCCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCt    <  2:3887972/93‑1 (MQ=255)
                                                                                      aaCCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCtt   >  1:1628832/1‑93 (MQ=255)
                                                                                      aaCCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCtt   >  1:2613962/1‑93 (MQ=255)
                                                                                      aaCCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCtt   >  1:4386299/1‑93 (MQ=255)
                                                                                      aaCCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCtt   <  2:2024819/93‑1 (MQ=255)
                                                                                       aCCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTTTAGTAATGGGTCACTCTTTTTTCTCGCttt  <  1:1806029/93‑1 (MQ=255)
                                                                                       aCCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCttt  <  1:2486085/93‑1 (MQ=255)
                                                                                          |                                                                                         
CTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTT  >  NC_000913/4621554‑4621733

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