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% 805.7 / ‑5.5 228A578T (GCC→ACC) feoBferrous iron transport protein B
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (109/118);  total (109/119)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

Only 100 of 236 total aligned reads displayed.

GCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGAATTT  >  NC_000913/3539684‑3539863
                                                                                          |                                                                                         
gCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTAcc                                                                                         <  1:3500372/93‑1 (MQ=255)
  gtccgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCAt                                                                                       <  2:353746/93‑1 (MQ=255)
   tccgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATg                                                                                      <  1:339114/93‑1 (MQ=255)
   tccgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATg                                                                                      >  1:2408755/1‑93 (MQ=255)
      gtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCg                                                                                   <  2:461032/93‑1 (MQ=255)
       tcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGa                                                                                  >  1:2722380/1‑93 (MQ=255)
         aGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGaaa                                                                                >  1:2537112/1‑93 (MQ=255)
         aGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGaaa                                                                                >  1:603159/1‑93 (MQ=255)
         aGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGaaa                                                                                <  1:4661957/93‑1 (MQ=255)
            cGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAaagaa                                                                             <  2:1955803/93‑1 (MQ=255)
            cGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAaagaa                                                                             >  1:2099989/1‑93 (MQ=255)
                 gATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGt                                                                        <  1:1215467/93‑1 (MQ=255)
                  aTCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTa                                                                       >  2:3474622/1‑93 (MQ=255)
                   tCACCCCGGTCTTCAAGCCAGTTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAg                                                                      <  2:3895313/93‑1 (MQ=255)
                    cACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGt                                                                     >  2:1081313/1‑93 (MQ=255)
                     cccccGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTg                                                                    <  1:170008/92‑1 (MQ=255)
                     aCCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTg                                                                    >  1:4011325/1‑93 (MQ=255)
                     aCCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTg                                                                    >  2:3731719/1‑93 (MQ=255)
                      ccccGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTgg                                                                   >  1:115365/1‑93 (MQ=255)
                      ccccGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTgg                                                                   <  2:2132251/93‑1 (MQ=255)
                       cccGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTggg                                                                  <  2:2279043/93‑1 (MQ=255)
                        ccGGTTTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGGTGGACTGTTTGCAGGCACCATGGCGAAAGAAGCGGAAGTGGGt                                                                 >  2:3424737/1‑93 (MQ=255)
                        ccGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTggga                                                                 >  1:1444188/1‑92 (MQ=255)
                           gTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACg                                                              >  1:2975341/1‑93 (MQ=255)
                            tCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGc                                                             <  1:814168/93‑1 (MQ=255)
                              ttCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCta                                                           >  1:2987101/1‑92 (MQ=255)
                              ttCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTc                                                           >  2:947793/1‑93 (MQ=255)
                              ttCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTc                                                           >  2:2350459/1‑93 (MQ=255)
                              ttCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTc                                                           >  1:1413670/1‑93 (MQ=255)
                              ttCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGCTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGCACGCTc                                                           >  1:394442/1‑93 (MQ=255)
                               tCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCa                                                          >  1:754306/1‑93 (MQ=255)
                                cAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCaa                                                         >  1:3918751/1‑93 (MQ=255)
                                 aaGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAac                                                        <  2:2969074/93‑1 (MQ=255)
                                  aGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAaca                                                       >  2:1531272/1‑93 (MQ=255)
                                  aGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAaca                                                       >  2:766689/1‑93 (MQ=255)
                                   gCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAacac                                                      <  2:204448/93‑1 (MQ=255)
                                   gCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAacac                                                      <  1:4076898/93‑1 (MQ=255)
                                     cAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACAccc                                                    <  1:2220903/93‑1 (MQ=255)
                                      aaTTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCct                                                   <  1:1052871/93‑1 (MQ=255)
                                      aaTTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCct                                                   >  1:3762375/1‑93 (MQ=255)
                                      aaTTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCct                                                   >  1:1267842/1‑93 (MQ=255)
                                      aaTTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGGACGCTCAACACCct                                                   >  2:4422719/1‑93 (MQ=255)
                                        ttGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCctct                                                 >  1:3055831/1‑93 (MQ=255)
                                         tGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTa                                                <  2:1790027/93‑1 (MQ=255)
                                          ggCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTac                                               <  1:2064647/93‑1 (MQ=255)
                                           gCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTaca                                              <  2:2664070/93‑1 (MQ=255)
                                            cGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTacac                                             <  2:174228/93‑1 (MQ=255)
                                             gTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACAcc                                            >  2:3741612/1‑93 (MQ=255)
                                             gTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGCGGTAGTGGGTACGCTCAACACCCTCTACAcc                                            >  1:4427351/1‑93 (MQ=255)
                                             gTGCATGAAGATAACTGGCAGGCAACGGGTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGGAGTGGGGACGCTCAACACCCTCTACAcc                                            >  2:111411/1‑93 (MQ=255)
                                              tGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCg                                           >  1:3311034/1‑93 (MQ=255)
                                                cATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCa                                         >  2:1464172/1‑93 (MQ=255)
                                                 aTGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAg                                        >  1:108084/1‑93 (MQ=255)
                                                 aTGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAg                                        >  2:380410/1‑93 (MQ=255)
                                                  agAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGa                                       <  2:4676441/92‑1 (MQ=255)
                                                   gAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGaa                                      >  1:3585943/1‑93 (MQ=255)
                                                    gaGATAACTGGCAGTCAACTGTTGACCTGTTTACAGCTGCCATGGCGAAAGAATTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGaaa                                     <  2:1600648/92‑1 (MQ=255)
                                                    aaGATAACTGGCAGGCATCGGTTGGCCTGTTTACAGGTACCATGTCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGaaa                                     <  1:2291961/93‑1 (MQ=255)
                                                     aGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGaaaa                                    <  1:4265314/93‑1 (MQ=255)
                                                     aGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGaaaa                                    >  1:3658592/1‑93 (MQ=255)
                                                     aGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGaaaa                                    <  2:2785221/93‑1 (MQ=255)
                                                     aGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGaaaa                                    >  2:2518033/1‑93 (MQ=255)
                                                     aGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGaaaa                                    <  2:2311501/93‑1 (MQ=255)
                                                     aGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCCCAACACCCTCTACACCGCAGaaaa                                    <  2:2931329/93‑1 (MQ=255)
                                                      gATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACATCCTCTACACCGCAGAAAat                                   <  2:4526880/93‑1 (MQ=255)
                                                      gATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAat                                   <  2:3874400/93‑1 (MQ=255)
                                                       aTAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAata                                  <  1:3882046/93‑1 (MQ=255)
                                                       aTAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAata                                  >  1:4393041/1‑93 (MQ=255)
                                                       aTAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAata                                  <  2:3157492/93‑1 (MQ=255)
                                                         aaCTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATAtt                                <  2:1322723/93‑1 (MQ=255)
                                                         aaCTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATAtt                                >  2:2859662/1‑93 (MQ=255)
                                                          aCTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACACTCTACACCGCAGAAAATAtta                               >  2:3268632/1‑92 (MQ=255)
                                                           ccggcaggcaACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCa                              <  1:1325738/91‑1 (MQ=255)
                                                           cTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCa                              >  1:3327221/1‑93 (MQ=255)
                                                             gccaggctACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAgg                            <  2:592701/91‑1 (MQ=255)
                                                              gcaggcaACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGa                           <  2:708263/93‑1 (MQ=255)
                                                              gcaggcaACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGa                           <  1:2149713/93‑1 (MQ=255)
                                                              gcaggcaACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGa                           <  1:1508657/93‑1 (MQ=255)
                                                              gcaggcaACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGa                           <  1:2848171/93‑1 (MQ=255)
                                                                  gcaACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGaa                       <  2:1774643/93‑1 (MQ=255)
                                                                   caACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAag                      >  2:1893793/1‑93 (MQ=255)
                                                                    aaCGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAaga                     <  2:2161030/93‑1 (MQ=255)
                                                                     aCGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAagag                    <  2:1028265/93‑1 (MQ=255)
                                                                       ggTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGtt                  >  2:1811526/1‑93 (MQ=255)
                                                                       ggTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGtt                  >  1:4593683/1‑93 (MQ=255)
                                                                        gTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTc                 >  2:4176193/1‑93 (MQ=255)
                                                                        gTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTc                 >  1:180051/1‑93 (MQ=255)
                                                                        gTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTc                 <  2:259492/93‑1 (MQ=255)
                                                                        gTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGGTc                 >  1:1877780/1‑93 (MQ=255)
                                                                          tGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCaa               <  1:2619742/93‑1 (MQ=255)
                                                                           ggCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAAt              <  2:240826/93‑1 (MQ=255)
                                                                            gCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATc             <  2:915500/93‑1 (MQ=255)
                                                                               tGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCgg          >  1:1618195/1‑93 (MQ=255)
                                                                                gTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGc         >  1:2597395/1‑93 (MQ=255)
                                                                                 tttACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCa        >  1:100823/1‑93 (MQ=255)
                                                                                 tttACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCa        >  1:3762574/1‑93 (MQ=255)
                                                                                     cAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGAAt    <  2:4121035/93‑1 (MQ=255)
                                                                                      aGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGAAtt   <  2:1695486/93‑1 (MQ=255)
                                                                                       ggTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGAAttt  >  2:3304678/1‑93 (MQ=255)
                                                                                       ggTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGAAttt  >  1:3556750/1‑93 (MQ=255)
                                                                                          |                                                                                         
GCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGAATTT  >  NC_000913/3539684‑3539863

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: