Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 556,642 A→C 8.1% intergenic (‑89/+19) 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,6420AC 8.1% 173.6 / 7.0 61intergenic (‑89/+19)ybcI/ybcJDUF457 family inner membrane protein/ribosome‑associated protein; putative RNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (25/31);  new base C (3/2);  total (28/33)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

CCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGATTTTCGCCTGC  >  NC_000913/556503‑556790
                                                                                                                                           |                                                                                                                                                    
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                 >  4:52824/1‑145 (MQ=255)
                                                                                                                                ccGGCCTGTATCAACAGGCCGGTAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATc                 <  3:52824/145‑1 (MQ=255)
                                                                                                                                 cGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGAt            >  2:234946/1‑149 (MQ=255)
                                                                                                                                   gCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGAttt          <  7:354957/149‑1 (MQ=255)
                                                                                                                                       gTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGATTTTCGc      >  3:72601/1‑149 (MQ=255)
                                                                                                                                       gTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGATTTTCGc      >  8:441826/1‑149 (MQ=255)
                                                                                                                                        ttGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGATTTTCGc      <  3:224281/148‑1 (MQ=255)
                                                                                                                                           aTACAGGCCGGGAAAGGGATCCGGCAACAACCTGTCAGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTGTCAACCTCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGATTTTCGCCTgc  >  5:339845/1‑149 (MQ=255)
                                                                                                                                           aTACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGATTTTCGCCTgc  >  1:108557/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                    
CCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGATTTTCGCCTGC  >  NC_000913/556503‑556790

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: