Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,535,071	A→G	5.8%	L130L (TTA→TTG)	PP_4022 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,535,071	0	A	G	5.8%	256.8 / 4.1	86	L130L (TTA→TTG)	PP_4022	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (40/41);  new base G (3/2);  total (43/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.99e-01

CTCATGAGATTCCACTCGCCTGACAGCGATAGCCCTTTCGACGGGGGTGGGCT‑CAACGCTTATGCATATTGCTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGAACAGTGTTATTAAAAAAGAACTGAAAACCCTTTCTGACACAGGTAAAGCCCT  >  NC_002947/4534931‑4535217                                                                                                                                              |                                                                                                                                                    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                                                                                                          >  8:250739/1‑148 (MQ=255)                                    tttCGACGGGGGTGGGCT‑CAACGCTTATGCATATTGCTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCtt                                                                                                          >  1:239897/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                                                                                    <  4:207259/99‑1 (MQ=255)                                                                                                            gTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCtg                                                                                    >  3:207259/1‑99 (MQ=255)                                                                                                            gTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGAACAGTGTTATTAAAAAAGa                                   >  4:20756/1‑148 (MQ=255)                                                                                                              cATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGAACAGTGTTATTAAAAAAGAACt                                >  8:168009/1‑149 (MQ=255)                                                                                                               aTGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGAACAGTGTTATTaaaa                                       >  5:100438/1‑141 (MQ=255)                                                                                                               aTGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGAACAGTGTTATTaaaa                                       <  6:100438/141‑1 (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                   <  1:80606/130‑1 (MQ=255)                                                                                                                                              aCCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGAACAGTGTTATTAAAAAAGAACTGAAAACCCTTTCt                   >  2:80606/1‑130 (MQ=255)                                                                                                                                              |                                                                                                                                                    CTCATGAGATTCCACTCGCCTGACAGCGATAGCCCTTTCGACGGGGGTGGGCT‑CAACGCTTATGCATATTGCTCAGGAGACCCCATTAACTACCATGACCCTTCAGGTCATGGTCCAGTTTTGTCCAAAAAATCTTTTTTACCAAAAAAGAAAACGATCATGCA‑GCGACGCGAGGCGGCGCTTAAAGCAATCAATGAGAACCCTGTGCTCGAAGCGTTCCAAGGCATTGCCGAGAACAGTGTTATTAAAAAAGAACTGAAAACCCTTTCTGACACAGGTAAAGCCCT  >  NC_002947/4534931‑4535217

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