Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 556,640 T→A 11.2% intergenic (‑87/+21) ybcI ← / ← ybcJ DUF457 family inner membrane protein/ribosome‑associated protein; putative RNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913556,6400TA11.2% 160.1 / 6.8 60intergenic (‑87/+21)ybcI/ybcJDUF457 family inner membrane protein/ribosome‑associated protein; putative RNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (24/29);  new base A (4/3);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGATTTTCGC  >  NC_000913/556503‑556786
                                                                                                                                         |                                                                                                                                                  
cccTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGGCTGAAAAGAGAACCGGCCTGTTGATACAGGCCg                                                                                                                                         <  5:421705/149‑1 (MQ=255)
 ccTAACCCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGCGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCgg                                                                                                                                        <  7:338083/149‑1 (MQ=255)
 ccTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCgg                                                                                                                                        <  8:379812/149‑1 (MQ=255)
 ccTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCgg                                                                                                                                        <  7:230479/149‑1 (MQ=255)
 agtAAACCAATGCATAGGGGAACTGCGGCGTGCGTAATATCGGTTGTCATGATGGTCGCTTCGGCAAAATTTCTATGCTATCAGCATGGATTAACGGGGCGTAGAGGGCAAAAGTCTGAAATGAGAACCGGCCTGTATCAACAGGCCgg                                                                                                                                        <  3:164617/147‑1 (MQ=255)
       ccAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCGTGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAgg                                                                                                                                  >  6:18212/1‑149 (MQ=255)
        cAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAgg                                                                                                                                  >  6:344591/1‑148 (MQ=255)
        cAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTATCAACAGGCCggttctctt                                                                                                                                 >  5:434615/1‑142 (MQ=255)
         aaTGCTTAGGGGAACGGCGGCGTGCGTAAGAACGGTTGGCATGATGGTCTCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACTGGGCGTAGAGGGCAAAAGTCTGAAAAGAGACCCGGCCTGTTGATACAGGCCGGGAAAGGGa                                                                                                                                <  6:348814/149‑1 (MQ=255)
             cATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCa                                                                                                                             <  1:9758/148‑1 (MQ=255)
              aTAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAgg                                                                                                                           <  3:280441/149‑1 (MQ=255)
                  gggAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGcaac                                                                                                                       <  4:2918/149‑1 (MQ=255)
                  gggAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGcaa                                                                                                                        >  3:2918/1‑148 (MQ=255)
                       cggcggCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCt                                                                                                                  >  4:183760/1‑149 (MQ=255)
                                gcgtAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCaacaac                                                                                                                    <  5:99771/138‑1 (MQ=255)
                                gcgtAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCaacaac                                                                                                                    >  6:99771/1‑138 (MQ=255)
                                  gtAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGa                                                                                                       >  5:4621/1‑149 (MQ=255)
                                   taataaCGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGAc                                                                                                      <  5:344591/149‑1 (MQ=255)
                                    aataaCGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATc                                                                                                                              >  7:328566/1‑124 (MQ=255)
                                    aataaCGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATc                                                                                                                              <  8:328566/124‑1 (MQ=255)
                                           gTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCaaaa                                                                                              <  1:383999/149‑1 (MQ=255)
                                               gCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCGTGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCa                                                                                          <  5:18212/149‑1 (MQ=255)
                                                cATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACACTGTGACCTGCAAAACTCAc                                                                                         >  2:13173/1‑149 (MQ=255)
                                                 atgatgGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACt                                                                                        >  4:176333/1‑149 (MQ=255)
                                                 atgatgCTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGAAGAGGGAAAAAGTCTGAAAACAGAACCGGCCTGTAGAACCAGGCCGGGAAAGGGAACAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACt                                                                                        <  3:269601/149‑1 (MQ=255)
                                                    atgGTCGCTTCGCGAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGACGGCAAAAGTCTGAAAAGAGAACCGGCCTGCTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGt                                                                                      <  6:13854/148‑1 (MQ=255)
                                                     tgGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGGGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCt                                                                                    <  3:249800/149‑1 (MQ=255)
                                                      ggTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCt                                                                                    >  4:224281/1‑148 (MQ=255)
                                                           cTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCgg                                                                              >  6:29666/1‑149 (MQ=255)
                                                              cGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGcga                                                                           >  3:299390/1‑149 (MQ=255)
                                                                 cAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGcgacga                                                                        <  2:162251/149‑1 (MQ=255)
                                                                  aaaaTGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTCTGACCTGCAAAACTCACTGTCTGACCGGCGACGAt                                                                       >  2:69332/1‑149 (MQ=255)
                                                                    aaTGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACACCCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGAttt                                                                     >  8:354957/1‑149 (MQ=255)
                                                                        tCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGa                                                                                  <  6:392600/132‑1 (MQ=255)
                                                                        tCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGa                                                                                  >  5:392600/1‑132 (MQ=255)
                                                                                          aTGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAAc                                               <  6:200995/149‑1 (MQ=255)
                                                                                            gAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCg                                             >  4:345699/1‑149 (MQ=255)
                                                                                              aCGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGtt                                                    >  6:104792/1‑140 (MQ=255)
                                                                                              aCGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGtt                                                    <  5:104792/140‑1 (MQ=255)
                                                                                              aCGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCa                                           >  5:189286/1‑149 (MQ=255)
                                                                                                        gagGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACt                                  <  3:183760/148‑1 (MQ=255)
                                                                                                          gggCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTc                               <  3:176333/149‑1 (MQ=255)
                                                                                                           ggCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGAc                                   <  1:469190/144‑1 (MQ=255)
                                                                                                           ggCAAAAGTCTGAAAACAGAACCGGCATGTAGATACAGGCCGGGAAAGGGAACAGGCAACAACCTGAACGCTGTGACCTGAAAAAATCAATGTCTGACCGGCGACGATTTTGAAGCGTTAGCGCGATTCAACCGCACCGTCGAcaacat                              >  2:469190/1‑144 (MQ=255)
                                                                                                            gCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCAc                             <  4:299390/149‑1 (MQ=255)
                                                                                                              aaaaGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCt                           >  6:239144/1‑149 (MQ=255)
                                                                                                                aaGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGCCCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTgg                         >  2:334013/1‑149 (MQ=255)
                                                                                                                      gAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTCAACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGCCGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCg                   >  6:328638/1‑149 (MQ=255)
                                                                                                                       aaaaGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCg                   <  5:29666/148‑1 (MQ=255)
                                                                                                                       aaaaGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACACTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCg                   <  1:13173/148‑1 (MQ=255)
                                                                                                                           gagaACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAAt              <  5:239144/149‑1 (MQ=255)
                                                                                                                           ctttcccGGCCTGTATCAACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAAt              >  4:34328/6‑149 (MQ=255)
                                                                                                                             gaACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATc             <  6:4621/148‑1 (MQ=255)
                                                                                                                               gatatcgtGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATcgc           <  3:140402/141‑1 (MQ=255)
                                                                                                                                ccGGCCTGTATCAACAGGCCGGTAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATc             <  3:52824/145‑1 (MQ=255)
                                                                                                                                ccGGCCTGTATCAACAGGCCGGTAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATc             >  4:52824/1‑145 (MQ=255)
                                                                                                                                 cGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGAt        >  2:234946/1‑149 (MQ=255)
                                                                                                                                   gCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGAttt      <  7:354957/149‑1 (MQ=255)
                                                                                                                                       gTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGATTTTCGc  >  8:441826/1‑149 (MQ=255)
                                                                                                                                       gTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGATTTTCGc  >  3:72601/1‑149 (MQ=255)
                                                                                                                                        ttGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGATTTTCGc  <  3:224281/148‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                  
CCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGATTTTCGC  >  NC_000913/556503‑556786

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