Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,982,841	A→C	9.1%	*120C (TGA→TGC)	PP_3510 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,982,841	0	A	C	9.1%	201.4 / 10.0	66	*120C (TGA→TGC)	PP_3510	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (29/31);  new base C (2/4);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

TATGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCAT  >  NC_002947/3982700‑3982988                                                                                                                                              |                                                                                                                                                    tATGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGTGGACTGGGGAAACCGCCTGCAGGGTCAATGCCCCTGca                                                                                                                                              <  4:315339/149‑2 (MQ=255)   tGCTGGCACTGACCTGCCAGGCATTGGACCTCGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCAGGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGTGTCAATGCCCCTGcag                                                                                                                                             <  3:57849/148‑3 (MQ=255)     cTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGTCAATGCCCCTGcaggc                                                                                                                                           <  2:114457/148‑5 (MQ=255)       ggCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGa                                                                                                                                        <  3:161570/149‑1 (MQ=255)        gCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGa                                                                                                                                        <  2:63910/148‑1 (MQ=255)        gCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGa                                                                                                                                        <  2:91194/148‑1 (MQ=255)         cACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCCGCAGGGTCAATGCCCCTGcag                                                                                                                                             >  6:381807/1‑140 (MQ=255)         cACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCCGCAGGGTCAATGCCCCTGcag                                                                                                                                             <  5:381807/142‑3 (MQ=255)          aCTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGa                                                                                                                                     <  1:231886/149‑1 (MQ=255)              aCCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATa                                                                                                                                   <  6:361876/147‑1 (MQ=255)                cTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTAt                                                                                                                               <  1:204373/149‑1 (MQ=255)                cTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTAt                                                                                                                               <  3:364662/149‑1 (MQ=255)                         aTTGGACCACGCTGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCa                                                                                                                      >  1:327034/1‑149 (MQ=255)                         aTTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACCCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCa                                                                                                                      >  5:294901/1‑149 (MQ=255)                                cACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTAt                                                                                                                >  5:335034/1‑148 (MQ=255)                                    ccGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAAtt                                                                                                           <  1:3266/149‑1 (MQ=255)                                         aTACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGa                                                                                                      <  1:41881/149‑1 (MQ=255)                                             aGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGa                                                                                                                                        <  6:281268/111‑1 (MQ=255)                                             aGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGa                                                                                                                                        >  5:281268/1‑111 (MQ=255)                                                   aCCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCTGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTc                                                                                            <  6:294901/149‑1 (MQ=255)                                                    ccATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTcc                                                                                           <  2:327034/149‑1 (MQ=255)                                                          cGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATcc                                                                                     <  7:99305/149‑1 (MQ=255)                                                           gCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCg                                                                                    >  4:358729/1‑149 (MQ=255)                                                                gCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAg                                                                               >  2:308902/1‑149 (MQ=255)                                                                    gTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGt                                                                                                        <  2:135199/120‑1 (MQ=255)                                                                            tGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAgg                                                                   >  6:86953/1‑149 (MQ=255)                                                                            tGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAg                                                                    <  3:358729/148‑1 (MQ=255)                                                                              aCAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTcc                                                                                           <  5:63413/123‑1 (MQ=255)                                                                              aCAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTcc                                                                                           >  6:63413/1‑123 (MQ=255)                                                                                aGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCaa                                                                                >  5:326058/1‑132 (MQ=255)                                                                                aGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCaa                                                                                <  6:326058/132‑1 (MQ=255)                                                                                  cGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCAt                                                             >  5:196066/1‑149 (MQ=255)                                                                                       agGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTgc                                                        <  6:335034/149‑1 (MQ=255)                                                                                           cTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATGGTATCGCTTCGCACACTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCa                                                     >  8:30165/1‑148 (MQ=255)                                                                                           cTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCa                                                     >  3:312234/1‑148 (MQ=255)                                                                                             gtgCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCata                                                   >  6:50711/1‑148 (MQ=255)                                                                                              tgCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCatat                                                  >  8:260129/1‑148 (MQ=255)                                                                                               gCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATaa                                                <  5:86953/149‑1 (MQ=255)                                                                                                 aCGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAAca                                              <  6:196066/149‑1 (MQ=255)                                                                                                          tGGCGTGGACTTGGGAAAACGCCTGCAGGTGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATaa                                                <  6:21793/138‑1 (MQ=255)                                                                                                          tGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATaa                                                >  5:21793/1‑138 (MQ=255)                                                                                                          tGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGa                                     <  1:147835/149‑1 (MQ=255)                                                                                                           ggCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAg                                    >  8:153030/1‑149 (MQ=255)                                                                                                                  aCTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCgg                             <  5:50711/149‑1 (MQ=255)                                                                                                                    tGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAg                                                                               <  6:71068/97‑1 (MQ=255)                                                                                                                    tGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAg                                                                               >  5:71068/1‑97 (MQ=255)                                                                                                                    tGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCggg                            >  1:116267/1‑148 (MQ=255)                                                                                                                      gggAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAAca                                              <  5:236782/128‑1 (MQ=255)                                                                                                                      gggAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAAca                                              >  6:236782/1‑128 (MQ=255)                                                                                                                      gggAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCAc                         >  7:54476/1‑149 (MQ=255)                                                                                                                      gggAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCAc                         >  4:138813/1‑149 (MQ=255)                                                                                                                        gAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATc                                                                                      <  5:254946/86‑1 (MQ=255)                                                                                                                        gAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATc                                                                                      >  6:254946/1‑86 (MQ=255)                                                                                                                          aaaCGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAAca                                              <  1:52797/124‑1 (MQ=255)                                                                                                                          aaaCGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAAca                                              >  2:52797/1‑124 (MQ=255)                                                                                                                           tcaaaccgCAGGGTCAATGCCCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGGATTAAGGTGGAATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGc                     >  5:166708/8‑148 (MQ=255)                                                                                                                            aCGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATGGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCaa                   <  7:30165/149‑1 (MQ=255)                                                                                                                              gCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGa                 >  1:14419/1‑149 (MQ=255)                                                                                                                              gCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGa                 >  4:224845/1‑149 (MQ=255)                                                                                                                                  gCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGGATGGTGCGCATATAACAGAAGCACGAGCTATCGGGCACTTGCATGACAAt             >  3:30536/1‑149 (MQ=255)                                                                                                                                   cAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACaa              <  2:116267/147‑1 (MQ=255)                                                                                                                                     ggggCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAg          <  7:260129/149‑1 (MQ=255)                                                                                                                                      gggCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAgg         <  4:312234/149‑1 (MQ=255)                                                                                                                                        gCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGgt       >  3:9496/1‑149 (MQ=255)                                                                                                                                          aTTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGgtgt     >  6:359907/1‑149 (MQ=255)                                                                                                                                          aTTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGgtg      >  4:200624/1‑148 (MQ=255)                                                                                                                                              aCCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCAt  >  8:136486/1‑148 (MQ=255)                                                                                                                                              |                                                                                                                                                    TATGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCAT  >  NC_002947/3982700‑3982988

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