Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,198,118 A→G D87G (GAC→GGC)  hupA → DNA‑binding protein HU‑alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,198,1180AG99.5% 684.4 / ‑6.4 212D87G (GAC→GGC) hupADNA‑binding protein HU‑alpha
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (111/100);  total (112/100)
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 223 total aligned reads displayed.

GCGTACTGGCCGCAACCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGC  >  NC_000913/4198029‑4198207
                                                                                         |                                                                                         
gCGTACTGGCCGCAACCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGgcgc                                                                                        <  1:2692046/93‑1 (MQ=255)
gCGTACTGGCCGCAACCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGgcgc                                                                                        <  2:3923534/93‑1 (MQ=255)
 cGTACTGGCCGCAACCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCa                                                                                       <  1:3902785/93‑1 (MQ=255)
  gTACTGGCCGCAACCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAg                                                                                      <  2:630421/93‑1 (MQ=255)
   tACTGGCCGCAACCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGt                                                                                     >  2:88418/1‑93 (MQ=255)
     cTGGCCGCAACCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTta                                                                                   >  2:4703341/1‑93 (MQ=255)
      tGGCCGCAACCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTtaa                                                                                  <  2:1369229/93‑1 (MQ=255)
       ggCCGCAACCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTtaag                                                                                 >  1:761599/1‑93 (MQ=255)
        gCCGCAACCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTtaagt                                                                                <  1:2217430/93‑1 (MQ=255)
        gCCGCAACCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTtaagt                                                                                <  1:345301/93‑1 (MQ=255)
        gCCGCAACCCGCAGACCGGTAAAGAAATCAAAATTACCGCAGCTAACGTACCGGCATTTGTTTCTGCCAAGGCACTGAAAGGCGCAGTtaaga                                                                                >  1:1635501/1‑92 (MQ=38)
           gCAACCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTtaagtaag                                                                             >  1:2781425/1‑93 (MQ=255)
           gCAACCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTtaagtaag                                                                             >  1:3081494/1‑93 (MQ=255)
             aaCCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGAt                                                                           >  1:1011858/1‑93 (MQ=255)
              aCCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGAtt                                                                          >  2:3883475/1‑93 (MQ=255)
                ccGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCTTTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGc                                                                        <  2:3236345/93‑1 (MQ=255)
                 cGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCg                                                                       <  1:3491459/93‑1 (MQ=255)
                  gCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGt                                                                      <  2:1489933/93‑1 (MQ=255)
                  gCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGt                                                                      >  2:4488775/1‑93 (MQ=255)
                   cAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTg                                                                     <  2:2387050/93‑1 (MQ=255)
                   cAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTg                                                                     >  2:3553772/1‑93 (MQ=255)
                    aGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTgg                                                                    >  2:4133987/1‑93 (MQ=255)
                       ccGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAg                                                                 <  1:4578891/93‑1 (MQ=255)
                       ccGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAg                                                                 <  1:2421273/93‑1 (MQ=255)
                        cGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGt                                                                <  2:1151751/93‑1 (MQ=255)
                         ggTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTg                                                               <  2:509924/93‑1 (MQ=255)
                         ggTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTg                                                               <  2:456999/93‑1 (MQ=255)
                          gTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGAAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGa                                                              <  2:4334857/93‑1 (MQ=255)
                            aaaGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAAc                                                            >  1:608311/1‑93 (MQ=255)
                             aaGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACa                                                           >  1:531008/1‑93 (MQ=255)
                              aGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAg                                                          <  1:152632/93‑1 (MQ=255)
                               gAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGt                                                         <  2:2825079/93‑1 (MQ=255)
                                aaaTCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGtt                                                        <  1:1864811/93‑1 (MQ=255)
                                 aaTCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGttt                                                       <  2:2373243/93‑1 (MQ=255)
                                  aTCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGtttt                                                      <  1:2591252/93‑1 (MQ=255)
                                  aTCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGtttt                                                      <  1:4369958/93‑1 (MQ=255)
                                   tCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTa                                                     >  1:198624/1‑93 (MQ=255)
                                   tCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTa                                                     >  2:3824841/1‑93 (MQ=255)
                                     aaaaTTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAAc                                                   <  2:4405809/93‑1 (MQ=255)
                                       aaTTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGa                                                 <  2:934531/93‑1 (MQ=255)
                                       aaTTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGa                                                 >  2:1181432/1‑93 (MQ=255)
                                          tGCCGCAGCTAGCGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTCAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAgg                                              <  2:4212915/93‑1 (MQ=255)
                                          tGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAgg                                              <  2:3077051/93‑1 (MQ=255)
                                            ccGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAgggg                                            <  2:604142/93‑1 (MQ=255)
                                             cGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGggt                                           <  1:2357354/93‑1 (MQ=255)
                                              gCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGggtg                                          <  2:362287/93‑1 (MQ=255)
                                              gCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGggtg                                          >  1:3315838/1‑93 (MQ=255)
                                              gCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGggtg                                          >  1:705881/1‑93 (MQ=255)
                                                aGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGTggtggt                                        >  2:4354095/1‑93 (MQ=255)
                                                aGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGggtggt                                        >  2:1821171/1‑93 (MQ=255)
                                                 gCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGtt                                       <  2:782411/93‑1 (MQ=255)
                                                 gCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGtt                                       <  2:3482273/93‑1 (MQ=255)
                                                  cTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGttt                                      >  1:3152658/1‑93 (MQ=255)
                                                   tAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGAGAACAGTTTTAACGAAGGGGTGGGTTc                                     >  2:3695510/1‑93 (MQ=255)
                                                    aaCGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCa                                    >  1:511962/1‑93 (MQ=255)
                                                       gTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCAccc                                 <  1:3385122/93‑1 (MQ=255)
                                                       gTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGCTTTAACGAAGGGGCGGTTTCAccc                                 >  1:2295924/1‑93 (MQ=255)
                                                       gTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGACAGTGAACAGGTTTAACAAAGTGAGGGGTTCActc                                 >  2:3902612/1‑91 (MQ=255)
                                                        tACCGGCATTTGTTTCTGGCAAGGCTCTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCAcccc                                >  1:805845/1‑93 (MQ=255)
                                                        tACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCAcccc                                >  1:3157498/1‑93 (MQ=255)
                                                        tACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCAcccc                                >  1:4250164/1‑93 (MQ=255)
                                                          ccGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCtt                              <  1:739641/93‑1 (MQ=255)
                                                           cGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCttt                             <  2:2914477/93‑1 (MQ=255)
                                                            ggCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCtttt                            >  2:1992155/1‑93 (MQ=255)
                                                            ggCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCtttt                            <  1:1356635/93‑1 (MQ=255)
                                                            ggCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCtttt                            <  2:1677091/93‑1 (MQ=255)
                                                               aTTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTc                         >  2:2114977/1‑93 (MQ=255)
                                                               aTTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTc                         <  2:3106211/93‑1 (MQ=255)
                                                               aTTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTc                         >  2:1936526/1‑93 (MQ=255)
                                                                tttGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCt                        >  1:2251782/1‑93 (MQ=255)
                                                                tttGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCt                        <  1:2430626/93‑1 (MQ=255)
                                                                tttGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCt                        >  1:2729041/1‑93 (MQ=255)
                                                                 ttGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCtt                       >  2:1425441/1‑93 (MQ=255)
                                                                 ttGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCtt                       >  2:2566801/1‑93 (MQ=255)
                                                                    tttCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCt                    <  2:1472291/93‑1 (MQ=255)
                                                                    tttCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCt                    <  2:3275391/93‑1 (MQ=255)
                                                                    tttCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCt                    <  1:4088326/93‑1 (MQ=255)
                                                                    tttCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCt                    <  1:3859663/93‑1 (MQ=255)
                                                                    tttCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGATGGTTTCACCCCTTTTGTCTTTCt                    >  1:4024048/1‑93 (MQ=255)
                                                                      tCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTgg                  <  2:4185633/93‑1 (MQ=255)
                                                                       cTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGc                 <  2:1244120/93‑1 (MQ=255)
                                                                        tGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGGTTCACCCCTTTTGTCTTTCTGGCg                >  1:1771065/1‑93 (MQ=255)
                                                                            aaGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGa            <  2:1327170/93‑1 (MQ=255)
                                                                             aGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCAGCCCTTTTGTCTTTCTGGCGTCGAt           >  2:735682/1‑93 (MQ=255)
                                                                                cACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCAt        >  1:132880/1‑93 (MQ=255)
                                                                                cACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCAt        <  2:2330960/93‑1 (MQ=255)
                                                                                cACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCAt        >  2:1498635/1‑93 (MQ=255)
                                                                                cACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCAt        >  1:3156072/1‑93 (MQ=255)
                                                                                 aCTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCAtt       <  1:3698542/93‑1 (MQ=255)
                                                                                  cTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTg      >  1:4432995/1‑93 (MQ=255)
                                                                                  cTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTg      >  2:3158882/1‑93 (MQ=255)
                                                                                   tGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGa     >  1:2106990/1‑93 (MQ=255)
                                                                                    gAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGAt    <  1:3750916/93‑1 (MQ=255)
                                                                                    gAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTACTTCTGGCGTCGATCATTGAt    >  1:1747958/1‑93 (MQ=255)
                                                                                    gAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGAGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGAt    >  1:3467602/1‑93 (MQ=255)
                                                                                      aaGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATgc  <  1:3746498/93‑1 (MQ=255)
                                                                                      aaGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATgc  <  1:1492198/93‑1 (MQ=255)
                                                                                      aaGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATgc  <  2:4269936/93‑1 (MQ=255)
                                                                                      aaGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATgc  <  2:2804543/93‑1 (MQ=255)
                                                                                      aaGGCGCAGTTAAGTAAGATTGCGGGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATgc  >  2:2361564/1‑93 (MQ=255)
                                                                                         |                                                                                         
GCGTACTGGCCGCAACCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGC  >  NC_000913/4198029‑4198207

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: ‑