Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,679,430 Δ1 bp 7.8% intergenic (‑63/+34) yphF ← / ← yphG putative ABC transporter periplasmic binding protein YphF/DUF5107 domain‑containing protein YphG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,679,4300G. 7.8% 172.5 / 11.9 64intergenic (‑63/+34)yphF/yphGputative ABC transporter periplasmic binding protein YphF/DUF5107 domain‑containing protein YphG
Reads supporting (aligned to +/- strand):  ref base G (33/25);  new base . (2/3);  total (35/29)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.99e-01

CTGGCAAATAGTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATAT  >  NC_000913/2679294‑2679566
                                                                                                                                        |                                                                                                                                         
cTGGCAAATAGTGCACTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTg                                                                                                                                    <  1:230883/143‑1 (MQ=255)
cTGGCAAACATTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTTTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTg                                                                                                                                    <  2:561588/143‑1 (MQ=255)
 tGGCTAATAGTGCGCTACCTAATAGTGTAGCCATTAATAATAGATTACGGGTTGTTCTCATCTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGc                                                                                                                                   <  1:184203/143‑1 (MQ=255)
 tGGCAAATAGTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAGCGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGc                                                                                                                                   <  2:426728/143‑1 (MQ=255)
 tGGCAAATAGTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGc                                                                                                                                   >  2:444519/1‑143 (MQ=255)
     agaTAGTGCGCTGCCTACTCGCGTCGCCAGTAATAATAGGTTACGGGTTGTTCTCATTTTTGTAGGCATCGAGCCGCCTGTAGGGTTTTTATTAACAACGGCTTCTGCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCt                                                                                                                               <  2:1022994/141‑1 (MQ=255)
     aaaTAGTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCACGCGCTGCCGCt                                                                                                                               <  1:202925/143‑1 (MQ=255)
     aaaTAGTGCGCTACCTAATAGCGTAGCCATTAATAAGAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCACGCGCTGCCGCt                                                                                                                               <  1:37678/143‑1 (MQ=255)
       aTAGTGCGCTACCTGATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCACGCGCTGCCGCTCa                                                                                                                             <  2:921843/143‑1 (MQ=255)
              gcTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTa                                                                                                                      <  2:826704/143‑1 (MQ=255)
              gcTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGAGGCAGAGCGTGCCGCTCAAATATTa                                                                                                                      >  1:56307/1‑143 (MQ=255)
                   ctactaGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGtt                                                                                                                 <  2:789000/143‑1 (MQ=255)
                     aaTAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACCACGGCTTATTCTAATTATTTTGTGATTAGCGGCAGCGCGTGCGGCTCCGATCTTATAGtttt                                                                                                               >  2:789978/1‑141 (MQ=255)
                     aaTAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCTAATACTACCGttca                                                                                                               >  1:427031/1‑141 (MQ=255)
                     aaTAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTACTGTAGGTGTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGttca                                                                                                               >  1:1015243/1‑141 (MQ=255)
                         gCGTAGCCATTAATAATAGATTACGGGTTGTTCTCAGTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGGTAATAgg                                                                                                           >  2:21255/1‑143 (MQ=255)
                              gCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTa                                                                                                      <  2:1131748/143‑1 (MQ=255)
                                 aTTAAGAATGGACGACTGGCTGTTCTCATTTTTGTAGGCATCGCGCCTCCCGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTg                                                                                                   <  1:293762/143‑1 (MQ=255)
                                  ttAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTgc                                                                                                  >  2:784561/1‑142 (MQ=255)
                                        ataGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCTTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCgg                                                                                            <  2:69302/143‑1 (MQ=255)
                                                cGGGTTGTTCTCAGTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAAt                                                                                    <  1:733676/143‑1 (MQ=255)
                                                   gttgttCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTGTAGTGAACGCGGTTAAAAATaac                                                                                 >  1:161264/1‑142 (MQ=255)
                                                          tcATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAAAAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATc                                                                          >  1:981079/1‑143 (MQ=255)
                                                          tcATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATCAGCGGCAGCGAGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATc                                                                          >  1:344182/1‑143 (MQ=255)
                                                                     ggCATAGGTCCTACTGTTGGGTTGATATTGACAACGTCTTATTCTAATTATTTTGGGATGGGCGGCAGTGCGTGTCGTTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGaa                                                               <  2:17146/143‑1 (MQ=255)
                                                                     ggCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGaa                                                               >  2:1077375/1‑143 (MQ=255)
                                                                      gCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATCGGCTTAAGGTAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGaat                                                              >  2:797188/1‑143 (MQ=255)
                                                                      gCATAGAGCCTCCTGTAGGGGTTTTATTAACAACGGCTTATTCTCATTATTTTGTGATGAGAGGCAGCGCGTGCCGCTCAAATATTAAAGTGAATAGGGTATAGTGAACGCGGTGACAACTAATAGGAATCACGGTTGTGagt                                                              >  2:616175/1‑141 (MQ=255)
                                                                        aTAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGaataa                                                            >  2:425655/1‑143 (MQ=255)
                                                                           gagCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGc                                                         >  2:584265/1‑143 (MQ=255)
                                                                           gagCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTGTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGc                                                         >  2:814244/1‑143 (MQ=255)
                                                                                   tGTAGTTTTTTTATTTACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAAGAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCg                                                 <  1:152088/143‑1 (MQ=255)
                                                                                    gTAGGGTTTTTAGTAACCACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCgg                                                <  1:161257/143‑1 (MQ=255)
                                                                                       gggTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTATGACGGTGCACTGAACGCGGTTAAAAATAACAGGCATCACGGTGGTGAATGCTGCCCCGTACGGCCa                                             >  2:918210/1‑143 (MQ=255)
                                                                                           ttttATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGcc                                         >  1:624063/1‑143 (MQ=255)
                                                                                            tttAGTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTGTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCg                                        <  1:937268/143‑1 (MQ=255)
                                                                                                  aacaacGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCCAGCCGGATTaa                                  >  2:92714/1‑143 (MQ=255)
                                                                                                   acaacGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCACGCCGGATTAAt                                 >  1:1040842/1‑143 (MQ=255)
                                                                                                     accGGCTTTTTTTAATTATGTGGAGATGAGCGGCATCGCGTGCCGCTCAAAAATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCg                               <  2:161264/141‑1 (MQ=255)
                                                                                                     aacGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTACAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCg                               >  1:605248/1‑143 (MQ=255)
                                                                                                      acGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGc                              <  1:444519/143‑1 (MQ=255)
                                                                                                       cGGCTTATTCTAATTATTTTGTGAAGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGcc                             <  1:278011/143‑1 (MQ=255)
                                                                                                         gCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGt                           <  1:558034/143‑1 (MQ=255)
                                                                                                           ttATGCTAATTATTTTGTGATGAGCGGCAGCGCGTCCCGCTGAAATATTACAGTTAATAGGGTTTAGAGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCa                         <  2:902621/143‑1 (MQ=255)
                                                                                                            tATTCCAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCaa                        <  1:1068986/143‑1 (MQ=255)
                                                                                                              ttCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAAtt                      >  2:422087/1‑143 (MQ=255)
                                                                                                               tCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTCATCGCCGCCAATTc                     >  2:749873/1‑143 (MQ=255)
                                                                                                               tCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTc                     >  1:657450/1‑143 (MQ=255)
                                                                                                               tCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATACATCGCCGCCAATTc                     >  2:854952/1‑143 (MQ=255)
                                                                                                                   attattTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGCCCAAGCCGGATTAATCGCCGTCAATTCAgcg                 >  1:764971/1‑143 (MQ=255)
                                                                                                                     tattTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGcgcg               >  2:317606/1‑143 (MQ=255)
                                                                                                                      actttGTGATGACCTGCAGCGCGTGCCACTCAAATATGAAAGTTAATAGTGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGa              <  2:939418/141‑1 (MQ=255)
                                                                                                                           gtgATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGGTGAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGa         <  1:359671/143‑1 (MQ=255)
                                                                                                                           gtgATGAGCGGCAGCGCGTGCCGCTCAAACATTACAGTTAATAGGGGTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGa         <  1:254196/143‑1 (MQ=255)
                                                                                                                            tgatgaGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCTGTCAATTCAGCGCGAGCAgca        >  2:1005330/1‑141 (MQ=255)
                                                                                                                             gatgaGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGGGCGCGTAGACt       >  2:857255/1‑143 (MQ=255)
                                                                                                                              ttgaGCGGCACGCGCTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAGATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCAGGCTGAATCGCCGTCAATTCAGCGCGCGCAGACTc      >  1:450058/2‑143 (MQ=255)
                                                                                                                              atgaGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTc      <  1:353083/143‑1 (MQ=255)
                                                                                                                               tgaGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGCGCAGACTCa     >  1:388625/1‑143 (MQ=255)
                                                                                                                               tgaGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATCCCGCGCGAGCAGCCTCa     >  1:145701/1‑143 (MQ=255)
                                                                                                                               tgaGCGGCAGCGCGTGACGAAAAAATAGTACAGTTAATAGGGTTTAGTGAACACGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCa     <  2:776309/143‑1 (MQ=255)
                                                                                                                               tgaGCGGCACGCGCTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGAGTAATCGACGTCAATTCAGCGCGAGCAGACCCa     >  2:886610/1‑143 (MQ=255)
                                                                                                                                gaGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTcgt    >  1:491786/1‑141 (MQ=255)
                                                                                                                                  gCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCatat  >  1:769013/1‑143 (MQ=255)
                                                                                                                                        |                                                                                                                                         
CTGGCAAATAGTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATAT  >  NC_000913/2679294‑2679566

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