Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,374,263 T→C Q50R (CAA→CGA)  sspB ← stringent starvation protein B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,374,2630TC99.7% 990.6 / ‑5.4 313Q50R (CAA→CGA) sspBstringent starvation protein B
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (171/140);  total (171/142)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

Only 100 of 329 total aligned reads displayed.

GGAATGCCACCAAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTTGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATCCAGCAAC  >  NC_000913/3374173‑3374352
                                                                                          |                                                                                         
ggAATGCCACCAAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGc                                                                                         >  1:3846632/1‑93 (MQ=255)
ggAATGCCACCAAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGc                                                                                         >  2:4893941/1‑93 (MQ=255)
       caccaAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCAc                                                                                  >  1:6546561/1‑93 (MQ=255)
       caccaAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCAc                                                                                  >  2:1044653/1‑93 (MQ=255)
       caccaAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCAc                                                                                  >  2:1107015/1‑93 (MQ=255)
        cccaAAGCGCGCGTTAAAGCGCATCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCAcg                                                                                 <  2:4295485/92‑1 (MQ=255)
        accaAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCAcg                                                                                 >  1:3873631/1‑93 (MQ=255)
        accaAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCAcg                                                                                 <  1:2393413/93‑1 (MQ=255)
         ccaAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCAcgc                                                                                <  2:6926370/93‑1 (MQ=255)
           aaaGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACgcgc                                                                              >  2:7693769/1‑93 (MQ=255)
           aaaGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACgcgc                                                                              >  2:465298/1‑93 (MQ=255)
           aaaGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACgcgc                                                                              <  2:5519467/93‑1 (MQ=255)
              gcgcgcGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCata                                                                           <  2:1997120/93‑1 (MQ=255)
              gcgcgcGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCata                                                                           <  1:4320598/93‑1 (MQ=255)
               cgcgcgTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCatat                                                                          <  2:1375087/93‑1 (MQ=255)
                 cgcgTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTc                                                                        <  2:6224857/93‑1 (MQ=255)
                 cgcgTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTc                                                                        <  1:109070/93‑1 (MQ=255)
                  gcgTTAGAGCGCACCCCACCTTTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCACTGTTGAGTACGATTCCCCCGTCACGCGCATATTcc                                                                       <  2:6606317/93‑1 (MQ=255)
                  gcgTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTcc                                                                       <  2:2679116/93‑1 (MQ=255)
                    gTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCAt                                                                     <  1:2793433/93‑1 (MQ=255)
                    gTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCAt                                                                     >  2:1552385/1‑93 (MQ=255)
                     ttAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATa                                                                    <  1:7838008/93‑1 (MQ=255)
                          gcgcACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAAc                                                               <  2:5617598/93‑1 (MQ=255)
                          gcccACCTTATCATTCACCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTAGGATTCGCCCGTCACGCGCATATTCCATAGGAAc                                                               <  2:2311826/90‑1 (MQ=255)
                             cACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTg                                                            >  1:6256419/1‑93 (MQ=255)
                               ccTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGACCCTGCa                                                          <  2:4154905/93‑1 (MQ=255)
                                cTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCAc                                                         >  1:5308505/1‑93 (MQ=255)
                                 tcatcaTTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACg                                                        >  1:7897276/1‑93 (MQ=255)
                                   atcatTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGcc                                                      >  1:6178944/1‑93 (MQ=255)
                                   atcatTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGcc                                                      >  2:916058/1‑93 (MQ=255)
                                    tcatTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCa                                                     >  2:7234455/1‑93 (MQ=255)
                                    tcatTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCa                                                     >  2:7216538/1‑93 (MQ=255)
                                    tcatTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCa                                                     >  1:5790603/1‑93 (MQ=255)
                                    tcatTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCa                                                     >  1:4957034/1‑93 (MQ=255)
                                     catTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAg                                                    <  1:6423858/93‑1 (MQ=255)
                                      atTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAgg                                                   >  1:4095892/1‑93 (MQ=255)
                                         cGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAg                                                <  1:1784362/93‑1 (MQ=255)
                                         cGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAg                                                >  2:5203042/1‑93 (MQ=255)
                                         cGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAg                                                <  2:4060884/93‑1 (MQ=255)
                                          gCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGc                                               >  1:2216097/1‑93 (MQ=255)
                                          gCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGc                                               <  1:5993433/93‑1 (MQ=255)
                                           ccAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCATGCCAGGGAGCg                                              >  2:5780960/1‑93 (MQ=255)
                                             atgtCCATATTGCCAACGGCACGCGGAGAATTGTTGAGTACGATTCGACCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTc                                            <  2:4478369/90‑1 (MQ=255)
                                             aGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTc                                            >  1:757234/1‑93 (MQ=255)
                                             aGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTc                                            <  1:1940535/93‑1 (MQ=255)
                                              gTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTca                                           >  1:4843714/1‑93 (MQ=255)
                                              gTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTca                                           >  1:2097314/1‑93 (MQ=255)
                                               ttCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTcac                                          >  1:4059881/1‑93 (MQ=255)
                                                   aGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATcc                                      <  1:6481838/93‑1 (MQ=255)
                                                   aGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATcc                                      <  2:7867247/93‑1 (MQ=255)
                                                   aGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATcc                                      >  1:7615277/1‑93 (MQ=255)
                                                    gATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATccc                                     >  2:6051063/1‑92 (MQ=255)
                                                    gATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATcca                                     >  2:1219654/1‑93 (MQ=255)
                                                    gATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATcca                                     >  2:3057272/1‑93 (MQ=255)
                                                    gATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATcca                                     >  2:3293637/1‑93 (MQ=255)
                                                    gATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATcca                                     >  2:1731203/1‑93 (MQ=255)
                                                    gATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATcca                                     >  1:7143821/1‑93 (MQ=255)
                                                     aTTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATccac                                    >  2:3457280/1‑93 (MQ=255)
                                                     aTTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATccac                                    <  2:2244014/93‑1 (MQ=255)
                                                        gCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATccaccac                                 <  2:1287829/93‑1 (MQ=255)
                                                         ccGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATccaccacc                                <  1:3629856/93‑1 (MQ=255)
                                                         ccGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATccaccacc                                <  1:3430716/93‑1 (MQ=255)
                                                          cGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATccaccacca                               <  2:6297407/93‑1 (MQ=255)
                                                            aCAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCAGATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAgg                             <  1:4210749/93‑1 (MQ=255)
                                                             cAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGt                            <  2:6589928/93‑1 (MQ=255)
                                                             cAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGt                            <  1:1517199/93‑1 (MQ=255)
                                                             cAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGt                            <  1:3123345/93‑1 (MQ=255)
                                                             cAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGt                            <  1:4578807/93‑1 (MQ=255)
                                                              aGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTg                           >  1:1721515/1‑93 (MQ=255)
                                                               gCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTgc                          >  2:4782508/1‑93 (MQ=255)
                                                               gCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTgc                          <  1:4641689/93‑1 (MQ=255)
                                                                  cgcgCCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTgcggc                       <  1:3949269/93‑1 (MQ=255)
                                                                      gcgcAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGa                   <  1:1162504/93‑1 (MQ=255)
                                                                      gcgcAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGa                   <  2:3830735/93‑1 (MQ=255)
                                                                       cgcAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAg                  <  2:5410487/93‑1 (MQ=255)
                                                                        gcAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGc                 <  2:7455324/93‑1 (MQ=255)
                                                                        gcAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGc                 >  2:933489/1‑93 (MQ=255)
                                                                         cAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCt                >  1:1233899/1‑93 (MQ=255)
                                                                          aaTGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTg               >  2:5745278/1‑93 (MQ=255)
                                                                          aaTGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTg               <  1:7371266/93‑1 (MQ=255)
                                                                           aTGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTgg              >  1:325604/1‑93 (MQ=255)
                                                                           aTGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTgg              >  2:7329315/1‑93 (MQ=255)
                                                                            tGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGt             >  2:6479290/1‑93 (MQ=255)
                                                                            tGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGt             >  1:1027980/1‑93 (MQ=255)
                                                                            tGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGt             <  1:3787443/93‑1 (MQ=255)
                                                                            tGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGt             >  2:5478097/1‑93 (MQ=255)
                                                                            tGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGt             >  1:2912271/1‑93 (MQ=255)
                                                                            tGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGt             <  1:6888732/93‑1 (MQ=255)
                                                                            tGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGt             <  1:5181892/93‑1 (MQ=255)
                                                                              ttGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTa           <  1:742947/93‑1 (MQ=255)
                                                                                gAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATc         <  2:6188625/93‑1 (MQ=255)
                                                                                 aGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATcc        <  1:951619/93‑1 (MQ=255)
                                                                                  gTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATCCa       <  2:3681228/93‑1 (MQ=255)
                                                                                  gTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATCCa       >  2:1929175/1‑93 (MQ=255)
                                                                                  gTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATCCa       >  2:187226/1‑93 (MQ=255)
                                                                                   tACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACACCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATCCAg      >  2:3870711/1‑93 (MQ=255)
                                                                                     cGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATCCAGCa    >  1:7282157/1‑93 (MQ=255)
                                                                                     cGATTCGCCCGTCACGCGCATAGTCCATAGGAACCTGCACGCCAGGGAGCGTCACATTCACCACCAGGTGCGGCGTGAGCTGGGTATCCATCa    >  2:7037393/1‑93 (MQ=255)
                                                                                      gATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATCCAGCaa   >  1:56794/1‑93 (MQ=255)
                                                                                       aTTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATCCAGCAAc  >  1:3527457/1‑93 (MQ=255)
                                                                                          |                                                                                         
GGAATGCCACCAAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTTGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATCCAGCAAC  >  NC_000913/3374173‑3374352

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: