Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,242,049	T→C	7.3%	intergenic (+47/‑15)	PP_t28 → / → PP_t29	tRNA‑Glu/tRNA‑Ala

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,242,049	0	T	C	7.3%	157.0 / 2.2	56	intergenic (+47/‑15)	PP_t28/PP_t29	tRNA‑Glu/tRNA‑Ala
Reads supporting (aligned to +/- strand):  ref base T (24/27);  new base C (2/2);  total (27/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 = 9.96e-01

GCGTTAGCGGATCGGCACGAAGGTTACGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCACGCAAGTGGATCTTCGCGAAGGTTACGTCCCCTTCGTCTAG  >  NC_002947/2241900‑2242196                                                                                                                                                      |                                                                                                                                                    gcgTTAGCGGATCGGCACGAAGGTTACGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTT‑AAAAGGCTCGTCGATTATtgat                                                                                                                                                     <  5:79820/149‑1 (MQ=255)    ttAGCGGATCGGCACGAAGGTTACGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATtgatga                                                                                                                                                   <  6:2861/149‑1 (MQ=255)          tatgtgCACGCAGGTTACTTGCCTTTAGTATAGTTCCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCATCAATCATCGACGAGCCttt                                                                                                                                             <  7:346492/144‑1 (MQ=17)          gATCGGCACGAAGGTTACGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCttt                                                                                                                                             >  8:188610/1‑149 (MQ=255)           aTCGGCACGAAGGTTACGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCtttt                                                                                                                                            <  3:68975/149‑1 (MQ=255)                  cGAAGGTTACGTCCCCTTCGTCTAGTGTCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCgg                                                                                                                                     <  2:309011/149‑1 (MQ=255)                   gAAGGTTACGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCggg                                                                                                                                    <  1:364203/149‑1 (MQ=255)                     aGGTTACGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGc                                                                                                                                  <  7:125604/149‑1 (MQ=255)                        ttACGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCta                                                                                                                                >  5:152917/1‑148 (MQ=255)                         tACGTCCCCTTCGTTTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCATCAATAATCGACGAGCCtttttaa                                                                                                                                         >  2:15859/1‑135 (MQ=17)                         tACGTCCCCTTCGTTTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCATCAATAATCGACGAGCCtttttaa                                                                                                                                         <  1:15859/138‑4 (MQ=17)                            gTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCt                                                                                                                           <  7:111750/149‑1 (MQ=255)                               cccTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAg                                                                                                                        >  7:113981/1‑149 (MQ=255)                                 cTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCt                                                                                                                      >  2:86079/1‑149 (MQ=255)                                    cGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTgg                                                                                                                    <  6:152917/148‑1 (MQ=255)                                            ggCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACACCACTTTGCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAgcgc                                                                                                            <  7:262379/148‑1 (MQ=255)                                             gCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTgg                                                                                                                    <  4:89908/139‑1 (MQ=255)                                             gCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTgg                                                                                                                    >  3:89908/1‑139 (MQ=255)                                             gCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCt                                                                                                          >  2:181279/1‑149 (MQ=255)                                             gCCTAGGACACAGCCCTTTCACGGCGGTAACAGGCGTTCGAGTCCCCTAGGGGACGCCACTTTTCCTGCGTTTTGCATGGCTTAAAAAGGCTCGTCGATTATTGAAGAGACTTTTGTTTCGGGGCTATAGCTCACCTGGGAGAACGCCt                                                                                                          >  3:236954/1‑149 (MQ=25)                                                tAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTgca                                                                                                       >  1:290672/1‑149 (MQ=255)                                                  ggACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTgcat                                                                                                      <  8:113981/148‑1 (MQ=255)                                                  ggACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTgca                                                                                                       <  8:191808/147‑1 (MQ=255)                                                    acacCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTgcatggc                                                                                                   <  6:189563/149‑1 (MQ=255)                                                      acCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTgcatggcat                                                                                                 >  2:217086/1‑149 (MQ=255)                                                      acCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTgcatggcat                                                                                                 >  5:47909/1‑149 (MQ=255)                                                      acCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTgcatggcat                                                                                                 >  4:38563/1‑149 (MQ=255)                                                          cccTTTCACGGCGGTAACAGGGGTTCGAGGCCCCTAGGGGACGCCACTTTTCCTGCGTTTTGCAAGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCag                                                                                             <  7:173252/149‑1 (MQ=255)                                                               tCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCaggag                                                                                          <  2:114503/147‑1 (MQ=255)                                                                 aCGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGGCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCa                                                                                      <  1:366905/149‑1 (MQ=255)                                                                        tAACAGGGGTTCGAGTCCCCTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTc                                                                               <  1:181279/149‑1 (MQ=255)                                                                           cAGGGTTTCTAGTCCCCTATGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTAAAAAACGCTCGTCGATTATTGATGATCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGAt                                                                            <  4:236954/149‑1 (MQ=37)                                                                               ggTTCGAGTCCCCTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCg                                                                        <  7:188610/149‑1 (MQ=255)                                                                               ggTTCGAGTCCCCTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCg                                                                        <  3:221268/149‑1 (MQ=255)                                                                                 ttCGAGTCCCCTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCt                                                                      >  2:371426/1‑149 (MQ=255)                                                                                         cccTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCACAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCTCTTAGGTcca                                                              >  6:348637/1‑149 (MQ=255)                                                                                          ccTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTccac                                                             >  1:3695/1‑149 (MQ=255)                                                                                           cTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTccac                                                             >  3:48066/1‑148 (MQ=255)                                                                                           cTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTccac                                                             >  8:168442/1‑148 (MQ=255)                                                                                               gggACACCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAtt                                                         >  6:333677/1‑148 (MQ=255)                                                                                                    gCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGCTCCCGCTTAGCTCCACCATTTCTGcc                                                   <  1:217086/149‑1 (MQ=255)                                                                                                    gCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGc                                                    >  5:323236/1‑148 (MQ=255)                                                                                                      cACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAg                                                 <  6:47909/149‑1 (MQ=255)                                                                                                          tttCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGt                                              >  5:373750/1‑148 (MQ=255)                                                                                                               tGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCACg                                         >  3:112462/1‑148 (MQ=255)                                                                                                                        cAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCACGCAAGTGGAt                                <  1:371426/148‑1 (MQ=255)                                                                                                                            gaaacaaaaGGCTCATCAATAATCGACGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCACGCAAGTGGATCTTcg                           >  6:419575/6‑149 (MQ=21)                                                                                                                            gCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCACGCAAGTGGATCTTcg                           <  4:112462/149‑1 (MQ=255)                                                                                                                              ttAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCt                                                                      <  4:224902/104‑1 (MQ=255)                                                                                                                              ttAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCt                                                                      >  3:224902/1‑104 (MQ=255)                                                                                                                                 aaaaGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCACGCAAGTGGATCTTCGCGAAg                      <  6:323236/149‑1 (MQ=255)                                                                                                                                 aaaaGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCACGCAAGTGGATCTTCGCGAAg                      <  6:373750/149‑1 (MQ=255)                                                                                                                                 aaaaGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCACGCAAGTGGATCTTCGCGAAg                      <  1:86079/149‑1 (MQ=255)                                                                                                                                  aaaGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCACGCAAGTGGATCTTCGCGAAgg                     >  5:72489/1‑149 (MQ=255)                                                                                                                                    aGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCACGCAAGTGGATCTTCGCGAAGGtt                   <  3:38563/149‑1 (MQ=255)                                                                                                                                        tcgtcgATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCACGCAAGTGGATCTTCGCGAAGGTTACGt               >  1:62325/1‑149 (MQ=255)                                                                                                                                            cgATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCACGCAAGTGGATCTTCGCGAAGGTTACGTcccc           >  7:71625/1‑149 (MQ=255)                                                                                                                                              attactGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCACGCAGGTGGATCTTCGCGAAGGGGACGTAACCtt         >  5:364921/1‑149 (MQ=21)                                                                                                                                                    gatgaGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCACGCAAGTGGATCTTCGCGAAGGGGACGTAACCTTCgcga    <  6:364921/148‑4 (MQ=21)                                                                                                                                                     atgaGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCACGCAAGTGGATCTTCGCGAAGGTTACGTCCCCTTCGTCTAg  <  6:72489/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                    GCGTTAGCGGATCGGCACGAAGGTTACGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCACTTTTCCTGCGTTTTGCAGGGCTTAAAAAGGCTCGTCGATTATTGATGAGCCTTTTGTTTCGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCATTTCTGCCAGGGTTCACGCAAGTGGATCTTCGCGAAGGTTACGTCCCCTTCGTCTAG  >  NC_002947/2241900‑2242196

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: ‑