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,7740GA97.1% 470.4 / ‑0.4 139A578T (GCC→ACC) feoBferrous iron transport protein B
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base A (60/74);  total (62/77)
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.85e-01

Only 100 of 140 total aligned reads displayed.

GCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGAATTT  >  NC_000913/3539684‑3539863
                                                                                          |                                                                                         
gCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTAcc                                                                                         <  1:4723422/93‑1 (MQ=255)
  gtccgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCAt                                                                                       <  1:2477373/93‑1 (MQ=255)
  gtccgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCAt                                                                                       >  2:437681/1‑93 (MQ=255)
   tccgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATg                                                                                      >  1:119760/1‑93 (MQ=255)
   tccgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATg                                                                                      <  1:1276386/93‑1 (MQ=255)
   gccgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATg                                                                                      <  2:4577572/92‑1 (MQ=255)
   gccgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATg                                                                                      <  1:94784/92‑1 (MQ=255)
    ccgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATgg                                                                                     <  1:3142034/93‑1 (MQ=255)
    ccgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATgg                                                                                     >  1:4047743/1‑93 (MQ=255)
     cgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGTTAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGc                                                                                    <  1:1406170/93‑1 (MQ=255)
     cgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGc                                                                                    >  2:2160033/1‑93 (MQ=255)
       tcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGa                                                                                  <  1:1271346/93‑1 (MQ=255)
         aGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCggag                                                                                >  2:1477434/1‑90 (MQ=255)
         aGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGaaa                                                                                >  2:4498432/1‑93 (MQ=255)
          gCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAaag                                                                               <  1:3981852/93‑1 (MQ=255)
            cGGGTGATCCCCCCGGTCTGCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAaagaa                                                                             <  2:1981116/93‑1 (MQ=255)
             gggTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAaagaag                                                                            >  1:949351/1‑93 (MQ=255)
             gggTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAaagaag                                                                            >  1:4669130/1‑93 (MQ=255)
                tGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTgg                                                                         <  2:1387674/93‑1 (MQ=255)
                 gaccACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGt                                                                        <  1:1178295/90‑1 (MQ=255)
                 gATCCCCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGt                                                                        <  2:1888802/93‑1 (MQ=255)
                 gATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGt                                                                        <  2:2647882/93‑1 (MQ=255)
                 gATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGt                                                                        <  2:1888801/93‑1 (MQ=255)
                   tcgccccGGTCTACAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAg                                                                      <  1:3506682/90‑1 (MQ=255)
                    cACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGt                                                                     <  2:4719654/93‑1 (MQ=255)
                    cACCCCGGTCTTCAAGCCAATTGCCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGt                                                                     <  2:2589853/93‑1 (MQ=255)
                     aCCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTg                                                                    >  2:356936/1‑93 (MQ=255)
                       cccGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTggg                                                                  <  2:2059808/93‑1 (MQ=255)
                         cGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTa                                                                >  1:4527716/1‑93 (MQ=255)
                          ggTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTAc                                                               <  1:1979542/93‑1 (MQ=255)
                           gTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACg                                                              <  1:2976174/93‑1 (MQ=255)
                              ttCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTc                                                           >  1:4117923/1‑93 (MQ=255)
                                cAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCaa                                                         <  1:3497367/93‑1 (MQ=255)
                                 aaGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAac                                                        <  1:2682420/93‑1 (MQ=255)
                                   gCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAacac                                                      >  2:1365776/1‑93 (MQ=255)
                                   gCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGCGGTAGTGGGTACGCTCAacac                                                      >  2:371360/1‑93 (MQ=255)
                                    ccAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACAcc                                                     <  1:4040140/93‑1 (MQ=255)
                                     cAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACAccc                                                    <  2:1255484/93‑1 (MQ=255)
                                     cAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACAccc                                                    <  1:42435/93‑1 (MQ=255)
                                      aaTTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCct                                                   >  1:2336000/1‑93 (MQ=255)
                                       aTTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCctc                                                  >  1:1521063/1‑93 (MQ=255)
                                       aTTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCctc                                                  <  2:4668734/93‑1 (MQ=255)
                                        ttGGCGTGCATGAAGATAACTGGCATGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCctct                                                 >  1:4533159/1‑93 (MQ=255)
                                        ttGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCctct                                                 <  2:1089461/93‑1 (MQ=255)
                                        ttGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCctct                                                 >  1:1941784/1‑93 (MQ=255)
                                        ttGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCctct                                                 >  1:276337/1‑93 (MQ=255)
                                           gCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTaca                                              <  2:4355929/93‑1 (MQ=255)
                                            cGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTacac                                             <  2:3123997/93‑1 (MQ=255)
                                              tGCATGAAGATAACTGGCAGGCAGCGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCg                                           <  2:250385/93‑1 (MQ=255)
                                               gCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTCAACCga                                          >  1:4038661/1‑92 (MQ=255)
                                                 aTGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAg                                        >  1:685404/1‑93 (MQ=255)
                                                 aTGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAg                                        <  2:4582538/93‑1 (MQ=255)
                                                  tGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGa                                       >  1:4748276/1‑93 (MQ=255)
                                                  tGAAGATAACTGGCAGGCAACGGTCGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGa                                       >  1:1839048/1‑93 (MQ=255)
                                                   gAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGaa                                      >  1:4749496/1‑93 (MQ=255)
                                                   gAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGaa                                      >  2:1218758/1‑93 (MQ=255)
                                                       aTAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAata                                  <  1:2785708/93‑1 (MQ=255)
                                                        tAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAatat                                 >  1:2667248/1‑93 (MQ=255)
                                                          aCTTGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTc                               <  1:2325742/93‑1 (MQ=255)
                                                          aCTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATAtta                               >  1:3639421/1‑92 (MQ=255)
                                                          aCTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTc                               <  2:1054784/93‑1 (MQ=255)
                                                           cTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCa                              <  1:2294578/93‑1 (MQ=255)
                                                           cTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCa                              >  1:152232/1‑93 (MQ=255)
                                                           cTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCa                              <  2:4679736/93‑1 (MQ=255)
                                                           cTGGCAGGCAACGGTTGGCCTGTTTACAGGGCCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCa                              <  1:2294579/93‑1 (MQ=255)
                                                            tGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAg                             >  2:148858/1‑93 (MQ=255)
                                                             ggcaggcaACGGTTGGCCTGTTTACAGGTACCATGGGGAACGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCTCAGAAAATATTCAgg                            >  2:1063032/1‑93 (MQ=255)
                                                             ggcaggcaACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAgg                            >  1:4381433/1‑93 (MQ=255)
                                                               caggcaACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAggct                          >  2:3408346/1‑91 (MQ=255)
                                                               caggcaACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGAc                          <  2:743673/93‑1 (MQ=255)
                                                               caggcaACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGAc                          <  2:1788799/93‑1 (MQ=255)
                                                                 ggcaACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGa                        >  1:3628175/1‑93 (MQ=255)
                                                                 ggcaACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGa                        <  1:3829920/93‑1 (MQ=255)
                                                                 ggcaACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGa                        >  2:4267557/1‑93 (MQ=255)
                                                                   caACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAag                      <  2:3125494/93‑1 (MQ=255)
                                                                    aaCGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAaga                     <  2:1959099/93‑1 (MQ=255)
                                                                    aaCGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAaga                     <  2:1103230/93‑1 (MQ=255)
                                                                    aaCGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAaga                     >  1:741233/1‑93 (MQ=255)
                                                                    aaCGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAaga                     <  2:2830263/93‑1 (MQ=255)
                                                                     aCGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAagag                    >  1:4623651/1‑93 (MQ=255)
                                                                      cGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGt                   >  1:4655760/1‑93 (MQ=255)
                                                                        gTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTc                 >  1:1146007/1‑93 (MQ=255)
                                                                        gTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTc                 <  2:1784672/93‑1 (MQ=255)
                                                                        gTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTc                 <  1:3674809/93‑1 (MQ=255)
                                                                          tGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCaa               >  2:3114555/1‑93 (MQ=255)
                                                                          tGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCaa               <  1:1799333/93‑1 (MQ=255)
                                                                            gCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTGAACAGCCTCTACACCGCAGAAAATATTCAGGACGAAGACGTCAATc             >  2:2434130/1‑93 (MQ=255)
                                                                            gCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATc             <  2:588498/93‑1 (MQ=255)
                                                                            gCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATc             <  2:4242870/93‑1 (MQ=255)
                                                                             ccTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATcc            >  1:899089/1‑93 (MQ=255)
                                                                             ccTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATcc            <  1:1890420/93‑1 (MQ=255)
                                                                               tGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCgg          <  2:397348/93‑1 (MQ=255)
                                                                               tGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCgg          >  1:3551101/1‑93 (MQ=255)
                                                                               tGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCgg          >  2:1194894/1‑93 (MQ=255)
                                                                                 tttACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCa        <  1:4388879/93‑1 (MQ=255)
                                                                                 tttACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCa        <  1:3307999/93‑1 (MQ=255)
                                                                                 tttACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCa        <  2:3761565/93‑1 (MQ=255)
                                                                                    aCAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGaa     <  2:1749071/93‑1 (MQ=255)
                                                                                       ggTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCTGGACGAAGAGTTCAATCCGGCAGAAttt  >  2:1086292/1‑93 (MQ=255)
                                                                                       ggTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGAAttt  <  2:2701478/93‑1 (MQ=255)
                                                                                          |                                                                                         
GCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGAATTT  >  NC_000913/3539684‑3539863

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