Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 612,751 T→C A238A (GCT→GCC fes → enterochelin esterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913612,7510TC99.7% 1083.6 / ‑4.8 313A238A (GCT→GCCfesenterochelin esterase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (148/164);  total (149/164)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.20e-01

Only 100 of 320 total aligned reads displayed.

TATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTG  >  NC_000913/612661‑612840
                                                                                          |                                                                                         
tATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCAt                                                                                         <  1:3681141/93‑1 (MQ=255)
tATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCAt                                                                                         >  2:6544587/1‑93 (MQ=255)
 aTCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCAtt                                                                                        <  2:30480/93‑1 (MQ=255)
 aTCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCAtt                                                                                        <  1:185165/93‑1 (MQ=255)
 aTCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTAATACCCCCGGTGAAAGCCAtt                                                                                        >  1:831996/1‑93 (MQ=255)
  tCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTg                                                                                       >  2:8380160/1‑93 (MQ=255)
     acacCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGccc                                                                                    >  2:7605038/1‑93 (MQ=255)
     acacCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGccc                                                                                    <  1:8782796/93‑1 (MQ=255)
      caccacGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGcccc                                                                                   <  2:7093679/93‑1 (MQ=255)
       accacGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCt                                                                                  <  2:1867714/93‑1 (MQ=255)
       accacGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGACTCGCAGTACAGCAAGAGTTATTCCCCCTGGTGAAAGCCATTGCCCCt                                                                                  <  2:7372239/93‑1 (MQ=255)
         cacGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCttt                                                                                >  1:556818/1‑93 (MQ=255)
         cacGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCttt                                                                                <  1:1457465/93‑1 (MQ=255)
          acGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCtttt                                                                               >  1:5417259/1‑93 (MQ=255)
             cACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGc                                                                            <  1:8366701/93‑1 (MQ=255)
              aCCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCg                                                                           <  1:4985116/93‑1 (MQ=255)
               ccGCGCCCACGAACTGCCGTGTATTGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGa                                                                          <  2:4757020/93‑1 (MQ=255)
               ccGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGa                                                                          <  2:1739611/93‑1 (MQ=255)
                cgcgCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGAt                                                                         <  2:325790/93‑1 (MQ=255)
                 gcgcCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATc                                                                        <  1:2030126/93‑1 (MQ=255)
                 gcgcCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATc                                                                        >  1:1357724/1‑93 (MQ=255)
                 gcgcCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATc                                                                        <  1:2800371/93‑1 (MQ=255)
                  cgcCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCg                                                                       <  1:6933384/93‑1 (MQ=255)
                  cgcCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCg                                                                       <  1:8101682/93‑1 (MQ=255)
                  cgcCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCg                                                                       <  1:5436491/93‑1 (MQ=255)
                   gcCCCCGCACTGCCGGGTAATTCGGATTTCTGGCTCGCAGTACAGCAAGTGTTATCCCCCCGGGTGAAAGCCATTGCCCCTTTTAGCGATCGt                                                                      <  1:5927072/93‑1 (MQ=255)
                     ccACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGc                                                                    <  1:6527535/93‑1 (MQ=255)
                      cACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGcc                                                                   >  1:20737/1‑93 (MQ=255)
                       aCGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCg                                                                  >  2:1364892/1‑93 (MQ=255)
                         gAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGAt                                                                >  2:2164282/1‑93 (MQ=255)
                          aaCTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATc                                                               >  1:2977781/1‑93 (MQ=255)
                          aaCTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATc                                                               <  1:1588184/93‑1 (MQ=255)
                           aCTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCg                                                              <  2:1394658/93‑1 (MQ=255)
                            cTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGc                                                             >  1:4800555/1‑93 (MQ=255)
                            cTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGc                                                             >  1:35128/1‑93 (MQ=255)
                            cTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGc                                                             >  1:2230208/1‑93 (MQ=255)
                            cTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGc                                                             >  1:7907876/1‑93 (MQ=255)
                             tGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCa                                                            <  1:2513306/93‑1 (MQ=255)
                              gCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCAc                                                           <  1:5729219/93‑1 (MQ=255)
                                cGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCg                                                         <  1:6628187/93‑1 (MQ=255)
                                 gtgtAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGt                                                        >  1:6153531/1‑93 (MQ=255)
                                 gtgtAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGt                                                        >  2:5577330/1‑93 (MQ=255)
                                  tgtAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTg                                                       >  2:5479342/1‑93 (MQ=255)
                                  tgtAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTg                                                       >  1:1113721/1‑93 (MQ=255)
                                   gtAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTgg                                                      <  2:2567994/93‑1 (MQ=255)
                                    tAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGt                                                     <  1:855122/93‑1 (MQ=255)
                                      aTGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTg                                                   <  2:4081612/93‑1 (MQ=255)
                                           gATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCggg                                              <  2:7473010/93‑1 (MQ=255)
                                            aTTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGc                                             <  2:6278191/93‑1 (MQ=255)
                                                cTGGCTCGGAGTACAGCAAGAGTTATTCCCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCagag                                         <  2:2059492/93‑1 (MQ=255)
                                                cTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCagag                                         >  1:3832697/1‑93 (MQ=255)
                                                cTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCagag                                         >  2:3417080/1‑93 (MQ=255)
                                                 tGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGt                                        >  1:4106421/1‑93 (MQ=255)
                                                  ggCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGtt                                       >  2:4872607/1‑93 (MQ=255)
                                                    cTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGtttt                                     >  1:2768537/1‑93 (MQ=255)
                                                     tCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTtg                                    >  2:4398781/1‑93 (MQ=255)
                                                     tCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTtg                                    >  1:2135267/1‑93 (MQ=255)
                                                     ccGCAGTACAGCAAGAGTTATTTCCCCTGGTGAAATCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTtg                                    <  1:2567588/92‑1 (MQ=255)
                                                      cGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTtgg                                   <  2:8300956/93‑1 (MQ=255)
                                                       gCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTtggt                                  >  2:5387979/1‑93 (MQ=255)
                                                       gCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTtggt                                  <  2:8174668/93‑1 (MQ=255)
                                                       gCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGGTTtggt                                  >  1:4582732/1‑93 (MQ=255)
                                                         aGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTtggtgg                                <  2:5624772/93‑1 (MQ=255)
                                                          gTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTggg                               <  2:6856740/93‑1 (MQ=255)
                                                          gTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTggg                               >  2:7022630/1‑93 (MQ=255)
                                                             cAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCtt                            <  2:1560120/93‑1 (MQ=255)
                                                              aGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCttt                           <  2:6918547/93‑1 (MQ=255)
                                                               gCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTc                          <  2:6465668/93‑1 (MQ=255)
                                                               gCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTc                          <  1:7269285/93‑1 (MQ=255)
                                                                cAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTcc                         <  1:7392516/93‑1 (MQ=255)
                                                                cAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTcc                         <  2:5098358/93‑1 (MQ=255)
                                                                 aaGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCcg                        >  1:3382442/1‑93 (MQ=255)
                                                                    agTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGGACCGTGGCTGCCGGGCAGAGTTTTGGTGGGTCTTCCccgc                     >  1:1042814/1‑93 (MQ=255)
                                                                      ttattaCCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTg                   >  2:1308299/1‑93 (MQ=255)
                                                                       tatttCCCCTGGTGAAACCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGt                  <  2:3861640/93‑1 (MQ=255)
                                                                       tattaCCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGt                  >  2:794452/1‑93 (MQ=255)
                                                                        attaCCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTa                 <  2:3506848/93‑1 (MQ=255)
                                                                         ttaCCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTAt                <  1:6711995/93‑1 (MQ=255)
                                                                          taCCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATg               <  1:6301220/93‑1 (MQ=255)
                                                                          taCCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGAGCTTTCCGCGCTGTATg               <  1:4271145/93‑1 (MQ=255)
                                                                           tccccTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGc              <  1:6330382/92‑1 (MQ=255)
                                                                           aCCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGc              <  2:5752654/93‑1 (MQ=255)
                                                                            ccccTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGcc             >  2:4688401/1‑93 (MQ=255)
                                                                             cccTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCg            <  1:6374480/93‑1 (MQ=255)
                                                                             cccTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCg            <  1:7136824/93‑1 (MQ=255)
                                                                              ccTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCgg           >  2:1184294/1‑93 (MQ=255)
                                                                              ccTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCgg           >  1:7456747/1‑93 (MQ=255)
                                                                                tGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGAc         <  1:7577134/93‑1 (MQ=255)
                                                                                tGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGAc         <  1:5625836/93‑1 (MQ=255)
                                                                                tGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGAc         <  1:7655661/93‑1 (MQ=255)
                                                                                  gTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTg       >  2:5928730/1‑93 (MQ=255)
                                                                                  gTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTg       >  1:3113235/1‑93 (MQ=255)
                                                                                   tGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTGGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGc      <  2:7715100/93‑1 (MQ=255)
                                                                                    gAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCa     <  2:105243/93‑1 (MQ=255)
                                                                                    gAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCa     >  2:8463322/1‑93 (MQ=255)
                                                                                     aaaGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCAc    <  1:2657570/93‑1 (MQ=255)
                                                                                     aaaGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCAc    >  2:7756297/1‑93 (MQ=255)
                                                                                      aaGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACt   >  2:2214193/1‑93 (MQ=255)
                                                                                       aGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTg  >  1:8617575/1‑93 (MQ=255)
                                                                                       aGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTg  >  1:7213972/1‑93 (MQ=255)
                                                                                          |                                                                                         
TATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTG  >  NC_000913/612661‑612840

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: