Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,224,873	C→T	5.8%	E601E (GAG→GAA) ‡	PP_1963 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,224,873	0	C	T	5.8%	202.3 / 2.6	69	E601E (GAG→GAA) ‡	PP_1963	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (32/33);  new base T (2/2);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.61e-01

TAATAGCGGGGCGTTCCGCCGTCTATATCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGAT  >  NC_002947/2224726‑2225021                                                                                                                                                    |                                                                                                                                                     tAATAGCGGGGCGTTCCGCCGTCTATATCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTc                                                                                                                                                    >  1:239965/1‑150 (MQ=255) tAATAGCGGGGCGTTCCGCCGTCTATATCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTc                                                                                                                                                    <  1:138816/150‑1 (MQ=255)   aTAGCGGGGCGTTCCGCCGTCTATATCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCaa                                                                                                                                                  <  2:43151/150‑1 (MQ=255)         gggCGTTCCGCCGTCTATATCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAg                                                                                                                                            >  2:654873/1‑150 (MQ=255)           gCGTTCCGCCGTCTATATCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATATTGAGTATGAGCt                                                                                                                                          <  1:238581/150‑1 (MQ=255)             gTTCCGCCGTCTATATCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAg                                                                                                                                        <  2:239965/150‑1 (MQ=255)                         atatCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCaa                                                                                                                            >  2:1425644/1‑150 (MQ=255)                          tatCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCaaa                                                                                                                           <  1:1220076/150‑1 (MQ=255)                                     cTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCaa                                                                                                                                                  >  1:98176/1‑116 (MQ=255)                                     cTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCaa                                                                                                                                                  <  2:98176/116‑1 (MQ=255)                                      tGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAg                                                                                                               >  2:1176026/1‑150 (MQ=255)                                        cTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTc                                                                                                             <  1:654873/150‑1 (MQ=255)                                              ggTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGCGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTAGTa                                                                                                       >  1:166057/1‑150 (MQ=255)                                                  tCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATa                                                                                                                         <  1:125810/128‑1 (MQ=255)                                                  tCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATa                                                                                                                         >  2:125810/1‑128 (MQ=255)                                                       atcatGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCaaa                                                                                                                           <  2:378607/121‑1 (MQ=255)                                                       atcatGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCaaa                                                                                                                           >  1:378607/1‑121 (MQ=255)                                                       atcatGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAgg                                                                                                                       <  2:1274933/125‑1 (MQ=255)                                                       atcatGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAgg                                                                                                                       >  1:1274933/1‑125 (MQ=255)                                                        tcatGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAAcc                                                                                             >  1:658149/1‑150 (MQ=255)                                                           tGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCCgg                                                                                          <  1:1425644/150‑1 (MQ=255)                                                             gatctCGTTGATGAGAAACTCCAAAATTTTACCTTCCATTCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGa                                                                                        <  1:1007676/146‑1 (MQ=255)                                                             aTTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGa                                                                                        <  2:259457/150‑1 (MQ=255)                                                                   gTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAATTAGGGTAGGGAGGGTTCTTAAATGCGAACCTGGGagagata                                                                                  >  2:1007676/1‑145 (MQ=255)                                                                   gTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGcc                                                                                  >  1:1226499/1‑150 (MQ=255)                                                                       atgaGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGaa                                                                                               >  2:687197/1‑133 (MQ=255)                                                                       atgaGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGaa                                                                                               <  1:687197/133‑1 (MQ=255)                                                                         gagaAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCa                                                                            >  2:536114/1‑150 (MQ=255)                                                                           gaAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCAc                                                                           <  2:549857/149‑1 (MQ=255)                                                                           gaAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCAc                                                                           >  1:549857/1‑149 (MQ=255)                                                                                tCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGa                                                                                        <  2:1089537/131‑1 (MQ=255)                                                                                tCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGa                                                                                        >  1:1089537/1‑131 (MQ=255)                                                                                  cAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTcc                                                                                                                                 >  1:1077499/1‑88 (MQ=255)                                                                                  cAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTcc                                                                                                                                 <  2:1077499/88‑1 (MQ=255)                                                                                  cAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTAt                                                                   >  2:322293/1‑150 (MQ=255)                                                                                  cAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTAt                                                                   >  1:1014807/1‑150 (MQ=255)                                                                                   aaaaTTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCaa                                                                                                                                                  >  2:882875/1‑70 (MQ=255)                                                                                   aaaaTTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCaa                                                                                                                                                  <  1:882875/70‑1 (MQ=255)                                                                                   aaaaTTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTc                                                                                                  >  1:1276133/1‑118 (MQ=255)                                                                                   aaaaTTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTc                                                                                                  <  2:1276133/118‑1 (MQ=255)                                                                                   aaaaTTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATg                                                                  >  2:980543/1‑150 (MQ=255)                                                                                   aaaaTTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACATTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGCGCAGCCGATGCACTTGTCTATg                                                                  >  1:25110/1‑150 (MQ=255)                                                                                              ttCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCt                                                                <  2:1095562/141‑1 (MQ=255)                                                                                              ttCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATAAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCt                                                                >  1:1095562/1‑141 (MQ=255)                                                                                               tCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCAtttc                                                      <  1:536114/150‑1 (MQ=255)                                                                                               tCCCAGCGCAGGTCAGAAGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTc                                                                                                             <  2:1262067/95‑1 (MQ=255)                                                                                               tCCCAGCGCAGGTCAGAAGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTc                                                                                                             >  1:1262067/1‑95 (MQ=255)                                                                                                      gcAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGt                                                                                                                  >  1:557797/1‑83 (MQ=255)                                                                                                      gcAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGt                                                                                                                  <  2:557797/83‑1 (MQ=255)                                                                                                       cAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGAt                                                                               <  1:652189/117‑1 (MQ=255)                                                                                                       cAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGAt                                                                               >  2:652189/1‑117 (MQ=255)                                                                                                        aGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTa                                             <  2:1226499/150‑1 (MQ=255)                                                                                                        aGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTa                                             <  2:721107/150‑1 (MQ=255)                                                                                                        aGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTa                                             <  1:1176026/150‑1 (MQ=255)                                                                                                          gTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCAATTAgg                                           >  1:139895/1‑150 (MQ=255)                                                                                                                     gTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGa                                                                                                                                             >  1:881035/1‑41 (MQ=255)                                                                                                                     gTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGa                                                                                                                                             <  2:881035/41‑1 (MQ=255)                                                                                                                     gTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCt                                <  2:640568/150‑1 (MQ=255)                                                                                                                      tGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTc                               <  2:658149/150‑1 (MQ=255)                                                                                                                           cTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTcc                                                 <  2:371174/127‑1 (MQ=255)                                                                                                                           cTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTcc                                                 >  1:371174/1‑127 (MQ=255)                                                                                                                                 tttCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTc                                                                      >  2:457408/1‑100 (MQ=255)                                                                                                                                 tttCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTc                                                                      <  1:457408/100‑1 (MQ=255)                                                                                                                                          ccAGCTCATACTCAATATGAGCTAGACGCGCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAATCTGGGAGTAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGtt           <  2:166057/150‑1 (MQ=255)                                                                                                                                           taGCTCATATTGAGTATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTg          >  1:861606/2‑150 (MQ=255)                                                                                                                                           cAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTg          >  1:586833/1‑150 (MQ=255)                                                                                                                                            aGCTCATATTGAGTATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTc                                                                                                          <  1:107479/53‑1 (MQ=25)                                                                                                                                            aGCTCATATTGAGTATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTc                                                                                                          >  2:107479/1‑53 (MQ=25)                                                                                                                                                   aCTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGAt  <  1:980543/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     TAATAGCGGGGCGTTCCGCCGTCTATATCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGAT  >  NC_002947/2224726‑2225021

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