Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,535,073	C→T	5.7%	P131L (CCA→CTA) ‡	PP_4022 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,535,073	0	C	T	5.7%	259.8 / 3.5	87	P131L (CCA→CTA) ‡	PP_4022	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (40/42);  new base T (3/2);  total (43/44)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.63e-02

CTCATGAGATTCCACTCGCCTGACAGCGATAGCCCTTTCGACGGGGGTGGGCT‑CAACGCTTATGCATATTGCTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGAACAGTGTTATTAAAAAAGAACTGAAAACCCTTTCTGACACAGGTAAAGCCCTGAAC  >  NC_002947/4534931‑4535221                                                                                                                                                |                                                                                                                                                      cTCATGAGATTCCACTCGCCTGACAGCGATAGCCCTTTCGACGGGGGTGGGCT‑CAACGCTTATGCATATTGCTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTGGTaaaaaa                                                                                                                                                 >  3:340352/1‑149 (MQ=255)  tCATGAGATTCCACTCGCCTGACAGCGATAGCCCTTTCGACGGGGGTGGGCT‑CAACGCTTATGCATATTGCTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACaaata                                                                                                                                                   >  5:264327/1‑144 (MQ=255)   cATGAGATTCCACTCGCCTGACAGCGATAGCCCTTTCGACGGGGGTGGGCT‑CAACGCTTATGCATATTGCTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAg                                                                                                                                                >  7:4954/1‑148 (MQ=255)      gagaTTCCACTCGCCTGACAGCGATAGCCCTTTCGACGGGGGTGGGCT‑CAACGCTTATGCATATTGCTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAGGaaaa                                                                                                                                            >  2:170281/1‑149 (MQ=255)        gaTTCCACTCGCCTGACAGCGATAGCCCTTTCGACGGGGGTGGGCT‑CAACGCTTATGCATATTGCTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAAc                                                                                                                                           <  8:119352/148‑1 (MQ=255)         aTTCCACTCGCCTGACAGCGATAGCCCTTTCGACGGGGGTGGGCT‑CAACGCTTATGCATATTGCTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTGGTAAAAAAGatttttt                                                                                                                                         <  4:340352/149‑7 (MQ=255)          ttCCACTCGCCTGACAGCGATTGCCCTTTCGACGGGGGTGGGCT‑CAACGCTTATGCATATTGCTCAGGAGACCCCATTAACTACCATCACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGa                                                                                                                                         <  2:280086/148‑1 (MQ=255)             cACTCGCCTGACAGCGATAGCCCTTTCGACGGGGGTGGGCT‑CAACGCTTATGCATATTGCTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCAt                                                                                                                                     <  8:4954/149‑1 (MQ=255)              aCTCGCCTGACAGCGATAGCCCTTTCGACGGGGGTGGGCT‑CAACGCTTATGCATATTGCTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATg                                                                                                                                    >  4:21379/1‑149 (MQ=255)               cTCGCCTGACAGCGATAGCCCTTTCGACGGGGGTGGGCT‑CAACGCTTATGCATATTGCTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATATTTTTTACCAAAAAAGAAAACGATCATg                                                                                                                                    >  7:224882/1‑148 (MQ=255)                  gCCTGACAGCGATAGCCCTTTCGACGGGGGTGGGCT‑CAACGCTTATGCATATTGCTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑Gc                                                                                                                               >  3:297809/1‑149 (MQ=255)                     tGACAGCGATAGCCCTTTCGACGGGGGTGGGCT‑CAACGCTTATGCATATTGCTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGAc                                                                                                                            >  7:106032/1‑149 (MQ=255)                       aCAGCGATAGCCCTTTCGACGGGGGTGGGCT‑CAACGCTTATGCATATTGCTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGAcgc                                                                                                                          >  2:336403/1‑149 (MQ=255)                           cGATAGCCCTTTCGACGGGGGTGGGCT‑CAACGCTTATGCATATTGCTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTGGTAAAAAAGatttttt                                                                                                                                         >  1:7955/1‑125 (MQ=255)                           cGATAGCCCTTTCGACGGGGGTGGGCT‑CAACGCTTATGCATATTGCTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTGGTAAAAAAGatttttt                                                                                                                                         <  2:7955/131‑7 (MQ=255)                           cGATAGCCCTTTCGACGGGGGTGGGCT‑CAACGCTTATGCATATTGCTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAGGAAAACGATCATGCA‑GCGACGCGAgg                                                                                                                      <  1:170281/149‑1 (MQ=255)                            gATAGCCCTTTCGACGGGGGTGGGCTACAACGCTTATGCATATTGCTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAgg                                                                                                                      >  1:188018/1‑149 (MQ=255)                             aTAGCCCTTTCGACGGGGGTGGGCT‑CAACGCTTATGCATATTGCTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAggcg                                                                                                                    >  5:3238/1‑149 (MQ=255)                                gCCCTTTCGACGGGGGTGGGCT‑CAACGCTTATGCATATTTCTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGgcggcg                                                                                                                 <  8:224882/149‑1 (MQ=255)                                    tttCGACGGGGGTGGGCT‑CAACGCTTATGCATATTGCTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑TCGACGCGAGGCGTCGCtt                                                                                                              >  8:245744/1‑148 (MQ=255)                                    tttCGACGGGGGTGGGCT‑CAACGCTTATGCATATTGCTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCtt                                                                                                              >  1:239897/1‑148 (MQ=255)                                    tttCGACGGGGGTGGGCT‑CAACGCTTATGCATATTGCTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCtt                                                                                                              >  8:250739/1‑148 (MQ=255)                                       cGACGGGGGTGGGCT‑CAACGCTTATGCATATTGCTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTaaa                                                                                                           <  5:183649/148‑1 (MQ=255)                                        gACGGGGGTGGGCT‑CAACGCTTATGCATATTGCTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAg                                                                                                          >  4:203193/1‑148 (MQ=255)                                        gACGGGGGTGGGCT‑CAACGCTTATGCATATTGCTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAg                                                                                                          <  4:139957/148‑1 (MQ=255)                                         aCGGGGGTGGGCT‑CAACGCTTATGCATATTGCTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGca                                                                                                        <  3:21379/149‑1 (MQ=255)                                           gggggTGGGCT‑CAACGCTTATGCATATTGCTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGcaat                                                                                                      <  1:336403/149‑1 (MQ=255)                                           gggggTGGGCT‑CAACGCTTATGCATATTGCTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGcaat                                                                                                      <  8:274557/149‑1 (MQ=255)                                                ttggat‑CAACTATTATGCATATTGCTCAGGAGAAAACATTAACTACCATGACCCTTCAGGTCATGGTACAGTTTTGTACAAAAAATCTTTTTTACCAAAAAAGAAAACTATCATGCATGCGACGCGAGGCGGCGCTTAAAGcaatcaa                                                                                                   <  1:215204/143‑1 (MQ=255)                                                tGGGCT‑CAACGCTTATGCATATTGCTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATg                                                                                                 <  3:11932/149‑1 (MQ=255)                                                  ggCT‑CAACGCTTATGCATATTGCTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATga                                                                                                <  3:256791/148‑1 (MQ=255)                                                    cT‑CAACGCTTATGCATATTGCTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATg                                                                                                                                    >  3:276167/1‑111 (MQ=255)                                                    cT‑CAACGCTTATGCATATTGCTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATg                                                                                                                                    <  4:276167/111‑1 (MQ=255)                                                    cT‑CAACGCTTATGCATATTGCTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGaa                                                                                             >  3:269275/1‑149 (MQ=255)                                                     t‑CAACGCTTATGCATATTGCTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGaaaa                                                                                                                                            <  7:334590/102‑1 (MQ=255)                                                     t‑CAACGCTTATGCATATTGCTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGaaaa                                                                                                                                            >  8:334590/1‑102 (MQ=255)                                                       cAACGCTTATGCATATTGCTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAAcc                                                                                           >  2:200467/1‑149 (MQ=255)                                                                tGCATATTGCTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTc                                                                                   <  4:218425/148‑1 (MQ=255)                                                                 gCATATTGCTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTAACAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCg                                                                                  >  7:113659/1‑148 (MQ=255)                                                                   atatTGCTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGaa                                                                                <  2:188018/148‑1 (MQ=255)                                                                   atatTGCTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GAGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAg                                                                               >  1:345301/1‑149 (MQ=255)                                                                    tatTGCTCAGGAGACCCCATAAACTACCATGACCCTTCAGTTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGc                                                                              <  7:323065/149‑1 (MQ=255)                                                                         cTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAggc                                                                                                                     >  3:292850/1‑105 (MQ=255)                                                                         cTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAggc                                                                                                                     <  4:292850/105‑1 (MQ=255)                                                                         cTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCAAAGCGTTc                                                                          <  2:239897/148‑1 (MQ=255)                                                                              gagaCCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAgg                                                                     >  7:245849/1‑148 (MQ=255)                                                                                 aCCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCAtt                                                                 >  7:173226/1‑149 (MQ=255)                                                                                   cccATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGc                                                               >  8:317238/1‑149 (MQ=255)                                                                                     cATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGcc                                                              >  4:141158/1‑148 (MQ=255)                                                                                      aTTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCga                                                            <  1:200467/149‑1 (MQ=255)                                                                                       ttAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCgag                                                           <  3:203193/149‑1 (MQ=255)                                                                                           cTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCCGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGaa                                                         >  4:80630/1‑147 (MQ=255)                                                                                              ccATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGAACAgtg                                                    <  4:297809/149‑1 (MQ=255)                                                                                              ccATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGATCGAAGCGTTCCAAGGCATTGCCGAGAACAgtg                                                    >  7:230676/1‑149 (MQ=255)                                                                                               cATGACCCTACAGGTCATGGTCCAGTTTAGTCCAAAAAATCTTTTATACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGAACAgtgt                                                   <  6:264327/149‑1 (MQ=255)                                                                                                aTGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGAACAGTGtt                                                  <  6:3238/149‑1 (MQ=255)                                                                                                 tGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGAACAGTGtt                                                  >  4:302966/1‑148 (MQ=255)                                                                                                   aCCCTTCAGGTCATTGTCAAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCTACGCGAGGCGGCGCTTAAAGCAATCAATTAGAACGCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGAACAGTGttat                                                <  7:97365/148‑1 (MQ=255)                                                                                                      cTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGcaatcaat                                                                                                  >  1:26151/1‑95 (MQ=255)                                                                                                      cTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGcaatcaat                                                                                                  <  2:26151/95‑1 (MQ=255)                                                                                                         cAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGAACAGTGTTATTaaaaa                                          <  8:245849/148‑1 (MQ=255)                                                                                                          aGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGAACAGTGTTATTaaaaaa                                         <  7:250739/148‑1 (MQ=255)                                                                                                          aGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGAACAGTGTTATTaaaaaa                                         <  7:245744/148‑1 (MQ=255)                                                                                                            gTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCtg                                                                                        >  3:207259/1‑99 (MQ=255)                                                                                                            gTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCtg                                                                                        <  4:207259/99‑1 (MQ=255)                                                                                                            gTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGAACAGTGTTATTAAAAAAGa                                       >  4:20756/1‑148 (MQ=255)                                                                                                              cATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGAACAGTGTTATTAAAAAAGAACt                                    >  8:168009/1‑149 (MQ=255)                                                                                                               aTGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGAACAGTGTTATTaaaa                                           <  6:100438/141‑1 (MQ=255)                                                                                                               aTGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGAACAGTGTTATTaaaa                                           >  5:100438/1‑141 (MQ=255)                                                                                                               aTGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGAACAGTGTTATTAAAAAAGAACTg                                   <  3:141158/149‑1 (MQ=255)                                                                                                                  gTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGAACAGTGTTATTAAAAAAGAACTGaa                                 <  3:80630/148‑1 (MQ=255)                                                                                                                   tCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGAACAGTGTTATTAAAAAAGAACTGaaaa                               >  4:152979/1‑149 (MQ=255)                                                                                                                       gTTTTCTCCAAAAAATCATTTTTACCAAAAAAGAAAACCATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGAACAGTGTTATTAAAAAAGAACTGAAAACCCt                           <  8:230676/149‑1 (MQ=255)                                                                                                                             tCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGAACAGTGTTATTAAAAAAGAACTGAAAACCCTTTCTGa                     <  7:325249/149‑1 (MQ=255)                                                                                                                              ccAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGAACAGTGTTATTAAAAAAGAACTGAAAACCCTTTCTGa                     >  3:115545/1‑148 (MQ=255)                                                                                                                                aaaaaaTCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGAACAGTGTTATTAAAAAAGAACTGAAAACCCTTTCTGaca                   >  2:195087/1‑148 (MQ=255)                                                                                                                                   ttttCTTTTTTGGTAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGCGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGAACAGTGTTATTAAAAAAGAACTGAAAAACCTTTCTGACCCAgg               >  7:260544/4‑149 (MQ=255)                                                                                                                                       cTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGAACAGTGTTATTAAAAAAGAACTGaaaa                               >  5:196960/1‑129 (MQ=255)                                                                                                                                       cTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGAACAGTGTTATTAAAAAAGAACTGaaaa                               <  6:196960/129‑1 (MQ=255)                                                                                                                                        ttttttACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGAACAGTGTTATTAAAAAAGAACTGAAAACCCTTTCTGACACAGGTAAAg          >  3:328171/1‑149 (MQ=255)                                                                                                                                          ttttACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGAACAGTGTTATTAAAAAAGAACTGAAAACCCTTTCTGACACAGGTaaa           <  7:317238/146‑1 (MQ=255)                                                                                                                                            ttACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGAACAGTGTTATTAAAAAAGAACTGAAAACCCTTTCTGACACAGGTAAAGccc       <  6:335299/148‑1 (MQ=255)                                                                                                                                             tACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGAACAGTGTTATTAAAAAAGAACTGAAAACCCTTTCTg                      <  6:10178/132‑1 (MQ=255)                                                                                                                                             tACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGAACAGTGTTATTAAAAAAGAACTGAAAACCCTTTCTg                      >  5:10178/1‑132 (MQ=255)                                                                                                                                             tACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGAACAGTGTTATTAAAAAAGAACTGAAAACCCTTTCTGACACAGGTAAAGCCCt      >  8:27338/1‑148 (MQ=255)                                                                                                                                             tACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGAACAGTGTTATTAAAAAAGAACTGAAAACCCTTTCTGACACAGGTAAAGCCCt      <  1:195087/148‑1 (MQ=255)                                                                                                                                              aCCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGAACAGTGTTATTAAAAAAGAACTGAAAACCCTTTCt                       >  2:80606/1‑130 (MQ=255)                                                                                                                                              aCCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGAACAGTGTTATTAAAAAAGAACTGAAAACCCTTTCt                       <  1:80606/130‑1 (MQ=255)                                                                                                                                                cAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGAACAGTGTTATTaaaaaa                                         <  3:81294/110‑1 (MQ=255)                                                                                                                                                cAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGAACAGTGTTATTaaaaaa                                         >  4:81294/1‑110 (MQ=255)                                                                                                                                                cAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGAACAGTGTTATTAAAAAAGAACTGAAAACCCTTTCTGACACAGGTAAAGCCCTGAAc  >  6:62428/1‑149 (MQ=255)                                                                                                                                                cAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGAACAGTGTTATTAAAAAAGAACTGAAAACCCTTTCTGACACAGGTAAAGCCCTGAAc  >  8:4965/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                      CTCATGAGATTCCACTCGCCTGACAGCGATAGCCCTTTCGACGGGGGTGGGCT‑CAACGCTTATGCATATTGCTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGAACAGTGTTATTAAAAAAGAACTGAAAACCCTTTCTGACACAGGTAAAGCCCTGAAC  >  NC_002947/4534931‑4535221

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