Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000731	25,669	(A)7→6	coding (322/333 nt)	USA300HOU_pUSA300HOUMR0031 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000731	25,663	0	A	.	98.5%	312.0 / ‑2.7	66	coding (316/333 nt)	USA300HOU_pUSA300HOUMR0031	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base . (35/30);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTGTATGTTATTTACACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCAATAGAAAACAA  >  CP000731/25523‑25798                                                                                                                                             |                                                                                                                                        attgtatGTTATTTACACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATa                                                                                                                                         >  2:127365/1‑141 (MQ=255)  ttgtatGTTATTTACACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaa                                                                                                                                        >  1:11337/1‑141 (MQ=255)   tgtatGTTATTTACACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaa                                                                                                                                       <  1:97961/141‑1 (MQ=255)   tgtatGTTATTTACACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGAATCTCGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaa                                                                                                                                       <  1:212231/141‑1 (MQ=255)       tGTTATTTACACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaaat                                                                                                                                   >  2:240772/1‑140 (MQ=255)           aTTTACACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATAc                                                                                                                              >  1:99631/1‑141 (MQ=255)            tttACACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACa                                                                                                                             <  1:54903/141‑1 (MQ=255)             ttACACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACAt                                                                                                                            >  1:40265/1‑141 (MQ=255)             ttACACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACAt                                                                                                                            >  1:227377/1‑141 (MQ=255)             ttACACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACAt                                                                                                                            >  1:198465/1‑141 (MQ=255)                 acTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATCATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTaa                                                                                                                        >  2:161406/1‑141 (MQ=255)                 acTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTaa                                                                                                                        >  2:134525/1‑141 (MQ=255)                  cTAGTGATTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTaaa                                                                                                                       <  2:291570/141‑1 (MQ=255)                   tAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAg                                                                                                                      <  1:232988/141‑1 (MQ=255)                     gTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGt                                                                                                                    <  1:269264/141‑1 (MQ=255)                     gTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGt                                                                                                                    >  2:239613/1‑141 (MQ=255)                     ctccTTAAATTTACAGGTTTGTCACCAAGTAGTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGt                                                                                                                    <  1:216478/138‑1 (MQ=255)                           aaaTTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAAtt                                                                                                              >  2:150266/1‑141 (MQ=255)                            aaTTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTa                                                                                                             >  2:213329/1‑141 (MQ=255)                             aTTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTaa                                                                                                            <  2:185340/141‑1 (MQ=255)                              tttACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAAt                                                                                                           >  2:578/1‑141 (MQ=255)                                tACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGa                                                                                                         <  1:56949/141‑1 (MQ=255)                                tACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGa                                                                                                         <  1:238681/141‑1 (MQ=255)                                 aCAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGaa                                                                                                        >  1:134269/1‑141 (MQ=255)                                  cAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGaaa                                                                                                       >  1:125793/1‑141 (MQ=255)                                   aGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGaaaa                                                                                                      >  1:248911/1‑141 (MQ=255)                                   aGGTTTGACACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGaaaa                                                                                                      >  1:241685/1‑141 (MQ=255)                                     gTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAAt                                                                                                    <  1:82767/141‑1 (MQ=255)                                     gTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAAt                                                                                                    <  1:140253/141‑1 (MQ=255)                                      tttGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTCATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATa                                                                                                   >  1:151836/1‑141 (MQ=255)                                        tGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGa                                                                                                 >  1:2912/1‑141 (MQ=255)                                        tGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGa                                                                                                 <  2:8559/141‑1 (MQ=255)                                         gTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGaa                                                                                                >  2:210206/1‑141 (MQ=255)                                             ccAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTag                                                                                            >  2:236047/1‑141 (MQ=255)                                              cAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTaga                                                                                           >  1:137545/1‑141 (MQ=255)                                                     tGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAAtata                                                                                    >  1:223472/1‑141 (MQ=255)                                                     tGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAAtata                                                                                    >  1:19049/1‑141 (MQ=255)                                                      gAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAg                                                                                   >  1:232210/1‑141 (MQ=255)                                                       aaGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGa                                                                                  >  2:177084/1‑141 (MQ=255)                                                                ggTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACaaaaa                                                                         <  1:203815/141‑1 (MQ=255)                                                                     tctcAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAAt                                                                    >  2:43523/1‑141 (MQ=255)                                                                       tcAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATAAAGAAAAAAAAAAGGAAttt                                                                  <  2:99631/141‑1 (MQ=255)                                                                             tCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAg                                                            >  2:18103/1‑141 (MQ=255)                                                                             tCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAg                                                            <  1:173657/141‑1 (MQ=255)                                                                              ccAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGa                                                           >  1:10469/1‑141 (MQ=255)                                                                              ccAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGa                                                           >  1:49837/1‑141 (MQ=255)                                                                                  aaaTTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATAATAAAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACt                                                       <  2:42605/141‑1 (MQ=255)                                                                                   aaTTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACtt                                                      >  1:172298/1‑141 (MQ=255)                                                                                     ttttCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTa                                                    <  2:241685/141‑1 (MQ=255)                                                                                     ttttCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATAAAAAAAACAAAAAGGAATTTCTCAAGAAACTTTa                                                    <  2:248911/141‑1 (MQ=255)                                                                                     ttttCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTAAAAATAAAGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTa                                                    <  2:120511/141‑1 (MQ=255)                                                                                     aattatttATAATTGTAGAACTATTTCATATGCAATTAAATTACTAAACAACAAAT‑AAAAAATATACATCTAAAGGATAAATTAATGAAAAATAGAATTAGAGAAAAAAAAAAAAAAAAAGGAATTTCTCAAGAAACTTTa                                                    <  2:120549/136‑1 (MQ=255)                                                                                         cTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTa                                                <  1:94816/141‑1 (MQ=255)                                                                                         cTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTa                                                >  2:52657/1‑141 (MQ=255)                                                                                            tataATTGTAGCACTATTTCTTATGCTATTTAAATACTATACTACAAA‑AAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAAAATAAAAAAAAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAg                                             <  2:273649/141‑1 (MQ=255)                                                                                                aTTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAAtt                                         <  2:63509/141‑1 (MQ=255)                                                                                                aTTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATATAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAAtt                                         <  2:63502/141‑1 (MQ=255)                                                                                                   gTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTaaa                                      <  1:38770/141‑1 (MQ=255)                                                                                                                 tATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCaaa                        <  1:135400/141‑1 (MQ=255)                                                                                                                       gTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTa                  >  2:68119/1‑141 (MQ=255)                                                                                                                       gTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTa                  >  2:248851/1‑141 (MQ=255)                                                                                                                           aaTTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGc              <  1:150266/141‑1 (MQ=255)                                                                                                                               actatactaCAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCAATa          >  1:279145/1‑141 (MQ=255)                                                                                                                                 tatactaCAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCAATAGa        >  1:84735/1‑141 (MQ=255)                                                                                                                                 tatactaCAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCAATAGa        >  1:300346/1‑141 (MQ=255)                                                                                                                                    actacAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCAATAGaaaa     <  2:49837/141‑1 (MQ=255)                                                                                                                                    actacAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCAATAGaaaa     <  2:10469/141‑1 (MQ=255)                                                                                                                                    actacAAAT‑AAAAAATAGACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCAATAGaaaa     <  2:245458/141‑1 (MQ=255)                                                                                                                                     ctacAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCAATAGAAAAc    <  1:287659/141‑1 (MQ=255)                                                                                                                                     ctacAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCAATAGAAAAc    <  1:112842/141‑1 (MQ=255)                                                                                                                                       acAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCAATAGAAAACaa  >  2:161218/1‑141 (MQ=255)                                                                                                                                             |                                                                                                                                        ATTGTATGTTATTTACACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCAATAGAAAACAA  >  CP000731/25523‑25798

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