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	.	96.8%	265.0 / ‑0.2	62	coding (316/333 nt)	USA300HOU_pUSA300HOUMR0031	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base . (45/15);  total (45/17)
Fisher's exact test for biased strand distribution p-value = 7.19e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTTATTTACACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCAATAGAAAACAAT  >  CP000731/25529‑25799                                                                                                                                       |                                                                                                                                         tgttccttACAAAAGTGCTTAAATTTACAAGTTTGTCACCAAGTATTAAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaaat                                                                                                                                    <  1:56917/135‑2 (MQ=255) tGTTATTTACACTAGTGCTTAAATTTACAGGTCTGTCACGTTGTCCATAAAGATCTTGATTATAAAAGCATTCTGAAATTTTCTTTATAATTGTAGCACAATATCTTATGCTGTTTAATTACTATACTACAAATaaaaaat                                                                                                                                    >  1:208019/1‑140 (MQ=255)   ttATTTACACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaaatat                                                                                                                                  <  1:47257/141‑4 (MQ=255)          cacTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTa                                                                                                                          >  1:90940/1‑141 (MQ=255)          cacTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTa                                                                                                                          >  2:174592/1‑141 (MQ=255)            cTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTaaa                                                                                                                        >  1:29049/1‑141 (MQ=255)            cTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTaaa                                                                                                                        >  1:229135/1‑141 (MQ=255)             tAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTGGCGCTATTTCTTATGCTGTTTAATTACTATACTACAACT‑AAAAAATATACATCTAAAg                                                                                                                       >  2:71188/1‑141 (MQ=255)              aGTGCTTAAATTTACAGGTTTGTCACCAAGTAGTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAgg                                                                                                                      >  1:55052/1‑141 (MQ=255)               gTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGt                                                                                                                     <  2:191519/141‑1 (MQ=255)                      aaTTTACAGGTTTGTCACCAAGTATTGAAGGATCTTAGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTa                                                                                                              <  2:110220/141‑1 (MQ=255)                        tttACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAAt                                                                                                            >  1:29685/1‑141 (MQ=255)                           aCAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGaa                                                                                                         >  1:98295/1‑141 (MQ=255)                            cAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTGCTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGaaa                                                                                                        >  1:3845/1‑141 (MQ=255)                            cAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGaaa                                                                                                        >  1:151198/1‑141 (MQ=255)                            cAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGaaa                                                                                                        >  2:199069/1‑141 (MQ=255)                               gTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAAt                                                                                                     >  2:235737/1‑141 (MQ=255)                                    tCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAACTaaact                                                                                                >  2:204108/1‑137 (MQ=255)                                      aCCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAACTAGTATTa                                                                                              >  1:130873/1‑141 (MQ=255)                                       ccAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTag                                                                                             >  1:2490/1‑141 (MQ=255)                                        cAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAACTAGAATTaga                                                                                            >  1:125512/1‑141 (MQ=255)                                        cAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTaga                                                                                            >  2:117265/1‑141 (MQ=255)                                        cAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTaga                                                                                            >  1:165257/1‑141 (MQ=255)                                        cAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAAAATTaga                                                                                            >  1:222431/1‑141 (MQ=255)                                         aaGTATTGAAGGATCTTGGTTTTCTCAATATTCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGATGTTTAATTACTATACTAAAAAT‑AAAAAATATATATCTAAAGGTTCAATTAATGAAAAATAGAATTagag                                                                                           <  2:90961/141‑1 (MQ=255)                                            tATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAat                                                                                        >  2:249525/1‑141 (MQ=255)                                                gAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑CAAAAATATACATCTACAGGTTAACTTAATGAAAAATAGAATTAGAGAATATAg                                                                                    >  1:114556/1‑141 (MQ=255)                                                gAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAACTAGAATTAGAGAATATAg                                                                                    >  1:76816/1‑141 (MQ=255)                                                      tCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACCAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACa                                                                              >  2:57021/1‑141 (MQ=255)                                                                  aaaTATCAAGAAATTTTCTTTATAATTGTAGAAATATTTATTATGCTGTTTAATTACTATAAAACAAA‑AAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTtc                                                                  <  2:219933/140‑1 (MQ=255)                                                                    atatACAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTtctc                                                                <  2:240883/141‑1 (MQ=255)                                                                     tatCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCa                                                               >  2:43503/1‑141 (MQ=255)                                                                       tCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCacg                                                             >  1:97618/1‑139 (MQ=255)                                                                       tCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAg                                                             >  1:14109/1‑141 (MQ=255)                                                                        ccAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTCTACTACAAAT‑AAAAAATATACATCTAAAGATTACATTAATGACAAATAGAACTAGAGAATATAGAAAACAAAAAGGAATATCTCAAGa                                                            >  2:118691/1‑141 (MQ=255)                                                                        ccAGAAATTTTCTTTATAATTGTAGCACTATTTATTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATAAAAAATAGAATTAAAGAATATATAAAACAAATAGGAAACTCTAAAGa                                                            >  2:233488/1‑141 (MQ=255)                                                                           gAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGGAAACAAAAAGGAATTTCTCAAGAAAc                                                         >  1:90711/1‑141 (MQ=255)                                                                           gAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAAc                                                         >  2:228627/1‑141 (MQ=255)                                                                            aaaTTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACt                                                        <  1:144670/141‑1 (MQ=255)                                                                               ttttCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTAATATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTa                                                     <  2:13609/141‑1 (MQ=255)                                                                                   cTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTa                                                 >  2:29362/1‑141 (MQ=255)                                                                                   atttAAAATAGTAAAAATATTTCTTATGCTGTTTAAATACTATAATA‑AAAAACAAGAATAGACATCTAAAGGTTAAATGAAGGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTa                                                 <  2:229154/140‑1 (MQ=255)                                                                                    tttATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTaa                                                <  1:210099/141‑1 (MQ=255)                                                                                       ataATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGa                                             <  1:99501/141‑1 (MQ=255)                                                                                         aaTTGTAGCACTATTTCTTATGCTGTTTAATTACTATACCACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAAt                                           <  1:239658/141‑1 (MQ=255)                                                                                          aTTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAAtt                                          <  1:26740/141‑1 (MQ=255)                                                                                           ttatAGAACTATATATTATGCTGTTTAATTAATAAAGTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATAGAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTa                                         <  2:253204/138‑1 (MQ=255)                                                                                           ttGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTa                                         >  2:34034/1‑141 (MQ=255)                                                                                            tGTAGAAATATTTCTTATGCTATTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTaa                                        <  2:152129/141‑1 (MQ=255)                                                                                              tAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTaaaa                                      >  1:270791/1‑141 (MQ=255)                                                                                                  aCTATTGCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAgtgt                                  <  1:126932/141‑1 (MQ=255)                                                                                                      tttCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAg                              >  1:129579/1‑141 (MQ=255)                                                                                                       ttCTTATGCTGTTTAATTACAATAAAAAAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAgg                             <  2:278943/141‑1 (MQ=255)                                                                                                          ttATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCaa                          >  1:43051/1‑141 (MQ=255)                                                                                                          ttATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCaa                          >  1:183962/1‑141 (MQ=255)                                                                                                          ttATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCaa                          >  2:230327/1‑141 (MQ=255)                                                                                                          ttATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCaa                          >  2:224469/1‑141 (MQ=255)                                                                                                             tGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACt                       >  1:29080/1‑141 (MQ=255)                                                                                                                tGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTAtt                    <  2:2490/141‑1 (MQ=255)                                                                                                                  tttAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTaa                  >  1:131263/1‑141 (MQ=255)                                                                                                                   ttAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAAt                 >  1:96922/1‑141 (MQ=255)                                                                                                                     aaTTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATga               >  2:176738/1‑140 (MQ=255)                                                                                                                     aaTTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGCAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTATCTAAGCAAACTATTaaaga               >  2:71327/1‑138 (MQ=255)                                                                                                                              actacAAAT‑AAAAAATATACATCTAAAGATTTAATAAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAAAATTAAAAGTGTCTAGACAAACTATTAATGCAATATaaaa      >  2:230708/1‑141 (MQ=255)                                                                                                                              aacaaaaaaaaaaaaATATACATATAAAGGACAAATTAATGAAAACTAGAATTAGACAAAACACAAAACCAAACGCACCGTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTATAGGCAAACTATTAATGCAATAGaaa       <  2:155751/136‑1 (MQ=255)                                                                                                                                  cAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCAATAGAAAACaat  >  1:276748/1‑141 (MQ=255)                                                                                                                                       |                                                                                                                                         TGTTATTTACACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCAATAGAAAACAAT  >  CP000731/25529‑25799

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