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	.	97.1%	295.1 / 0.6	68	coding (316/333 nt)	USA300HOU_pUSA300HOUMR0031	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base . (40/26);  total (40/28)
Fisher's exact test for biased strand distribution p-value = 1.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGTATGTTATTTACACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCAATAGAAAACAATA  >  CP000731/25524‑25800                                                                                                                                            |                                                                                                                                          ttgtatGTTATTTACACTATTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATa                                                                                                                                           <  1:10495/140‑1 (MQ=255) ttgtatGTTATTTACACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATa                                                                                                                                           <  2:127334/140‑1 (MQ=255) ttgtatGTTATTTACACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATa                                                                                                                                           >  2:263210/1‑140 (MQ=255) ttgtatGTTATTTACACTAGTGCTTAAAATTACACGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATa                                                                                                                                           <  2:86171/140‑1 (MQ=255) gtgtatGTTATTTACACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATa                                                                                                                                           <  1:122791/139‑1 (MQ=255)   gtatGTTATTTACACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACACATaaa                                                                                                                                         >  2:79974/1‑140 (MQ=255)   gtatGTTATTTACACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaa                                                                                                                                         >  1:32394/1‑140 (MQ=255)            ttACACTAGTGCTTAGATGGACATGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACa                                                                                                                               >  2:233648/1‑140 (MQ=255)            ttACACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAACT‑AAAAAATATACa                                                                                                                               >  2:126258/1‑140 (MQ=255)             tACACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACAt                                                                                                                              >  1:121357/1‑140 (MQ=255)              acacTAGTGTTTAAAGTTACAGGGTTGTCACCAAGTATTGAAGGATATTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATc                                                                                                                             <  2:115662/140‑1 (MQ=255)                acTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTACACTACACCT‑ACAAAATATACATCtc                                                                                                                           >  1:213773/1‑139 (MQ=255)                 cTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTaa                                                                                                                          >  1:111679/1‑140 (MQ=255)                    gTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAgg                                                                                                                       >  2:154744/1‑140 (MQ=255)                     tGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATTCAGAAAGTTTCTTTATAAAAGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGt                                                                                                                      <  2:146351/140‑1 (MQ=255)                     tGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGt                                                                                                                      <  1:9639/140‑1 (MQ=255)                     tGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGt                                                                                                                      <  1:207285/140‑1 (MQ=255)                       cTTAAATTTACAGGTTTGTCACCAGGTATTGAAGGATCTTGTTTTTCTCAATATCCATAAATTTTCTTTATAATTTTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTa                                                                                                                    >  1:33489/1‑140 (MQ=255)                       cTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCACTATCCAGAAATTTTCTTTATAATTGTAGTACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAACTATACATCTAAATGTTa                                                                                                                    >  2:113256/1‑140 (MQ=255)                       cTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTa                                                                                                                    >  2:48968/1‑140 (MQ=255)                       cTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTAAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATTTAACGGTTa                                                                                                                    >  1:44365/1‑140 (MQ=255)                        ttAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGttca                                                                                                                   >  1:62852/1‑138 (MQ=255)                        ttAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTaa                                                                                                                   <  1:158067/140‑1 (MQ=255)                             tttACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTCTCTTTGTAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGTTTAAATTaa                                                                                                              >  1:168797/1‑140 (MQ=255)                             tttACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTATCAATATCAAGTATTTTTCTTTATAATGGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTaa                                                                                                              >  1:138742/1‑140 (MQ=255)                             tttACAGGTTTGTCACCAAGTATTGAAGGATCATGATTAACGATATATTCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTaa                                                                                                              <  1:251229/140‑1 (MQ=255)                               tACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATg                                                                                                            <  1:67902/140‑1 (MQ=255)                                aCAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAAtaa                                                                                                           >  1:73492/1‑138 (MQ=255)                                aCAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGa                                                                                                           >  2:206899/1‑140 (MQ=255)                                      ttGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATa                                                                                                     >  1:1694/1‑140 (MQ=255)                                        gTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGa                                                                                                   <  1:248349/140‑1 (MQ=255)                                          cACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAAt                                                                                                 >  1:203529/1‑140 (MQ=255)                                           atcAAATATTGAAGTATCTTTGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATATCTTATGCTGTTTAAATACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAAtt                                                                                                <  1:20817/138‑1 (MQ=255)                                           aCCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATACTTGTAGCACTATTTCTTATGATGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAATGTTAAATTAATTAAAAAAAGAAtt                                                                                                >  2:50086/1‑140 (MQ=255)                                            ccAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTa                                                                                               <  1:28200/140‑1 (MQ=255)                                                  aTTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAat                                                                                         >  1:221237/1‑140 (MQ=255)                                                    tGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATAAGAGAatat                                                                                       >  1:81003/1‑140 (MQ=255)                                                           tCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGACAACTAGAATTAGAGAATATAGAAAAc                                                                                >  1:111553/1‑140 (MQ=255)                                                           tCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTAATATAATACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAAc                                                                                <  2:165008/140‑1 (MQ=255)                                                                 ttttCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAg                                                                          <  2:218256/140‑1 (MQ=255)                                                                          tatCTATAATTTTGCTTTATAATTGTAACACTATTTCTTATGTGGTTTAATTACTATACTAAAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTtctc                                                                 <  2:131114/140‑1 (MQ=255)                                                                          tatCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTtctc                                                                 >  1:207145/1‑140 (MQ=255)                                                                          tatCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAAGAATTtctc                                                                 >  1:153316/1‑140 (MQ=255)                                                                             ccAGAAATTTTCTTGATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAg                                                              <  1:16092/140‑1 (MQ=255)                                                                              cAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGa                                                             >  2:244082/1‑140 (MQ=255)                                                                              cAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGa                                                             <  1:20486/140‑1 (MQ=255)                                                                              cAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGa                                                             >  1:261928/1‑140 (MQ=255)                                                                              cAGAAATTTTCTTTATAATTGTAACAATATTTCTTATGCTGTTTAATTACTAAAATA‑AAAAAAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGa                                                             <  2:253276/140‑1 (MQ=255)                                                                                  aaTTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACt                                                         <  2:225942/140‑1 (MQ=255)                                                                                  aaTTTTCTTTATAATTGTAGCACTATTTATTATGCTGTTTAATTACTATACTAAAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACt                                                         <  2:238998/140‑1 (MQ=255)                                                                                   aTTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACtt                                                        >  2:38633/1‑140 (MQ=255)                                                                                       tCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAATAGACTTTAAg                                                    >  2:139569/1‑140 (MQ=255)                                                                                            ataATTGTAGCACTATTTCTTATGTTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCAAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAg                                               <  1:129954/140‑1 (MQ=255)                                                                                             taATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTaaaca                                              >  1:222166/1‑138 (MQ=255)                                                                                             taATTATAGCAAAATATAATATGCTGTTTAAATACTATACTACAAAT‑AAAAAATATACATCTAAAGGTAAAATTAATGAAAAATAGAATTAGAGAAAAAAAAAAAAAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGa                                              <  2:155894/140‑1 (MQ=255)                                                                                                  gTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAACATATTCATATAATGGTTTTATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACATTAAGTAAAGAATTaa                                         >  2:232193/1‑140 (MQ=255)                                                                                                     gCACTATTTCTTATGCTGTTTAATTACTATACAACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAAAAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAg                                      <  1:66557/140‑1 (MQ=255)                                                                                                      cACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAgt                                     >  1:4273/1‑140 (MQ=255)                                                                                                        cTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAgtgt                                   <  1:37485/140‑1 (MQ=255)                                                                                                        cTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAAAAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAgtgt                                   >  2:167567/1‑140 (MQ=255)                                                                                                          aTTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCt                                 <  1:12314/140‑1 (MQ=255)                                                                                                                tATGCTGTTTAATTACTATAATACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCaa                           <  2:111679/140‑1 (MQ=255)                                                                                                                 aTGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCaaa                          >  2:26154/1‑140 (MQ=255)                                                                                                                 aTGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCaaa                          >  1:233794/1‑140 (MQ=255)                                                                                                                  tGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAAc                         <  1:154744/140‑1 (MQ=255)                                                                                                                       tttAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTAAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTa                    >  2:9257/1‑140 (MQ=255)                                                                                                                          aaTTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATg                 <  1:38633/140‑1 (MQ=255)                                                                                                                          aaTTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATAAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATg                 >  2:8732/1‑140 (MQ=255)                                                                                                                           aTTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGCATTAAAAGTGACTAGGCAAACTATTAATGc                >  1:60378/1‑140 (MQ=255)                                                                                                                             tactatactaCAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCaa              >  1:196539/1‑140 (MQ=255)                                                                                                                               ctatactaCAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCAATa            >  1:138769/1‑140 (MQ=255)                                                                                                                                tatactaCAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAACATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAACGAATTAAAAGTGTCTAGGCAAACTATTAATGCAATAg           >  1:213702/1‑140 (MQ=255)                                                                                                                                tatactaCAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCAATAg           >  1:98927/1‑140 (MQ=255)                                                                                                                                       cAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCAATAGAAAACaa    <  2:121357/140‑1 (MQ=255)                                                                                                                                         aaTAAAAAAA‑AAAAAGCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCAATAGAAAACaata  <  2:168797/140‑1 (MQ=255)                                                                                                                                            |                                                                                                                                          TTGTATGTTATTTACACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCAATAGAAAACAATA  >  CP000731/25524‑25800

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