Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,411,661	A→C	10.8%	N59H (AAT→CAT) ‡	PP_5551 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,411,661	0	A	C	10.8%	183.4 / 13.0	65	N59H (AAT→CAT) ‡	PP_5551	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (33/25);  new base C (3/4);  total (36/29)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

ATTGCACTTGTCCTGTCTGTTGTAACGCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAGGTAAAGGACCTTTCGGCAAGATGTCT  >  NC_002947/3411514‑3411809                                                                                                                                                    |                                                                                                                                                     aTTGCACTTGTCCTGTCTGTTGTAACGCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCaa                                                                                                                                                     >  1:158771/1‑149 (MQ=255)           tCCTGTCTGTTGTAACGCTGGGCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGATTGACCTTCGCAAGCtc                                                                                                                                           <  6:342574/149‑2 (MQ=255)             ctgtctgtTGTAACGCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGATTGACCTTCGCAAGCtct                                                                                                                                          <  3:249306/148‑3 (MQ=255)             ctgtctgtTGTAACGCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCg                                                                                                                                         >  4:142265/1‑149 (MQ=255)              tgtctgtTGTAACGCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGATTGACCTTCGCAAGc                                                                                                                                             >  2:94748/1‑144 (MQ=255)              tgtctgtTGTAACGCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGATTGACCTTCGCAAGc                                                                                                                                             <  1:94748/144‑1 (MQ=255)                tctgtTGTAACGCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCt                                                                                                                                       <  5:296572/148‑1 (MQ=255)                  tgttgtAACGCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCtt                                                                                                                                      <  8:309873/147‑1 (MQ=255)                  tgttgtAACGCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCtt                                                                                                                                      >  7:309873/1‑147 (MQ=255)                  tgttgtAACGCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTg                                                                                                                                     >  6:43514/1‑148 (MQ=255)                  tgttgtAACGCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTg                                                                                                                                     >  3:92412/1‑148 (MQ=255)                       tAACGCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGaa                                                                                                                                  >  2:259840/1‑146 (MQ=255)                         aCGCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGa                                                                                                                              <  7:283718/148‑1 (MQ=255)                          cGCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGAtt                                                                                                                            <  5:276229/149‑1 (MQ=255)                                  gCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTc                                                                                                                        >  3:139843/1‑145 (MQ=255)                                   cTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGa                                                                                                                                   <  6:236680/133‑1 (MQ=255)                                   cTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGa                                                                                                                                   >  5:236680/1‑133 (MQ=255)                                     cAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGAtt                                                                                                                  <  4:92412/148‑1 (MQ=255)                                        aaTAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAg                                                                                                                               >  3:91579/1‑132 (MQ=255)                                        aaTAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAg                                                                                                                               <  4:91579/132‑1 (MQ=255)                                                 agcaagcGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGa                                                                                                                              <  6:243795/124‑1 (MQ=255)                                                 agcaagcGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGa                                                                                                                              >  5:243795/1‑124 (MQ=255)                                                  gcaagcGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAAATCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTg                                                                                                    >  5:328877/1‑149 (MQ=255)                                                   caagcGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTgg                                                                                                   <  3:142265/149‑1 (MQ=255)                                                      gcGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTg                                                                                                >  1:133067/1‑149 (MQ=255)                                                         aTCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCa                                                                                             <  6:328877/149‑1 (MQ=255)                                                             gCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGAttt                                                                                                                 >  3:215727/1‑125 (MQ=255)                                                             gCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGAttt                                                                                                                 <  4:215727/125‑1 (MQ=255)                                                              cGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCtt                                                                                        <  5:43514/149‑1 (MQ=255)                                                               ggCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCtt                                                                                        >  2:191070/1‑148 (MQ=255)                                                               ggCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTc                                                                                       <  3:231367/149‑1 (MQ=255)                                                                            aaTTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAg                                                                          >  7:157049/1‑149 (MQ=255)                                                                             aTTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGt                                                                         <  1:191070/149‑1 (MQ=255)                                                                                tAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGa                                                                       <  5:191252/148‑1 (MQ=255)                                                                                      gTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGt                                                                >  4:58374/1‑149 (MQ=255)                                                                                       tCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGt                                                                >  5:64072/1‑148 (MQ=255)                                                                                         aGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAg                                                             >  7:169989/1‑149 (MQ=255)                                                                                              gcgAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGgcgc                                                        >  2:153147/1‑149 (MQ=255)                                                                                                gAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTa                                                      >  5:22536/1‑149 (MQ=255)                                                                                                   ggATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACtt                                                   <  3:58374/149‑1 (MQ=255)                                                                                                   ggATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACt                                                    <  2:158771/148‑1 (MQ=255)                                                                                                    gATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTc                                                  >  4:39697/1‑149 (MQ=255)                                                                                                      tCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCg                                                >  6:220273/1‑149 (MQ=255)                                                                                                         tGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCt                                             >  2:338365/1‑149 (MQ=255)                                                                                                          gAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCtt                                            >  4:46918/1‑149 (MQ=255)                                                                                                           aaGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCtt                                            <  6:22536/148‑1 (MQ=255)                                                                                                             gATCTAATGAATTCTCCAAAAGAAGAGCTTGCGAAGGTCATTCGAAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACTCGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGAt                                         <  6:64072/149‑1 (MQ=255)                                                                                                                 tAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGc                                               >  4:301526/1‑139 (MQ=255)                                                                                                                 tAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGc                                               <  3:301526/139‑1 (MQ=255)                                                                                                                    aGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTaaaa                                  <  1:153147/149‑1 (MQ=255)                                                                                                                    aGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTaaaa                                  <  5:220273/149‑1 (MQ=255)                                                                                                                       aTTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGt                                >  6:318241/1‑148 (MQ=255)                                                                                                                           tGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAgg                           >  8:302896/1‑149 (MQ=255)                                                                                                                           tGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAgg                           >  7:20497/1‑149 (MQ=255)                                                                                                                              aaaaGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGc                                       <  6:362398/134‑1 (MQ=255)                                                                                                                               caagcgaagCTTGCGATTGACCTTCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAGGTa                         >  4:173456/8‑147 (MQ=255)                                                                                                                                aagaagAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTaaaa                                  <  1:261069/137‑1 (MQ=255)                                                                                                                                aagaagAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTaaaa                                  >  2:261069/1‑137 (MQ=255)                                                                                                                                  gaagaGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAGGTAAAGGa                    >  2:109180/1‑149 (MQ=255)                                                                                                                                    gaagCTTGCGATTGACCTTCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAg                            <  3:143369/137‑1 (MQ=255)                                                                                                                                    gaagCTTGCGATTGACCTTCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAg                            >  4:143369/3‑139 (MQ=255)                                                                                                                                        cTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAGGTAAAGGACCTTTc              >  2:168218/1‑149 (MQ=255)                                                                                                                                        cTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAGGTAAAGGACCTTTc              >  6:32525/1‑149 (MQ=255)                                                                                                                                         ttGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAGGTAAAGGACCTTTCg             >  6:140385/1‑149 (MQ=255)                                                                                                                                           gCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAGGTAAAGGACCTTTCg             <  8:169989/147‑1 (MQ=255)                                                                                                                                             gAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAGGTAAAGGACCTTTCGGc           >  2:315836/1‑147 (MQ=255)                                                                                                                                                    aaTCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGa                                                                              >  6:40079/1‑73 (MQ=255)                                                                                                                                                    aaTCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGa                                                                              <  5:40079/73‑1 (MQ=255)                                                                                                                                                    aaTCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAGGTAAAGGACCTTTCGGCAAGATGTCt  >  7:238600/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     ATTGCACTTGTCCTGTCTGTTGTAACGCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAGGTAAAGGACCTTTCGGCAAGATGTCT  >  NC_002947/3411514‑3411809

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