Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,520,464	A→G	14.8%	intergenic (+55/+17)	cspD → / ← icd	DNA replication inhibitor/NADP(+)‑specific isocitrate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,520,464	0	A	G	14.8%	169.6 / 19.7	74	intergenic (+55/+17)	cspD/icd	DNA replication inhibitor/NADP(+)‑specific isocitrate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (35/28);  new base G (6/5);  total (41/33)
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.17e-01

GTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGTGCCCTTGATGATCAG  >  NC_002947/4520326‑4520609                                                                                                                                           |                                                                                                                                                   gTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCt                                                                                                                                          <  3:9335/149‑1 (MQ=255)  ttGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAACCGGTCGCTTCTTTTAAAGAAGCGGCCGAttttttt                                                                                                                                         <  7:149490/149‑4 (MQ=38)   tGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGGTTTCTGt                                                                                                                                        >  7:238630/1‑149 (MQ=255)     aaaTCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTg                                                                                                                                      >  3:80211/1‑149 (MQ=255)       aTCAAAGATGCCGCGGCCGCTGCTGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGcc                                                                                                                                    <  8:191579/149‑1 (MQ=255)        tCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCt                                                                                                                                   <  2:116533/149‑1 (MQ=255)           cctctGCCGCAGCGGCCGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTa                                                                                                                                 >  5:43939/5‑148 (MQ=255)               tGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAACCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACAtgt                                                                                                                            <  5:170196/149‑1 (MQ=255)                gCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGAc                                                                                                                                                   <  5:76566/125‑1 (MQ=255)                gCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGAc                                                                                                                                                   >  6:76566/1‑125 (MQ=255)                 ccGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATgtgt                                                                                                                          >  5:131151/1‑149 (MQ=255)                    cgGCCGCTGCGTCAGATGGTGCAAGAGCCTTTGGCGAGCCGCAAATCGAAACTGCAAAAGCCTGACTGCATTCGATGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTg                                                                                                                       <  1:311988/149‑1 (MQ=255)                               gCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATc                                                                                                            >  8:196067/1‑149 (MQ=255)                                     gCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCa                                                                                                        >  5:73863/1‑147 (MQ=255)                                     gCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCa                                                                                                        >  5:244078/1‑147 (MQ=255)                                       tGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGc                                                                                                    >  4:113309/1‑149 (MQ=255)                                        gCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGcc                                                                                                   >  3:243866/1‑149 (MQ=255)                                                     gCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTAc                                                                                       >  7:201425/1‑148 (MQ=255)                                                     gCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTAc                                                                                       <  3:113309/148‑1 (MQ=255)                                                       gAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAg                                                                                    >  1:108268/1‑149 (MQ=255)                                                           cGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTg                                                                                <  2:108268/149‑1 (MQ=255)                                                            gCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTgg                                                                               <  6:244078/149‑1 (MQ=255)                                                             cAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGTCGCTTCTTTTAAAGAAGCGGCCGAttttttt                                                                                                                                         <  2:202283/90‑4 (MQ=25)                                                             cAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGTCGCTTCTTTTAAAGAAGCGGCCGAttttttt                                                                                                                                         >  1:202283/1‑87 (MQ=25)                                                                 gCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGccc                                                                          >  8:102728/1‑149 (MQ=255)                                                                  cGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCt                                                                         <  4:80211/149‑1 (MQ=255)                                                                   gAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTc                                                                        >  1:195647/1‑149 (MQ=255)                                                                          cAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCaa                                                                 <  2:195647/148‑1 (MQ=255)                                                                           aaaaGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCAATc                                                               >  1:214585/1‑149 (MQ=255)                                                                               gCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCAATCGCTc                                                           >  1:340505/1‑149 (MQ=255)                                                                                  tGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCAATCGCTCGa                                                         >  1:4872/1‑148 (MQ=255)                                                                                       gcctgcGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCAATCGCTCGAAATCg                                                    >  5:59531/1‑148 (MQ=255)                                                                                        cctgcGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCAATCGCTCGAAATc                                                     <  3:227977/146‑1 (MQ=255)                                                                                        cctgcGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCAATCGCTCGAAATc                                                     >  4:227977/1‑146 (MQ=255)                                                                                        cctgcGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCAATCGCTCGAAATCg                                                    >  7:45726/1‑147 (MQ=255)                                                                                        cctgcGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCAATCGCTCGAAATCGTa                                                  <  4:243866/149‑1 (MQ=255)                                                                                        cctgcGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCAATCGCTCGAAATCGTa                                                  <  7:102728/149‑1 (MQ=255)                                                                                         ctgcGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCAATCGCTCGAAATCGTAg                                                 <  2:214585/149‑1 (MQ=255)                                                                                          tgcGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCAATCGCTCGAAATCGTAgg                                                >  2:157611/1‑149 (MQ=255)                                                                                                          aCAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCg                                 >  7:112949/1‑148 (MQ=255)                                                                                                            aaaaaaaTCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCAATCGCTCGAAATCGTa                                                  <  5:300213/129‑1 (MQ=255)                                                                                                            aaaaaaaTCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCAATCGCTCGAAATCGTa                                                  >  6:300213/1‑129 (MQ=255)                                                                                                            aaaaaaaTCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCaa                              >  2:275881/1‑149 (MQ=255)                                                                                                            aaaaaaaTCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCaa                              <  2:4872/149‑1 (MQ=255)                                                                                                               aaaaTCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTACTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATgg                           >  1:229947/1‑149 (MQ=255)                                                                                                                 aaTCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGgc                                                                              >  6:163806/1‑97 (MQ=255)                                                                                                                 aaTCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGgc                                                                              <  5:163806/97‑1 (MQ=255)                                                                                                                 aaTCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGgcg                         >  7:244701/1‑149 (MQ=255)                                                                                                                    aggtcGCTTGTTTTAAAGAAGCGGCCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATGGCCAAAGCGGGAGCTGCTTACCAGCTTGGCGCCCTC‑AATCAAGCGCTCGAAATCGTAGGTCACGCGCTTGGCCGCAATGGCGCCa                      >  8:320342/5‑149 (MQ=255)                                                                                                                     ggtcGCTTCTTTTAAAGAAGCGGCCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCAt                     >  7:276040/4‑149 (MQ=255)                                                                                                                     ggtcGCTTCTTTTAAAGAAGCGGCCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCAt                     >  5:160390/4‑149 (MQ=255)                                                                                                                     ggCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCAt                     <  6:43939/149‑1 (MQ=255)                                                                                                                      gtcGCTTCTTTTAAAGAAGCGGCCGGTTTTCTGTTGCCTTACATGTGTTTGATCTTTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCAATCGCTCGAAATCGTa                                                  <  1:340478/117‑1 (MQ=255)                                                                                                                      gtcGCTTCTTTTAAAGAAGCGGCCGGTTTTCTGTTGCCTTACATGTGTTTGATCTTTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCAATCGCTCGAAATCGTa                                                  >  2:340478/3‑119 (MQ=255)                                                                                                                       tcGCTTCTTTTAAAGAAGCGGCCGATTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTg                   >  4:75471/2‑149 (MQ=255)                                                                                                                       ccGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTg                   <  8:244701/149‑1 (MQ=255)                                                                                                                       ccGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTg                   <  2:340505/149‑1 (MQ=255)                                                                                                                             cttTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCAATCGCTCGAAATCGTAGGTCACg                                           >  6:49293/1‑119 (MQ=255)                                                                                                                             cttTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCAATCGCTCGAAATCGTAGGTCACg                                           <  5:49293/119‑1 (MQ=255)                                                                                                                             cttTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGTGcc              <  8:45726/148‑1 (MQ=255)                                                                                                                              tttAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGTGCCCt            >  1:268168/1‑149 (MQ=255)                                                                                                                               ttAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGTGCCCtt           >  4:307082/1‑149 (MQ=255)                                                                                                                               ttAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGTGCCCtt           >  5:304773/1‑149 (MQ=255)                                                                                                                                tAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGTGCCCtt           <  7:196067/148‑1 (MQ=255)                                                                                                                                  aaaGAAGCGGCCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGTGCCCTtga         <  8:276040/148‑1 (MQ=255)                                                                                                                                  aaaGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCTCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGTGCCCTtga         <  1:275881/149‑1 (MQ=255)                                                                                                                                  aaaGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGTGCCCTtgat        >  5:146292/1‑149 (MQ=255)                                                                                                                                  aaaGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGTGCCCTtga         <  6:304773/148‑1 (MQ=255)                                                                                                                                   aagaagCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGTGCCCTtgat        >  5:135201/1‑148 (MQ=255)                                                                                                                                     gaagCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGTGCCCTTgatgat     >  5:6069/1‑149 (MQ=255)                                                                                                                                      aagCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGTGCCCTTGATGATc    <  6:131151/149‑1 (MQ=255)                                                                                                                                        gCGGCCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGTGCCCTTGATGATCAg  <  6:160390/149‑1 (MQ=255)                                                                                                                                        gCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTg                   <  4:136207/132‑1 (MQ=255)                                                                                                                                        gCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTg                   >  3:136207/1‑132 (MQ=255)                                                                                                                                           |                                                                                                                                                   GTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTC‑CATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGTGCCCTTGATGATCAG  >  NC_002947/4520326‑4520609

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