Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,505,538	G→T	8.7%	G217V (GGC→GTC)	petB →	cytochrome b

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,505,538	0	G	T	8.7%	214.1 / 5.1	72	G217V (GGC→GTC)	petB	cytochrome b
Reads supporting (aligned to +/- strand):  ref base G (33/32);  new base T (3/4);  total (36/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGTGGATTCGTGGTGACTACCTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGAT‑TCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGTTCTTCTTCCCTGAAATGGG  >  NC_002947/1505390‑1505685                                                                                                                                                      |                                                                                                                                                    cAGTGGATTCGTGGTGACTACCTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGAT‑TCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATgg                                                                                                                                                     >  3:238686/1‑149 (MQ=255)      gATTCGTGGTGACTACCTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGAT‑TCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTc                                                                                                                                                 >  1:146365/1‑148 (MQ=255)      gATTCGGTGTGACTAACTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTTCATGTGGTTGCCCTGCCAATCGGTATCTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGt                                                                                                                                                  <  1:120207/148‑1 (MQ=255)        ttCGTGGTGAGTACCTGATCTCGGGCCTGAGCCTGAACCGCTTCTACGCCCTGCATGTGGTTGCCCTGCCAATCGTGAG‑GGTCGGCCTGGTGGTATTGCACATCCTCGCCTGGAGCGAAGTGGGGTCGCGCAACCCGGATGGCGTCGa                                                                                                                                               >  8:239383/1‑148 (MQ=255)         tCGTGGTGACTACCTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGAT‑TCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACa                                                                                                                                             >  3:377764/1‑149 (MQ=255)              gtgACTACCTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGAT‑TCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCaa                                                                                                                                         <  4:238686/148‑1 (MQ=255)                   tACCTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGAT‑TCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACGTCAagaaga                                                                                                                                    >  1:393914/1‑148 (MQ=255)                   tACCTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGAT‑TCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAagaaga                                                                                                                                    >  6:118133/1‑148 (MQ=255)                     ccTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGAT‑TCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGaaa                                                                                                                                  >  5:273262/1‑148 (MQ=255)                      cTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGAT‑TCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGaaaaa                                                                                                                                <  2:364814/149‑1 (MQ=255)                      cTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGAT‑TCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGaaa                                                                                                                                  >  3:315142/1‑147 (MQ=255)                       tGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGAT‑TCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGTCGACGCCATCagg                                                                                                                                        <  1:102157/140‑3 (MQ=255)                       tGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGAT‑TCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGTCGACGCCATCagg                                                                                                                                        >  2:102157/1‑138 (MQ=255)                         aTCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGAT‑TCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAgg                                                                                                                              <  4:315142/148‑1 (MQ=255)                          tctcGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGAT‑TCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGaaaaggac                                                                                                                             <  2:255015/148‑5 (MQ=255)                              gggCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGAT‑TCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGaaa                                                                                                                        <  1:384456/149‑1 (MQ=255)                                gCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGAT‑TCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAAc                                                                                                                      <  5:16608/149‑1 (MQ=255)                                gCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGAT‑TCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAAc                                                                                                                      >  1:393568/1‑149 (MQ=255)                                       ccTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGAT‑TCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTc                                                                                                               >  3:272175/1‑149 (MQ=255)                                             ccGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATAGTGAT‑TCTCGGCATGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGTCGTCGACATCAAGAAGAAAAAGGACGAAAACGTTATTAAGCTg                                                                                                          >  7:251379/1‑148 (MQ=255)                                              cGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGAT‑TCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGa                                                                                                        <  4:377764/149‑1 (MQ=255)                                                        cccTGCATGTGGTTGCCCTGCCAATCGTGAT‑TCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCg                                                                                              >  3:402864/1‑149 (MQ=255)                                                         ccTGCATGTGGTTGCCCTGCCAATCGTGAT‑TCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGTCGACGCCATCag                                                                                                                                         <  2:386170/105‑2 (MQ=255)                                                         ccTGCATGTGGTTGCCCTGCCAATCGTGAT‑TCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGTCGACGCCATCag                                                                                                                                         >  1:386170/1‑104 (MQ=255)                                                          cTGCATGTGGTTGCCCTGCCAATCGTGAT‑TCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGtt                                                                                            >  1:380861/1‑149 (MQ=255)                                                             gaTGTTGGTGCCGTGCCAATCTTTAT‑TCTGGGCCTGGTTGTATTGCACATCCTCGGCTTGCACGAAGTGGGTTCTAACAACCCTGATGTCGTCGACATCAAGAAGAAAAAGGAGGAAAACGGTTTTCCGCTGGATGGTATTCCGTTc                                                                                           <  8:251379/146‑1 (MQ=255)                                                              aTGTGGTTGCCCTGCCAATCGTGAT‑TCTCGGCCTGGTGGGATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCAc                                                                                        <  8:363310/149‑1 (MQ=255)                                                                   gTTGCCCTGCCAATCGTGAT‑TCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGt                                                                                    <  4:272175/148‑1 (MQ=255)                                                                    ttGCCCTGCCAATCGTGAT‑TCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGta                                                                                   <  1:336630/148‑1 (MQ=255)                                                                    ttGCCCTGCCAATCGTGAT‑TCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACAAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGtac                                                                                  >  2:106454/1‑149 (MQ=255)                                                                        ccTGCCAATCTTGCT‑TCTGGGCATGGTGGTATTGCAGATCCTCGCCTTGCTCGAATTGGGTTCGAACAACGCTGATGGCGTCGACATCAAGAAGAAAAAGGACTAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGtactac                                                                               <  6:273262/148‑1 (MQ=255)                                                                         cTGCCAATCGTGAT‑TCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTacac                                                                             <  1:295518/149‑1 (MQ=255)                                                                            ccAATCGTGAT‑TCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACGTCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACAcc                                                                            <  2:393914/147‑1 (MQ=255)                                                                             cAATCGTGAT‑TCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGtacta                                                                                >  4:246194/1‑142 (MQ=255)                                                                             cAATCGTGAT‑TCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGtacta                                                                                <  3:246194/142‑1 (MQ=255)                                                                                tCGTGAT‑TCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTaaa                                                                      <  2:146365/149‑1 (MQ=255)                                                                                   tGAT‑TCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATg                                                                                                      <  7:301684/114‑1 (MQ=255)                                                                                   tGAT‑TCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATg                                                                                                      >  8:301684/1‑114 (MQ=255)                                                                                    gAT‑TCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACa                                                                  >  4:160441/1‑149 (MQ=255)                                                                                     aT‑TCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACa                                                                  <  5:297151/148‑1 (MQ=255)                                                                                          tcGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTc                                                             >  3:57777/1‑149 (MQ=255)                                                                                              ccTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCg                                                         >  2:364628/1‑149 (MQ=255)                                                                                                 ggtggtATTGCACATCCTCGCCTTTCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCgtgg                                                      <  2:393568/149‑1 (MQ=255)                                                                                                  gtggtATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGtt                                                                         >  8:384822/1‑129 (MQ=255)                                                                                                  gtggtATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGtt                                                                         <  7:384822/129‑1 (MQ=255)                                                                                                  gtggtATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCgtggt                                                     >  2:353722/1‑149 (MQ=255)                                                                                                   tggtATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCgtggtg                                                    <  2:100767/149‑1 (MQ=255)                                                                                                       aTTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCTTGGTGgtgt                                                >  6:52469/1‑149 (MQ=255)                                                                                                               tCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACGCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTc                                        <  4:402864/149‑1 (MQ=255)                                                                                                                  tCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCgtggtg                                                    >  3:298541/1‑134 (MQ=255)                                                                                                                  tCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCgtggtg                                                    <  4:298541/134‑1 (MQ=255)                                                                                                                  tCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCgtg                                     >  1:384757/1‑149 (MQ=255)                                                                                                                      cTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCt                                 <  3:16605/149‑1 (MQ=255)                                                                                                                       ttGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTg                                >  6:162320/1‑149 (MQ=255)                                                                                                                         gCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTgcg                              >  6:347101/1‑149 (MQ=255)                                                                                                                         gCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTgc                               <  2:380861/148‑1 (MQ=255)                                                                                                                           aCGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTgcgtc                            >  7:352942/1‑147 (MQ=255)                                                                                                                           aCGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTgcg                              >  2:231925/1‑147 (MQ=255)                                                                                                                           aCGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGcc                            >  4:369817/1‑149 (MQ=255)                                                                                                                             gAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGcgtcgt                          <  8:352942/149‑1 (MQ=255)                                                                                                                             gAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCgt                          >  2:239753/1‑149 (MQ=255)                                                                                                                              aaGTGGGTTCGAACAACCCTGATGTCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCgtg                         <  3:192687/149‑1 (MQ=255)                                                                                                                                   ggTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGgtgt                     <  1:364628/148‑1 (MQ=255)                                                                                                                                    gTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGttc                   >  5:393698/1‑149 (MQ=255)                                                                                                                                    gTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGttc                   <  6:393698/149‑1 (MQ=255)                                                                                                                                    gTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGttc                   <  4:57777/149‑1 (MQ=255)                                                                                                                                         aacaacCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTc                                        >  5:356304/1‑123 (MQ=255)                                                                                                                                         aacaacCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTc                                        <  6:356304/123‑1 (MQ=255)                                                                                                                                             acCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGTTCTTCTTCCCt          <  5:347101/149‑1 (MQ=255)                                                                                                                                              cccTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGTTCTTCTTCCCt          >  6:316415/1‑148 (MQ=255)                                                                                                                                               ccTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGTTCTTCTTCCCTg         >  1:341739/1‑148 (MQ=255)                                                                                                                                                 tGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGttcttcttc             >  1:137064/1‑142 (MQ=255)                                                                                                                                                 tGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGttcttcttc             <  2:137064/142‑1 (MQ=255)                                                                                                                                                      gCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGTTCTTCTTCCCTGAAATggg  <  5:162320/148‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                    CAGTGGATTCGTGGTGACTACCTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGAT‑TCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGTTCTTCTTCCCTGAAATGGG  >  NC_002947/1505390‑1505685

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