Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,639,252 A→T 6.3% V11E (GTA→GAA)  yfgM ← ancillary SecYEG translocon subunit; putative anti‑RcsB factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,639,2520AT 6.3% 194.1 / 3.8 64V11E (GTA→GAA) yfgMancillary SecYEG translocon subunit; putative anti‑RcsB factor
Reads supporting (aligned to +/- strand):  ref base A (30/30);  new base T (2/2);  total (32/32)
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.75e-01

GATAGGCAAGAGAAGCGGAGCGTGCAGAATCAACCTGATGGCTGTTCCAGTAGCGCCAGCCAATCAGTGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTA  >  NC_000913/2639112‑2639397
                                                                                                                                            |                                                                                                                                                 
gATAGGCAAGAGAAGCGGACCGTGCAGAATCAACCTGATGGCTGTTCCACTAGAGCCAGCCAATCAGTGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTTCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTc                                                                                                                                            <  8:71553/148‑1 (MQ=255)
 aTAGGCAAGAGAAGCGGAGCGTGCAGAATCAACCTGATGGCTGTTCCAGTAGCGCCAGCCAATCAGTGCGCCACCGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGt                                                                                                                                          >  7:328390/1‑149 (MQ=255)
 aTAGGCAAGAGAAGCGGAGCGTGCAGAATCAACCTGATGGCTGTTCCAGTAGCGCCAGCCAATCAGTGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGt                                                                                                                                          >  5:1465/1‑149 (MQ=255)
 aTAGGCAAGAGAAGCGGAGCGTGCAGAATCAACCTGATGGCTGTTCCAGTAGCGCCAGCCAATCAGTGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGt                                                                                                                                          >  6:362373/1‑149 (MQ=255)
        agagAAGCGGAGCGTGCAGAATCAACCTGATGGCTGTTCCAGTAGCGCCAGCCAATCAGTGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGtt                                                                                                                                   >  1:104215/1‑149 (MQ=255)
        agagAAGCGGAGCGTGCAGAATCAACCTGATGGCTGTTCCAGTAGCGCCAGCCAATCAGTGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGtt                                                                                                                                   <  2:360056/149‑1 (MQ=255)
         gagaAGCGGAGCGTGCAGAATCAACCTGATGGCTGTTCCAGTAGCGCCAGCCAATCAGTGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTtc                                                                                                                                  >  6:176843/1‑149 (MQ=255)
                   gCGTGCAGAATCAACCTGATGGCTGTTCCAGTAGCGCCAGCCAATCAGTGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAAttt                                                                                                                        <  6:1465/149‑1 (MQ=255)
                         aGAATCAACCTGATGGCTGTTCCAGTAGCGCCAGCCAATCAGTGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGc                                                                                                                  >  7:345892/1‑149 (MQ=255)
                             tCAACCTGATGGCTGTTCCAGTAGCGCCAGCCAATCAGTGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGtc                                                                                                              >  1:168396/1‑149 (MQ=255)
                                 ccTGATGGCTGTTCCAGTAGCGCCAGCCAATCAGTGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGtccttc                                                                                                          <  8:345892/149‑1 (MQ=255)
                                 ccTGATGGCTGTTCCAGGAGCGCCAGCCAATCAGTGCGCCAACTCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGtccttc                                                                                                          <  8:328390/149‑1 (MQ=255)
                                   tGATGGCTGTTCCAGTAGCGCCAGCCAATCAGTGCGCCAACGCCCAAAACCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGtccttctc                                                                                                        >  5:108147/1‑149 (MQ=255)
                                     aTGGCTGTTCCAGTAGCGCCAGCCAATCAGTGCGCCAACGCCCAAAATCACCCCAACAGCCAGGGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGtccttctc                                                                                                        <  1:316073/147‑1 (MQ=255)
                                      tGGCTGTTCCAGTAGCGCCAGCCAATCAGTGCGCAAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGATTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGtccttctcctt                                                                                                     <  3:280009/149‑1 (MQ=255)
                                            ttCCAGTAGCGCCAGCCAATCAGTGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTaa                                                                                               >  6:79725/1‑149 (MQ=255)
                                               cAGTAGCGCCAGCCAATCAGTGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAAcc                                                                                             >  2:334743/1‑148 (MQ=255)
                                               cAGTAGCGCCAGCCAATCAGAGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAAcc                                                                                             <  4:37090/148‑1 (MQ=255)
                                                 gTAGCGCCAGCCAATCAGTGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAg                                                                                          <  2:168396/149‑1 (MQ=255)
                                                             aaTCAGTGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTTCCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCa                                                                               >  1:269988/1‑148 (MQ=255)
                                                             aaTCAGTGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTGTAACAGCTTCTTCCTGGGCGTGTTCGTTCTCGTAAATTTCCACGCTGTCCTACTCATTCCTTAACCCAGTAACGTGCGCa                                                                               >  1:269992/1‑148 (MQ=255)
                                                              aTCAGTGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCaaa                                                                             >  5:253688/1‑149 (MQ=255)
                                                              aTCAGTGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCaaa                                                                             >  8:367087/1‑149 (MQ=255)
                                                               tCAGTGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAAt                                                                            <  8:264490/149‑1 (MQ=255)
                                                                     cgcCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCgg                                                                      <  7:182200/149‑1 (MQ=255)
                                                                        cAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTa                                                                   >  8:164554/1‑149 (MQ=255)
                                                                         aaCGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGtccttctccttc                                                                                                    >  3:291729/1‑115 (MQ=255)
                                                                         aaCGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGtccttctccttc                                                                                                    <  4:291729/115‑1 (MQ=255)
                                                                           cGCCCAAAACCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGc                                                                <  6:108147/149‑1 (MQ=255)
                                                                            gCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGc                                                                >  3:273036/1‑148 (MQ=255)
                                                                            gCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGc                                                                >  8:18541/1‑148 (MQ=255)
                                                                             cccAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTa                                                              <  5:176843/149‑1 (MQ=255)
                                                                              ccAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATgcgc                                                                        <  4:229730/138‑1 (MQ=255)
                                                                              ccAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGATGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTAt                                                             >  1:309151/1‑149 (MQ=255)
                                                                               cAAAATCACCCCTACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATc                                                            <  5:79725/149‑1 (MQ=255)
                                                                                aaaaTCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATcc                                                           >  8:8269/1‑149 (MQ=255)
                                                                                aaaaTCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATcc                                                           >  8:367651/1‑149 (MQ=255)
                                                                                aaaaTCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATcc                                                           <  1:334743/149‑1 (MQ=255)
                                                                                 aaaTCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCt                                                          >  1:7970/1‑149 (MQ=255)
                                                                                      aCCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTgcgc                                                      <  6:253688/148‑1 (MQ=255)
                                                                                        cccAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAAc                                                   >  8:220929/1‑149 (MQ=255)
                                                                                               gCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGtt                                            <  2:309151/149‑1 (MQ=255)
                                                                                                 cAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGttt                                           >  1:249342/1‑148 (MQ=255)
                                                                                                 cAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGttt                                           <  2:249342/148‑1 (MQ=255)
                                                                                                  aGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTg                                          >  6:367605/1‑148 (MQ=255)
                                                                                                  aGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGc                                         >  4:340487/1‑149 (MQ=255)
                                                                                                      cTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCAc                                     <  7:367651/149‑1 (MQ=255)
                                                                                                         tGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCa                                   <  5:362373/148‑1 (MQ=255)
                                                                                                         tGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCa                                   <  2:7970/148‑1 (MQ=255)
                                                                                                           ccATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCagag                                <  4:282303/149‑1 (MQ=255)
                                                                                                             aTTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTgc                                                                                  <  1:136360/97‑1 (MQ=255)
                                                                                                             aTTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTgc                                                                                  >  2:136360/1‑97 (MQ=255)
                                                                                                                   aGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGttt                                           <  2:337943/130‑1 (MQ=255)
                                                                                                                   aGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGttt                                           >  1:337943/1‑130 (MQ=255)
                                                                                                                   aGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATc                        >  8:212974/1‑149 (MQ=255)
                                                                                                                      aaaaaaGCGTTTAACCGCTTCTTCCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCtt                     <  2:269988/149‑1 (MQ=255)
                                                                                                                      aaaaaaGCGTTTAACCGCTTCTTACTTGTCGTTTTAGTTCTCGTAAATTTCCACGCTGTCATGCTCCTTCCTTACCCCAGTAAAGTGCGCAACTGCGCGGCTACGCTATCCTGCGCAACTGCCGGTTGCTCACCAGAGCGCAAATCCtt                     <  2:269992/149‑1 (MQ=255)
                                                                                                                       aaaaaGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTc                    >  7:106548/1‑149 (MQ=255)
                                                                                                                               tttAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTActgc            >  8:284382/1‑149 (MQ=255)
                                                                                                                                  aaCCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCATTCCTTAACCCAGTAACGTGCGCAAATGCGCTGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCtgt         >  3:278134/1‑149 (MQ=255)
                                                                                                                                   aCCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCtgt         <  8:106548/148‑1 (MQ=255)
                                                                                                                                          cTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTa  <  7:8269/148‑1 (MQ=255)
                                                                                                                                          cTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTa  <  7:284382/148‑1 (MQ=255)
                                                                                                                                          cTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTa  <  7:164554/148‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                 
GATAGGCAAGAGAAGCGGAGCGTGCAGAATCAACCTGATGGCTGTTCCAGTAGCGCCAGCCAATCAGTGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTA  >  NC_000913/2639112‑2639397

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑