Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,639,254 C→T 6.7% Q10Q (CAG→CAA)  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,2540CT 6.7% 186.3 / 3.4 60Q10Q (CAG→CAA) yfgMancillary SecYEG translocon subunit; putative anti‑RcsB factor
Reads supporting (aligned to +/- strand):  ref base C (29/27);  new base T (2/2);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-01

AGGCAAGAGAAGCGGAGCGTGCAGAATCAACCTGATGGCTGTTCCAGTAGCGCCAGCCAATCAGTGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCAC  >  NC_000913/2639115‑2639402
                                                                                                                                           |                                                                                                                                                    
aGGCAAGAGAAGCGGAGCGTGCAGAATCAACCTGATGGCTGTTCCAGTAGCGCCAGCCAATCAGTGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGttt                                                                                                                                             <  3:400308/149‑1 (MQ=255)
  gCAAGAGAAGCGGAGCGTGCAGAATCAACCTGATGGCTGTTCCAGTAGCGCCAGCCAATCAGTGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTc                                                                                                                                           <  2:369701/149‑1 (MQ=255)
       agaAGCGGAGCGTGCAGAATCAACCTGATGGCTGTTCCAGTAGCGCCAGCCAATCAGTGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTtct                                                                                                                                      >  3:492115/1‑149 (MQ=255)
            cGGAGCGTGCAGAATCAACCTGATGGCTGTTCCAGTAGCGCCAGCCAATCAGTGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCg                                                                                                                                    >  1:196900/1‑146 (MQ=255)
            cGGAGCGTGCAGAATCAACCTGATGGCTGTTCCAGTAGCGCCAGCCAATCAGTGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCg                                                                                                                                    <  2:196900/146‑1 (MQ=255)
                          tCAACCTGATGGCTGTTCCAGTAGCGCCAGCCAATCAGTGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGtc                                                                                                                   <  4:264194/149‑1 (MQ=255)
                                   tGGCTGTTCCAGTAGCGCCAGCCAATCAGTGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGtccttctcctt                                                                                                          <  3:198867/149‑1 (MQ=255)
                                             aGTAGCGCCAGCCAATCAGTGCGCCAACGCCCAAAATCACACCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCa                                                                                                <  1:853443/149‑1 (MQ=255)
                                               tAGCGCCAGCCAATCAGTGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGt                                                                                              >  3:45288/1‑149 (MQ=255)
                                                    ccagccaATCAGTGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGt                                                                                         >  2:157134/1‑149 (MQ=255)
                                                           aTCAGCGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCaaa                                                                                  <  1:356301/149‑1 (MQ=255)
                                                              aGTGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATgc                                                                               <  3:377070/149‑1 (MQ=255)
                                                              aGTGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATgc                                                                               >  3:830510/1‑149 (MQ=255)
                                                              aGTGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATgc                                                                               >  4:377070/1‑149 (MQ=255)
                                                                tGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATgcgc                                                                             <  3:229542/149‑1 (MQ=255)
                                                                            cAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATc                                                                 >  1:461233/1‑149 (MQ=255)
                                                                              aaaTCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATgcgc                                                                             <  4:136013/135‑1 (MQ=255)
                                                                              aaaTCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATgcgc                                                                             >  3:136013/1‑135 (MQ=255)
                                                                               aaTCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTg                                                              <  1:588675/149‑1 (MQ=255)
                                                                                aTCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTgc                                                             >  2:1360/1‑149 (MQ=255)
                                                                                   aCCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCa                                                          >  1:865525/1‑149 (MQ=255)
                                                                                    ccccAACAGCCAGTGCTTTGCCATTTTCAGCGAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCaa                                                         >  2:207430/1‑149 (MQ=255)
                                                                                    ccccAACAGCCAGTGCTTTGCCATTTTCAGCAATAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCaa                                                         >  3:750495/1‑149 (MQ=255)
                                                                                     cccAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAAttt                                                                                                                             <  2:623953/80‑1 (MQ=255)
                                                                                     cccAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAAttt                                                                                                                             >  1:623953/1‑80 (MQ=255)
                                                                                     cccAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAAc                                                        <  1:170874/149‑1 (MQ=255)
                                                                                      ccAACAACCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACt                                                       >  3:419175/1‑149 (MQ=255)
                                                                                       cAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTAGCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTg                                                      >  4:478601/1‑149 (MQ=255)
                                                                                       cAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTg                                                      >  2:231070/1‑149 (MQ=255)
                                                                                           aGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGt                                                  >  3:64676/1‑149 (MQ=255)
                                                                                             ccAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGttt                                                <  4:853875/149‑1 (MQ=255)
                                                                                               aGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGc                                              >  1:375929/1‑149 (MQ=255)
                                                                                               aGTGCTTTGCCATTTTCAACAAAAAAGCGTTTGACCGCTTCTACCTGGTCCTTTTCGTTCTCGTAAATTTCCACGCTGACCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGc                                              <  2:565328/149‑1 (MQ=255)
                                                                                                    tttGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACTTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATgcgc                                                                             >  3:341259/1‑113 (MQ=255)
                                                                                                    tttGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACTTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATgcgc                                                                             <  4:341259/113‑1 (MQ=255)
                                                                                                    tttGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACTTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATgcgc                                                                             <  4:120850/113‑1 (MQ=255)
                                                                                                    tttGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACTTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATgcgc                                                                             >  3:120850/1‑113 (MQ=255)
                                                                                                    tttGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCAcc                                         >  2:142092/1‑149 (MQ=255)
                                                                                                      tGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCag                                       >  3:121439/1‑149 (MQ=255)
                                                                                                        ccATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCagag                                     >  1:729175/1‑149 (MQ=255)
                                                                                                         cATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAgc                                    <  1:446081/149‑1 (MQ=255)
                                                                                                          aTTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAgcg                                   >  1:44963/1‑149 (MQ=255)
                                                                                                           ttttCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAgcgc                                  >  3:292271/1‑149 (MQ=255)
                                                                                                                aGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATc                             >  2:480781/1‑149 (MQ=255)
                                                                                                                 gCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATcc                            >  3:363378/1‑149 (MQ=255)
                                                                                                                 gCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATcc                            >  4:371309/1‑149 (MQ=255)
                                                                                                                   aaaaaaGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCtt                          >  3:478998/1‑149 (MQ=255)
                                                                                                                     aaaaGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCa                        <  3:693394/149‑1 (MQ=255)
                                                                                                                     aaaaGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCa                        <  3:789114/149‑1 (MQ=255)
                                                                                                                           gTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTActg                  >  4:329645/1‑149 (MQ=255)
                                                                                                                               aaCCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCtgt              >  1:540460/1‑149 (MQ=255)
                                                                                                                               aaCCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCtgt              <  4:154443/149‑1 (MQ=255)
                                                                                                                               aaCCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCtgt              <  4:492115/149‑1 (MQ=255)
                                                                                                                               aaCCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCtgt              <  4:587832/149‑1 (MQ=255)
                                                                                                                                     ttCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGtt        <  4:292271/149‑1 (MQ=255)
                                                                                                                                     ttCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTACTCACCAGAGCGCAAATCCTTCACTGCTGCTGTGCCGtt        <  4:766980/149‑1 (MQ=255)
                                                                                                                                      tCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTa       <  3:371309/149‑1 (MQ=255)
                                                                                                                                      tCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTa       <  4:45288/149‑1 (MQ=255)
                                                                                                                                        tACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGc     <  3:478601/149‑1 (MQ=255)
                                                                                                                                        tACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGc     <  3:564967/149‑1 (MQ=255)
                                                                                                                                           cTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCAc  >  4:528232/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                    
AGGCAAGAGAAGCGGAGCGTGCAGAATCAACCTGATGGCTGTTCCAGTAGCGCCAGCCAATCAGTGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCAC  >  NC_000913/2639115‑2639402

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑