Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,923,417	A→G	7.4%	P88P (CCT→CCC)	rsmG ←	16S rRNA m(7)G527 methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,923,417	0	A	G	7.4%	180.1 / 4.9	70	P88P (CCT→CCC)	rsmG	16S rRNA m(7)G527 methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (39/23);  new base G (2/3);  total (42/28)
Fisher's exact test for biased strand distribution p-value = 3.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.24e-02

ATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTC‑CAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAG  >  NC_000913/3923279‑3923559                                                                                                                                            |                                                                                                                                               aTCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTC‑CAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGAc                                                                                                                                            <  1:770810/143‑1 (MQ=255) aTCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATACTCTC‑CAGTTTAAGCTCCTGTTGCATGGGAAGCCCTTGTTGCACGTGTTTACCAAGGCTATCCAGCAGAGTGGAATGGGCGTCCGGAc                                                                                                                                            >  2:406577/1‑143 (MQ=255)         gcggcTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTC‑CAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATa                                                                                                                                    <  2:798709/143‑1 (MQ=255)         gcggcTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATCTTCGCACCGCCA‑CTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATa                                                                                                                                    <  1:668777/143‑1 (MQ=255)          ctgcgctGAACGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTC‑CAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATag                                                                                                                                   <  2:830491/138‑1 (MQ=255)           ggcTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTC‑CAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGGAATGGGCTTCAGGACGCACGATcga                                                                                                                                  >  1:503569/1‑143 (MQ=255)             ctctGACGGCACCTCTTCTACCCTGCTCTGGACTGGTTCAACATGCTC‑CAGTGTAAGCTCATGTTGCACCTGCCGAAGGAAACGCACGCGTTGACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATagaga                                                                                                                                <  2:333456/143‑1 (MQ=255)             ccctGAAGGAAACTCTTCCCCCCTGCTCTGTACTGGTTCAATATTCCC‑CAGTTTAAGCTCATGTTGCACCTGGCGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATagaga                                                                                                                                <  1:506058/141‑1 (MQ=255)              tctGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTC‑CAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAATGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATagagag                                                                                                                               >  2:652709/1‑143 (MQ=255)                  aaGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTC‑CAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGa                                                                                                                           >  1:215022/1‑143 (MQ=255)                    ggAAACTCTTCTACCCTGCTCTGTACTGGTTCAATCTTCTC‑CAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCCACCGAGTGAAATGGCCTTCAGGCCGCACGATAGCGAGTGGAAt                                                                                                                         >  1:917635/1‑143 (MQ=255)                        cctctCCTACCCTGCTCTGTACTGGTTCAATATTCTC‑CCGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCGGGACGCACGATAGAGAGTGGAATGCCt                                                                                                                     <  1:226050/142‑1 (MQ=255)                        aCTCTTCTACCCTGCTCTGTACTGGTTCAATAATCTC‑CAGTTTAAGCTCATCTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGCACGCACGATCGAGAGTGGAATGCCt                                                                                                                     >  2:482170/1‑143 (MQ=255)                         ctctTCTACCCTGCTCTGTACTGGTTCAATATTCTC‑CAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATCGAGAGAGGAATACCTg                                                                                                                    >  2:340794/1‑143 (MQ=255)                            ttatACCCTGCACTGTAGTGGTTCAATATTCTC‑CAGTTTAAGCGCATGTTGGACCTGACGAAGGAAGCGCCGGCGATTACCAAGGCTATACAACAGAGTGAAAGGGGCTCCAGGACGCACGATGGAGAGAGGGATGACTGTCa                                                                                                                 >  2:254007/4‑143 (MQ=255)                              ctACCCTGCTCTGTACTGGTTCAATATTCTC‑CAGTTTAAGCTCATGTTGCACCTGACGAAAGAAACGCATGCGATTACCAAGGCTCTCCAACAGAGTAAAATGAGCTTGAGGACGCACGCTAGAGGGTGGAATGCCTGGCAgg                                                                                                               >  2:298137/1‑142 (MQ=255)                               tACCCTGCTCTGTACTGGTTCAATATTCTC‑CAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTGCCAAGGCTATCCAACAGAGTGAAATGGGCTTCGGGACGCACGATAGAGAGTGGAATGCCTGGCAGTc                                                                                                              >  1:583954/1‑143 (MQ=255)                               tACCCTGCTCTGTACTGGTTCAATATTCTC‑CAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCGGGACGCACGATAGAGAGTGGAATGCCTGGCAGTc                                                                                                              >  1:141573/1‑143 (MQ=255)                                aCCCTGCTCTGTACTGGTTCAATATTCTC‑CAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTcc                                                                                                             >  1:233540/1‑143 (MQ=255)                                 cccTGCTCTGTACTGGTTCAATATTCTC‑CAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTccc                                                                                                            >  2:355952/1‑142 (MQ=255)                                  ccTGCTCTGTACTGGTTCAATATTCTC‑CAGTTTAAGGTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTg                                                                                                           <  1:677486/143‑1 (MQ=255)                                     gCTCTGTACTGGTTCAATATTCTC‑CAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTGCCAAGGCTATCCAACAGAGTGAAATGGGCTTCGGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTc                                                                                                        <  2:583954/143‑1 (MQ=255)                                     gCTCTGTACTGGTTCAATATTCTC‑CAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCGGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTc                                                                                                        <  2:141573/143‑1 (MQ=255)                                        ctGTACTGGTTCAATATTCTC‑CAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAAAGCGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGGCCTGGGCCgg                                                                                                     >  2:48167/1‑143 (MQ=255)                                         tGTACTGGTTCAATATTCTC‑CAGTTTACGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCACCAGAGTGCAATGGGCTTCAGGACGCACGGTAGAGAGTGGACTGCCCGGCAGTCCTGGACCGGt                                                                                                    >  1:476800/1‑143 (MQ=255)                                         tGTACTGGTTCAATATTCTC‑CAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGt                                                                                                    <  1:451290/143‑1 (MQ=255)                                           tACTGGTTCAATATTCTC‑CAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAACGGGCTTCAGGACGCACGATAGAGAGTGGAATGCGTGGCAGTCATTGTcccgcgc                                                                                                  >  2:85022/1‑138 (MQ=255)                                            aCTGGTTCAATATTCTC‑CAGTTTAAGCGCATGTGGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGcc                                                                                                 <  1:601650/143‑1 (MQ=255)                                                    aaTATTCTC‑CAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGa                                                                                         <  1:407061/143‑1 (MQ=255)                                                           tc‑CAGTTTAAGCTCATGTTGCACCTGACGACGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGCGTGGTATGCCTGGCAGGCCTGATCCGGTGCCGACATCGATAAACCg                                                                                  >  2:317213/1‑143 (MQ=255)                                                           tc‑CAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAAcgg                                                                                  >  1:814904/1‑141 (MQ=255)                                                           tc‑CAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCg                                                                                  >  2:22353/1‑143 (MQ=255)                                                                 tttAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACCACGATAAACCGTTcgc                                                                             >  2:282189/1‑141 (MQ=255)                                                                       cTCATGTTGCACCTGACGAAGGAAACGCACGCGTTGACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCa                                                                       <  1:645363/143‑1 (MQ=255)                                                                        tCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAg                                                                      >  1:874337/1‑143 (MQ=255)                                                                        tCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCGTGCAg                                                                      >  1:184476/1‑143 (MQ=255)                                                                         cATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGGATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTGCACCTTGCAGa                                                                     >  1:641707/1‑143 (MQ=255)                                                                         cATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGAGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTACCGACATCGATAAACCGTTCACCCCGCAGa                                                                     >  2:552968/1‑143 (MQ=255)                                                                             ttGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCATTCACCTTGCAGATACg                                                                 >  2:143365/1‑143 (MQ=255)                                                                              tccACCTGACGAAGGAAACGCACGCGTTGACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACgg                                                                <  1:536462/141‑1 (MQ=255)                                                                               gCACGTGAGGGAGGAAACGCACGCGTTGAACAAGGCTACCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCGGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGt                                                               <  1:106174/143‑1 (MQ=255)                                                                               gCACCTGGCGAAGGAAACGCACGCGTTGACCAAGGGTATCCAACAGAGTGAAATGGGCTGCAGGACGCACGAGAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGt                                                               <  1:181517/143‑1 (MQ=255)                                                                                   cTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCTACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGcca                                                           >  2:499525/1‑143 (MQ=255)                                                                                    ggaggaaggaagCGCCTGCGGTGACCAAGGCTAGCCAACAGAGTGAAGTGAGCTTCAGGACGCACGTGATACAGCGGAACGCGTGGCAGACCTGGGCCGGTGCCGACATCGATAAACCGTCCACCTTGCAGATACGGTGccac                                                          <  2:720018/139‑1 (MQ=255)                                                                                      aCGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGccacca                                                        >  2:821330/1‑143 (MQ=255)                                                                                        gacggcaaCGCACGCGTGGACCAAGGGTATCCCACAGAGTGGAAGGGGCTTGAGGACGCACGATTGAGAGTGGAATGCCTGGCAGTCCTGGTGCGGTGCCTACATCGATAAACCGTTCACCTTGCTGATACGGTGCCACcaca                                                      <  2:740856/137‑1 (MQ=255)                                                                                        gAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGAGGAATGCCTGGCAGTCATGGTCCGGTGCCGACATCGATAAACCGATCACCTTGCAGATCCGGTGCCACcaca                                                      >  2:275304/1‑143 (MQ=255)                                                                                         aaGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACaa                                                     <  1:539614/143‑1 (MQ=255)                                                                                           ggAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATg                                                   >  2:573318/1‑143 (MQ=255)                                                                                            gAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCCCCACAATGc                                                  >  2:253733/1‑143 (MQ=255)                                                                                              aaCGCACGCGTTTACCAAGGCTATCCAACAGAGAGAAATAGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTGCGGTGCCGACAACGATAAACCGTTCACCTTGCAGATACGGTGCGACCACAATGATa                                                >  2:917470/1‑143 (MQ=255)                                                                                                cgcacgcGGTTAGCAAGGCTATCCAACAGAGGGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATc                                              <  1:140263/143‑1 (MQ=255)                                                                                                        tttACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACGATGCTGTCCAGAatat                                      >  1:392591/1‑143 (MQ=255)                                                                                                             cAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTCGCAGATACGGCGCCACCACAATGCTATCGAGACTATGGCGt                                 >  2:424294/1‑143 (MQ=255)                                                                                                                  cTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAg                            >  1:791109/1‑143 (MQ=255)                                                                                                                  cTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCAACTTGTAGAGACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAg                            <  1:388514/143‑1 (MQ=255)                                                                                                                         aCAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTCTCGCGAATATGGCGTACCAGCACCTCa                     >  2:562292/1‑143 (MQ=255)                                                                                                                           agagCGAAATGGGCTTCATGACGCCCGATCGAGAGTGGAATGCCTCCCCGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCAtt                   <  1:95062/143‑1 (MQ=255)                                                                                                                            gagTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTa                  <  1:70663/143‑1 (MQ=255)                                                                                                                             ggggAAATGGGCGTCTGGAGGCGCGATAGAGCGGGGAATGCCGGGCAGTGCTGGTCCGGGGCCGACATCGATAGACCGTTCGCCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAg                 <  2:341563/140‑1 (MQ=255)                                                                                                                                 aaaTGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATAGCGTACCAGCATCTCATTAggctt             >  2:120810/1‑140 (MQ=255)                                                                                                                                  aaTGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATcg            <  2:374550/143‑1 (MQ=255)                                                                                                                                    tGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAATCCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATcgcg          <  2:107220/143‑1 (MQ=255)                                                                                                                                     gggCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAGTATGGCGTACCACCATCTCATTAGGCTCGCgg         >  2:701744/1‑143 (MQ=255)                                                                                                                                       tcTTCCGGACCCCCGATAGAGAGTCGAAGGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTGGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGAc       <  1:355952/142‑1 (MQ=255)                                                                                                                                       gCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGGTAGCCGTGCAGATACGGTGCCACCACCATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGag       >  2:553146/1‑142 (MQ=255)                                                                                                                                        cTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATCCGGTGCCACCACAATGCTATCGAGCATATCGCGTACCAGCATCCCATTAGGATCGCGGAcc      >  2:385722/1‑143 (MQ=255)                                                                                                                                        cTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTAGCGAGGATATGGCGTACCAGCATCTCATTAGGATCGCGGAcc      >  1:170409/1‑143 (MQ=255)                                                                                                                                          tCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCGGGTGCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGCGCCACCAAAATGATATCAAAAATATGGCGTACCAGCATATAATCAGGAGCGCGGACCGa    >  2:7387/1‑143 (MQ=255)                                                                                                                                           cAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTACGATCGCGGACCGaa   >  1:374734/1‑143 (MQ=255)                                                                                                                                            aGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGCATATAGCATACCAGCAGCTCGTTAGGATCGCGGACCGAAg  >  1:619860/1‑143 (MQ=255)                                                                                                                                            |                                                                                                                                               ATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTC‑CAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAG  >  NC_000913/3923279‑3923559

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