Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,411,662	A→T	6.7%	N59I (AAT→ATT) ‡	PP_5551 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,411,662	0	A	T	6.7%	178.8 / 3.2	60	N59I (AAT→ATT) ‡	PP_5551	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (29/27);  new base T (2/2);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

ACTTGTCCTGTCTGTTGTAACGCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAGGTAAAGGACCTTTCGGCAAGATGTCTT  >  NC_002947/3411519‑3411810                                                                                                                                                |                                                                                                                                                     aCTTGTCCTGTCTGTTGTAACGCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGc                                                                                                                                                  >  2:39630/1‑148 (MQ=255)      tCCTGTCTGTTGTAACGCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCtt                                                                                                                                            >  8:286655/1‑149 (MQ=255)        ctgtctgtTGTAACGCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTc                                                                                                                                           <  6:15631/148‑1 (MQ=255)           tctgtTGTAACGCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGATTGACCTTCGCAAGCtcttc                                                                                                                                         <  4:348517/147‑5 (MQ=255)           tctgtTGTAACGCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCtt                                                                                                                                       <  1:76839/149‑1 (MQ=255)              gttgtAACGCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTgg                                                                                                                                     >  3:263154/1‑148 (MQ=255)                 gtAACGCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGaa                                                                                                                                   >  4:133903/1‑147 (MQ=255)                  tAACGCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGaaaa                                                                                                                                 <  7:286655/148‑1 (MQ=255)                  tAACGCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAg                                                                                                                                <  4:34982/149‑1 (MQ=255)                   aaCGCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGaa                                                                                                                                   >  8:92271/1‑145 (MQ=255)                      gCTGGTCTCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTTATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTa                                                                                                                            <  5:31072/149‑1 (MQ=255)                      gCTGGTCGCTCAGAATAAGCAAAGCACGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTa                                                                                                                            <  8:86002/149‑1 (MQ=255)                      gCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGATTGACCTTCGCAAGCtc                                                                                                                                            <  8:267200/133‑2 (MQ=255)                      gCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGATTGACCTTCGCAAGCtc                                                                                                                                            >  7:267200/1‑132 (MQ=255)                          gTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCTCAAGCTTCGCTTGGAAAAGAttcatc                                                                                                                         >  7:125462/1‑144 (MQ=255)                          gTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACt                                                                                                                          <  5:233167/147‑1 (MQ=255)                               tCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGAtt                                                                                                                   >  6:71657/1‑149 (MQ=255)                                cAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGAttt                                                                                                                  <  6:116530/149‑1 (MQ=255)                                 aGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTa                                                                                                                 <  5:133073/149‑1 (MQ=255)                                    aTAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGAc                                                                                                              <  5:263348/149‑1 (MQ=255)                                       agcaaagcaaGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCTCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGa                                                                                                            <  3:133903/148‑1 (MQ=255)                                        gcaaagcaaGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGc                                                                                                          >  7:256813/1‑149 (MQ=255)                                          aaagcaaGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCtt                                                                                                        >  1:328266/1‑149 (MQ=255)                                          aaagcaaGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCtt                                                                                                        >  4:331619/1‑149 (MQ=255)                                              caagcGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTgg                                                                                                    <  4:263154/149‑1 (MQ=255)                                                agcGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGa                                                                                                                     >  6:335987/1‑130 (MQ=255)                                                agcGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGa                                                                                                                     <  5:335987/130‑1 (MQ=255)                                                  cGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGt                                                                                                <  3:331619/149‑1 (MQ=255)                                                          ggCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTc                                                                                        <  7:92271/149‑1 (MQ=255)                                                           gCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCa                                                                                       <  2:328266/149‑1 (MQ=255)                                                           gCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCa                                                                                       <  5:78326/149‑1 (MQ=255)                                                           gCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTATTGGCCGTCACTCGTcc                                                                                       <  2:750/149‑2 (MQ=255)                                                                   ggCTAATTTAGCAAGTCAGGACGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTAGCTTGTAAAAGATTAATCAAGATTTAGACGAGCTTTTGGCTGTCACTCTACAGTATTgt                                                                                >  5:45223/1‑147 (MQ=255)                                                                         tttAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGt                                                                                                >  7:229217/1‑126 (MQ=255)                                                                         tttAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGt                                                                                                <  8:229217/126‑1 (MQ=255)                                                                              cAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGAc                                                                                                              <  3:223435/107‑1 (MQ=255)                                                                              cAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGAc                                                                                                              >  4:223435/1‑107 (MQ=255)                                                                                    aGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGa                                                               <  5:294590/148‑1 (MQ=255)                                                                                               gATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGACGTGAGGGCGATACTTc                                                   >  5:250523/1‑149 (MQ=255)                                                                                                    tGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCt                                              >  4:137922/1‑149 (MQ=255)                                                                                                    tGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCt                                              >  2:69572/1‑149 (MQ=255)                                                                                                       aGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGa                                           <  3:137922/149‑1 (MQ=255)                                                                                                            tAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCt                                       >  3:112736/1‑148 (MQ=255)                                                                                                                gAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAg                                  >  3:203408/1‑149 (MQ=255)                                                                                                                 aaTTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAg                                  >  1:135127/1‑148 (MQ=255)                                                                                                                      tGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAgg                            >  3:203134/1‑149 (MQ=255)                                                                                                                        gaaatGAAGAGGTTGCGAAGGTCAATCTCAAGCTTCGCTTGGAATATATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGGTACTTCCGGCTTGATGCTAAAAGTCAAgg                            <  6:45223/146‑1 (MQ=255)                                                                                                                        cAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAGGt                           >  3:305230/1‑148 (MQ=255)                                                                                                                         aaaaGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAGGTa                          >  2:170544/1‑148 (MQ=255)                                                                                                                         aaaaGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAGGTa                          >  8:4362/1‑148 (MQ=255)                                                                                                                          aaaGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAGGTaaa                        <  4:203408/149‑1 (MQ=255)                                                                                                                          aaaGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAGGTa                          >  4:340065/1‑147 (MQ=255)                                                                                                                           aagaagAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTaaaa                                   <  6:366776/137‑1 (MQ=255)                                                                                                                           aagaagAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTaaaa                                   >  5:366776/1‑137 (MQ=255)                                                                                                                             gcgaagCTTGCGATTGACCTTCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAGGTAAAGGa                     >  2:359808/5‑149 (MQ=255)                                                                                                                             gaagaGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAGGTAAAGGa                     <  7:4362/149‑1 (MQ=255)                                                                                                                                        gAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAGGTAAAGGACCTTTCGGCaa          <  4:305230/149‑1 (MQ=255)                                                                                                                                          aGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAGGTAAAGGACCTTTCGGc            >  8:100371/1‑145 (MQ=255)                                                                                                                                              cAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAGGTAAAGGACCTTTCGGCAAGATGTc    >  8:214699/1‑149 (MQ=255)                                                                                                                                              cAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAGGTAAAGGACCTTTCGGCAAGATGTc    >  8:206761/1‑149 (MQ=255)                                                                                                                                               aaTCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAGGTAAAGGACCTTTCGGCAAGATGTc    <  1:170544/148‑1 (MQ=255)                                                                                                                                               aaTCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAGGTAAAGGACCTTTCGGCAAGATGTCt   >  8:131980/1‑149 (MQ=255)                                                                                                                                                aTCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAGGTAAAGGACCTTTCGGCAAGATGTCtt  <  2:135127/149‑1 (MQ=255)                                                                                                                                                aTCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAGGTAAAGGACCTTTCGGCAAGATGTCtt  <  1:271256/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                     ACTTGTCCTGTCTGTTGTAACGCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAGGTAAAGGACCTTTCGGCAAGATGTCTT  >  NC_002947/3411519‑3411810

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