Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,539,774 G→A A578T (GCC→ACC)  feoB → ferrous iron transport protein B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,539,7740GA99.6% 841.7 / NA 234A578T (GCC→ACC) feoBferrous iron transport protein B
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (107/125);  total (107/127)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 238 total aligned reads displayed.

GCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGAATTT  >  NC_000913/3539684‑3539863
                                                                                          |                                                                                         
 cgtcggtcAGCCGGGTGATCACCCCGGTCTGCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCa                                                                                        <  2:2748105/93‑1 (MQ=255)
  gtccgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCAt                                                                                       <  2:5564336/93‑1 (MQ=255)
   tccgtcAGCCGGGTGACCCCCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATg                                                                                      <  1:1657142/93‑1 (MQ=255)
   cccgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATg                                                                                      <  2:206909/92‑1 (MQ=255)
       tcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGa                                                                                  >  2:168614/1‑93 (MQ=255)
       tcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGa                                                                                  >  1:1585414/1‑93 (MQ=255)
         cgCCGGGTGATCCCCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGaaa                                                                                <  1:2906189/92‑1 (MQ=255)
         aGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGaaa                                                                                <  1:3365347/93‑1 (MQ=255)
            cGGGTGATCCCCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAaagaa                                                                             <  1:3923827/93‑1 (MQ=255)
            cGGGGGATCACCCCGGTCTTCAAGACAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAaagaa                                                                             >  1:1430992/1‑93 (MQ=255)
             gggTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAaagaag                                                                            >  2:6126964/1‑93 (MQ=255)
              ggTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGt                                                                           >  1:3204496/1‑93 (MQ=255)
              ggTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGt                                                                           <  1:3825738/93‑1 (MQ=255)
               gTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTg                                                                          >  2:5004962/1‑93 (MQ=255)
                tGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAgggg                                                                         >  2:5846871/1‑93 (MQ=255)
                tGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTgg                                                                         >  1:3676322/1‑93 (MQ=255)
                tGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTgg                                                                         <  1:4877771/93‑1 (MQ=255)
                 gATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGt                                                                        >  1:1971746/1‑93 (MQ=255)
                 gATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTAACATGGCGAAAGAAGTGGt                                                                        >  2:1252779/1‑93 (MQ=255)
                  aTCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTa                                                                       >  2:3004865/1‑93 (MQ=255)
                    cACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGt                                                                     <  1:3871146/93‑1 (MQ=255)
                     aCCCCGGTCTTCAAGCCAATTGGCGTGCATGATGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTg                                                                    <  1:5685038/93‑1 (MQ=255)
                      ccccGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTgg                                                                   >  2:950790/1‑93 (MQ=255)
                        ccGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGt                                                                 <  1:6133846/93‑1 (MQ=255)
                        ccGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGt                                                                 >  1:3717854/1‑93 (MQ=255)
                         cGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTa                                                                >  1:4623975/1‑93 (MQ=255)
                         cGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTa                                                                <  1:4834689/93‑1 (MQ=255)
                           gTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACg                                                              >  1:1125894/1‑93 (MQ=255)
                            tCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGc                                                             >  2:5333960/1‑93 (MQ=255)
                              ttCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTc                                                           >  2:5595481/1‑93 (MQ=255)
                                cAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCaa                                                         <  1:1342999/93‑1 (MQ=255)
                                 aaGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAac                                                        <  1:1048761/93‑1 (MQ=255)
                                 aaGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAac                                                        >  2:5115455/1‑93 (MQ=255)
                                  aGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAaca                                                       <  2:5923730/93‑1 (MQ=255)
                                  aGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAaca                                                       <  2:4839013/93‑1 (MQ=255)
                                  aGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAaca                                                       <  2:623982/93‑1 (MQ=255)
                                   gCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAacac                                                      <  2:2295833/93‑1 (MQ=255)
                                   gCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAacac                                                      >  1:976627/1‑93 (MQ=255)
                                   gCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAacac                                                      >  1:3455170/1‑93 (MQ=255)
                                    ccAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACAcc                                                     >  1:4329453/1‑93 (MQ=255)
                                    ccAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACAcc                                                     <  1:2425474/93‑1 (MQ=255)
                                     cAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACAccc                                                    >  1:3226513/1‑93 (MQ=255)
                                      aaTTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCct                                                   >  1:2326805/1‑93 (MQ=255)
                                      aaTTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCct                                                   <  2:4308881/93‑1 (MQ=255)
                                       aTTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCctc                                                  >  1:4703488/1‑93 (MQ=255)
                                          ggCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTac                                               <  2:5026821/93‑1 (MQ=255)
                                              tGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGCAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCg                                           >  2:5291993/1‑93 (MQ=255)
                                               gCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGc                                          <  2:2464878/93‑1 (MQ=255)
                                                cATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCa                                         >  1:188214/1‑93 (MQ=255)
                                                 aTGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAg                                        >  1:78089/1‑93 (MQ=255)
                                                     aGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGaaaa                                    >  1:2211504/1‑93 (MQ=255)
                                                      gATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAat                                   <  2:4155016/93‑1 (MQ=255)
                                                       aTAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAata                                  <  1:3723839/93‑1 (MQ=255)
                                                        tAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAatat                                 >  2:91062/1‑93 (MQ=255)
                                                         aaCTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATAtt                                <  2:1971746/93‑1 (MQ=255)
                                                         aaCTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATAtt                                <  2:1318235/93‑1 (MQ=255)
                                                          aCTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGAGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTc                               >  1:2523319/1‑93 (MQ=255)
                                                           cTGGCAGGCAACGGTTGTCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCa                              <  2:1542839/93‑1 (MQ=255)
                                                           cTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCa                              >  2:4234993/1‑93 (MQ=255)
                                                             ggcaggcaACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAgg                            >  2:2501289/1‑93 (MQ=255)
                                                             ggcaggcaACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAgg                            >  1:5311507/1‑93 (MQ=255)
                                                             ggcaggcaACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAgg                            <  2:5432945/93‑1 (MQ=255)
                                                             ggcaggcaACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAgg                            <  2:1127400/93‑1 (MQ=255)
                                                             ggcaggcaACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAgg                            <  2:2025784/93‑1 (MQ=255)
                                                             ggcaggcaACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAgg                            >  2:444616/1‑93 (MQ=255)
                                                              gcaggcaACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGa                           <  1:5838505/93‑1 (MQ=255)
                                                                aggcaACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACg                         >  2:6053920/1‑93 (MQ=255)
                                                                aggcaACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACg                         >  1:5971195/1‑93 (MQ=255)
                                                                 ggcaACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGa                        >  2:1134337/1‑93 (MQ=255)
                                                                  gcaACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGaa                       >  2:4864567/1‑93 (MQ=255)
                                                                  gcaACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGaa                       <  2:4011142/93‑1 (MQ=255)
                                                                   caACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAag                      >  2:3721456/1‑93 (MQ=255)
                                                                    aaCGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAaga                     >  1:1805393/1‑93 (MQ=255)
                                                                    aaCGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAaga                     >  1:5361401/1‑93 (MQ=255)
                                                                     aCGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAagag                    >  1:4998630/1‑93 (MQ=255)
                                                                      cGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGt                   <  2:4685644/93‑1 (MQ=255)
                                                                       ggTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGtt                  <  2:1467709/93‑1 (MQ=255)
                                                                        gTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTc                 >  1:1177072/1‑93 (MQ=255)
                                                                        gTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTc                 >  1:2835031/1‑93 (MQ=255)
                                                                        gTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTc                 <  1:2605629/93‑1 (MQ=255)
                                                                          tGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCaa               <  2:1378701/93‑1 (MQ=255)
                                                                          tGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCaa               <  2:32800/93‑1 (MQ=255)
                                                                          tGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCaa               <  2:2342941/93‑1 (MQ=255)
                                                                            gCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATc             <  1:5945753/93‑1 (MQ=255)
                                                                            gCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATc             <  2:495917/93‑1 (MQ=255)
                                                                            gCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATc             <  1:2805601/93‑1 (MQ=255)
                                                                              cTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCg           <  1:3522088/93‑1 (MQ=255)
                                                                              cTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCg           <  2:5133257/93‑1 (MQ=255)
                                                                               tGTTTACAGTTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCgg          <  1:1748429/93‑1 (MQ=255)
                                                                               tGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCgg          >  2:3304569/1‑93 (MQ=255)
                                                                                 tttACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCa        >  1:860925/1‑93 (MQ=255)
                                                                                    aCAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACCCCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGaa     <  1:2797055/93‑1 (MQ=255)
                                                                                    aCAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGaa     <  2:5452693/93‑1 (MQ=255)
                                                                                    aCAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGaa     <  1:3044036/93‑1 (MQ=255)
                                                                                    aCAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGaa     <  2:381342/93‑1 (MQ=255)
                                                                                     cAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGAAt    <  1:428428/93‑1 (MQ=255)
                                                                                      aGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGAAtt   >  1:1516316/1‑93 (MQ=255)
                                                                                       ggTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGAAttt  >  1:1268023/1‑93 (MQ=255)
                                                                                       ggTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGAAttt  <  1:1922581/93‑1 (MQ=255)
                                                                                       ggTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGAAttt  <  1:6262966/93‑1 (MQ=255)
                                                                                          |                                                                                         
GCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGAATTT  >  NC_000913/3539684‑3539863

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