Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,158,718 T→C Q686R (CAG→CGG)  fhuE ← outer membrane receptor for ferric iron uptake

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,158,7180TC83.4% 501.5 / 100.9 229Q686R (CAG→CGG) fhuEouter membrane receptor for ferric iron uptake
Reads supporting (aligned to +/- strand):  ref base T (24/14);  new base C (116/75);  total (140/89)
Fisher's exact test for biased strand distribution p-value = 8.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-02

Only 100 of 229 total aligned reads displayed.

CCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGA  >  NC_000913/1158628‑1158807
                                                                                          |                                                                                         
  gTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTa                                                                                       >  2:1226566/1‑93 (MQ=255)
  gTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTa                                                                                       <  2:5354890/93‑1 (MQ=255)
   tAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAg                                                                                      >  1:2252176/1‑93 (MQ=255)
     gACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCg                                                                                    >  1:4741084/1‑93 (MQ=255)
      aCGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCgg                                                                                   >  2:3017586/1‑93 (MQ=255)
      aCGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCgg                                                                                   <  2:3331403/93‑1 (MQ=255)
      aCGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCgg                                                                                   <  2:4157840/93‑1 (MQ=255)
      aCGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCgg                                                                                   <  2:5016690/93‑1 (MQ=255)
         aTAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGTGGGTg                                                                                >  2:3419187/1‑93 (MQ=255)
         aTAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTg                                                                                >  2:1974827/1‑93 (MQ=255)
         aTAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTg                                                                                >  2:4760235/1‑93 (MQ=255)
         aTAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTg                                                                                >  2:4769278/1‑93 (MQ=255)
         aTAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTg                                                                                >  2:5140202/1‑93 (MQ=255)
          tAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGa                                                                               >  2:886741/1‑93 (MQ=255)
              aCCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAagag                                                                           >  1:6562357/1‑93 (MQ=255)
              aCCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAagag                                                                           >  1:5438147/1‑93 (MQ=255)
              aCCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAagag                                                                           >  2:4436681/1‑93 (MQ=255)
               ccTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAgaga                                                                          >  1:7472714/1‑93 (MQ=255)
               ccTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAgaga                                                                          >  2:2008058/1‑93 (MQ=255)
                 ttCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATc                                                                        >  1:3892287/1‑93 (MQ=255)
                 ttCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATc                                                                        <  2:5755766/93‑1 (MQ=255)
                  tCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATcc                                                                       >  1:6796661/1‑93 (MQ=255)
                  tCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATcc                                                                       <  1:5512212/93‑1 (MQ=255)
                   ccACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATcca                                                                      <  2:5258937/93‑1 (MQ=255)
                    cACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATccac                                                                     >  2:1014180/1‑93 (MQ=255)
                    cACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATccac                                                                     >  1:686614/1‑93 (MQ=255)
                     aCGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATccacc                                                                    >  2:6290504/1‑93 (MQ=255)
                      cGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATccacca                                                                   <  1:3042249/93‑1 (MQ=255)
                       gTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAg                                                                  <  2:6024049/93‑1 (MQ=255)
                       gTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAg                                                                  >  1:2915043/1‑93 (MQ=255)
                          ggTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAgcgc                                                               <  2:4084784/93‑1 (MQ=255)
                           gTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGcgcg                                                              <  1:132621/93‑1 (MQ=255)
                           gTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGcgcg                                                              >  1:3309149/1‑93 (MQ=255)
                              tCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAg                                                           >  1:4190904/1‑93 (MQ=255)
                               cGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGc                                                          <  2:4944580/93‑1 (MQ=255)
                               cGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGc                                                          <  1:1747565/93‑1 (MQ=255)
                                 tAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTa                                                        >  1:4904130/1‑93 (MQ=255)
                                 tAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTa                                                        >  1:5611714/1‑93 (MQ=255)
                                 tAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTa                                                        >  2:7166706/1‑93 (MQ=255)
                                 tAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTa                                                        >  1:3184126/1‑93 (MQ=255)
                                   ggTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTAcc                                                      >  1:170360/1‑93 (MQ=255)
                                   ggTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTAcc                                                      >  2:5462825/1‑93 (MQ=255)
                                    gTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAATGAGAAATTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCt                                                     >  2:1533616/1‑93 (MQ=255)
                                    gTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCt                                                     >  2:356914/1‑93 (MQ=255)
                                    gTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCt                                                     >  1:433686/1‑93 (MQ=255)
                                     ttttGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCtt                                                    <  1:4539239/93‑1 (MQ=255)
                                     ttttGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCtt                                                    <  2:5726400/93‑1 (MQ=255)
                                       ttGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGc                                                  >  1:414879/1‑93 (MQ=255)
                                       ttGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGGTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGc                                                  >  1:6215382/1‑93 (MQ=255)
                                           cGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCgg                                              <  1:3797594/93‑1 (MQ=255)
                                            gAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGgc                                             >  2:6779970/1‑93 (MQ=255)
                                            gAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGgc                                             <  2:4007503/93‑1 (MQ=255)
                                            gAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGgc                                             >  2:1599267/1‑93 (MQ=255)
                                               cAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGcgcg                                          >  2:2707201/1‑93 (MQ=255)
                                               cAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGcgcg                                          >  2:5201419/1‑93 (MQ=255)
                                                aGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCgg                                         >  1:4105904/1‑93 (MQ=255)
                                                aGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCgg                                         >  1:5341331/1‑93 (MQ=255)
                                                aGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCgg                                         >  2:3867260/1‑93 (MQ=255)
                                                  gTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGaa                                       <  1:2594428/93‑1 (MQ=255)
                                                  gTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGaa                                       <  1:2193940/93‑1 (MQ=255)
                                                   ttATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAg                                      >  1:2262819/1‑93 (MQ=255)
                                                     aTTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGt                                    >  2:5596547/1‑93 (MQ=255)
                                                      ttGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTg                                   >  2:6638803/1‑93 (MQ=255)
                                                       tGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGc                                  <  2:6163313/93‑1 (MQ=255)
                                                         aCGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCa                                <  2:6531646/93‑1 (MQ=255)
                                                         aCGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCa                                >  2:7469938/1‑93 (MQ=255)
                                                         aCGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCa                                <  2:2878874/93‑1 (MQ=255)
                                                         aCGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCa                                <  1:1497547/93‑1 (MQ=255)
                                                           gTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATa                              >  1:1171923/1‑93 (MQ=255)
                                                           gTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATa                              >  1:1932835/1‑93 (MQ=255)
                                                              ccccTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACgg                           >  2:3620477/1‑93 (MQ=255)
                                                                 cTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGt                        >  2:4865904/1‑93 (MQ=255)
                                                                  tGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTc                       >  2:5930484/1‑93 (MQ=255)
                                                                  tGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTc                       >  2:173005/1‑93 (MQ=255)
                                                                  tGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTc                       >  1:6497232/1‑93 (MQ=255)
                                                                   gTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCa                      <  2:4335739/93‑1 (MQ=255)
                                                                   gTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCa                      <  1:1027654/93‑1 (MQ=255)
                                                                   gTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCa                      >  1:2547424/1‑93 (MQ=255)
                                                                    tAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCCTCAc                     <  2:7370391/93‑1 (MQ=255)
                                                                     aaCGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACg                    >  2:668173/1‑93 (MQ=255)
                                                                     aaCGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACg                    >  1:2063425/1‑93 (MQ=255)
                                                                      aCGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACgg                   <  1:4425405/93‑1 (MQ=255)
                                                                      aCGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACgg                   >  2:223052/1‑93 (MQ=255)
                                                                      aCGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACgg                   <  2:1331398/93‑1 (MQ=255)
                                                                      aCGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACgg                   >  1:2826846/1‑93 (MQ=255)
                                                                       cGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGgt                  <  2:4549802/93‑1 (MQ=255)
                                                                         agaAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGgtgt                >  1:513903/1‑93 (MQ=255)
                                                                          gaAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTc               <  1:625964/93‑1 (MQ=255)
                                                                          gaAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTc               >  2:3398893/1‑93 (MQ=255)
                                                                           aaGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCg              >  1:4038345/1‑93 (MQ=255)
                                                                              tttttCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGta           <  1:5897526/93‑1 (MQ=255)
                                                                              tttttCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGta           <  2:6098723/93‑1 (MQ=255)
                                                                                tttCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGtata         <  2:2482226/93‑1 (MQ=255)
                                                                                tttCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGtata         >  1:3574200/1‑93 (MQ=255)
                                                                                 ttCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATac        <  2:3438685/93‑1 (MQ=255)
                                                                                  tCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATaca       <  2:7029986/93‑1 (MQ=255)
                                                                                   cGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATacac      >  1:716840/1‑93 (MQ=255)
                                                                                     tCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACAcgc    <  2:3841720/93‑1 (MQ=255)
                                                                                       aCCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGa  >  1:5909233/1‑93 (MQ=255)
                                                                                       aCCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGa  >  1:4958791/1‑93 (MQ=255)
                                                                                          |                                                                                         
CCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGA  >  NC_000913/1158628‑1158807

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