Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	2,772,254	T→G	5.8%	S351S (TCT→TCG)	phoB →	alkaline phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	2,772,254	0	T	G	5.8%	240.3 / 8.9	87	S351S (TCT→TCG)	phoB	alkaline phosphatase
Reads supporting (aligned to +/- strand):  ref base T (47/34);  new base G (2/3);  total (50/37)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.25e-01

TCGATGACGCCATTCAATATGCTAAGAAACATAAAGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAGATGATATGAAAAAGGTTAAAAAAGAGAATAAAAAGCTAA  >  CP000730/2772124‑2772390                                                                                                                                   |                                                                                                                                         tCGATGACGCCATTCAATATGCTAAGAAACATAAAGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACTAATGGAATCCTCAACCGATTAAATCTATGAAACAct                                                                                                                                >  2:1084983/1‑141 (MQ=255) tCGATGACGCCATTCAATATGCTAAGAAACATAAAGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGAAAAGATAAAGGCTACGAATGAAATCCTAAACCGATTAAATCAAAGAAACAct                                                                                                                                >  1:209495/1‑141 (MQ=255) tCGATGACGCCATTCAAGATGCTAAGAAACATAAAGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACAct                                                                                                                                <  1:604440/141‑1 (MQ=255)  cGATGACGCCATTCAATATGCTAAGAAACATAAAGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAATGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACActc                                                                                                                               >  2:274408/1‑141 (MQ=255)  cGATGACGCCATTCAATATGCTAAGAAACATAAAGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAAAAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAGAGCTATGAAACActc                                                                                                                               >  1:227222/1‑141 (MQ=255)   gATGACGCCATTCAATATTCTAAGAAACATAAAGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACActct                                                                                                                              <  2:364363/141‑1 (MQ=255)   gATGACGCCATTCAATATGCTAAGAAACATAAAGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACActct                                                                                                                              <  2:531000/141‑1 (MQ=255)     tGACGCCATTCAATATGCTAAGAAACATAAAGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTgg                                                                                                                            <  2:31606/141‑1 (MQ=255)     tGACGCCATTCAATATGCTAAGAAACATAAAGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTgg                                                                                                                            <  1:726608/141‑1 (MQ=255)       aCGCCATTCAATATGCTAAGAAACATAAAGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGAt                                                                                                                          >  2:134972/1‑141 (MQ=255)        cGCCATTCAATATGCTAAGAAACATAAAGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATc                                                                                                                         >  2:546371/1‑141 (MQ=255)        cGCCATTCAATATGCTAAGAAACATAAAGATACACGTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGATACACTCTGGATc                                                                                                                         <  2:116967/141‑1 (MQ=255)         gCCATTCAATATGCTAAGAAACATAAAGATACACTTTTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCa                                                                                                                        <  1:341958/141‑1 (MQ=255)         gCCATTCAATATGCTAAGAAACATAAAGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCa                                                                                                                        <  2:444084/141‑1 (MQ=255)         gCCATTCAATATGCTAAGAAACATAAAGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCGATGAAACACTCTGGATCa                                                                                                                        <  1:410429/141‑1 (MQ=255)         gCCATTCAATATGCTAAGAAACATAAAGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCGATGAAACACTCTGGATCa                                                                                                                        <  1:1053644/141‑1 (MQ=255)              tCAATATGCTAAGAAACATAAAGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCatatat                                                                                                                   <  2:1080039/141‑1 (MQ=255)                     gCTAAGAAACATAAAGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGaa                                                                                                            >  2:262231/1‑141 (MQ=255)                         aGAAACATAAAGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGaaaaaa                                                                                                        >  1:889253/1‑141 (MQ=255)                          gAAACATAAAGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAAt                                                                                                       >  1:356679/1‑141 (MQ=255)                            aaCATAAAGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATgg                                                                                                     <  1:160428/141‑1 (MQ=255)                                 aaaGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTaaa                                                                                                <  1:864083/141‑1 (MQ=255)                                  aaGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAATGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATACATGACAGAAAAAATGGCTAAAg                                                                                               >  2:22367/1‑141 (MQ=255)                                   aGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAgg                                                                                              >  2:57041/1‑141 (MQ=255)                                   aGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAgg                                                                                              >  2:883223/1‑141 (MQ=255)                                      tACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTGTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGa                                                                                           >  2:934716/1‑141 (MQ=255)                                         acTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGa                                                                                        >  2:511993/1‑141 (MQ=255)                                          cTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGAc                                                                                       <  2:65860/141‑1 (MQ=255)                                            tGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGAccc                                                                                     <  1:211690/141‑1 (MQ=255)                                             gTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCt                                                                                    <  1:134972/141‑1 (MQ=255)                                               cGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGa                                                                                  <  2:113939/141‑1 (MQ=255)                                                 ttGCAACAGCAGATCACTCTACTGGTGGTCTAACAATGGGTAAAGATAAAGGCTACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGaaa                                                                                >  1:23110/1‑141 (MQ=255)                                                  tGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGaaaa                                                                               >  2:410680/1‑141 (MQ=255)                                                  tGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGaaaa                                                                               >  1:245484/1‑141 (MQ=255)                                                  tGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGaaaa                                                                               >  1:271064/1‑141 (MQ=255)                                                  tGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGaaaa                                                                               <  2:971337/141‑1 (MQ=255)                                                   gCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGaaaaa                                                                              >  2:1024897/1‑141 (MQ=255)                                                    cAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAg                                                                             <  1:380754/141‑1 (MQ=255)                                                     aaCAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGt                                                                            <  2:566682/141‑1 (MQ=255)                                                      aCAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGtt                                                                           >  2:107883/1‑141 (MQ=255)                                                      aCAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATAGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGAAAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGtt                                                                           <  2:980633/141‑1 (MQ=255)                                                              tCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGa                                                                   >  1:1027791/1‑141 (MQ=255)                                                              tCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGa                                                                   >  1:91724/1‑141 (MQ=255)                                                              tCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGa                                                                   >  2:520232/1‑141 (MQ=255)                                                               cACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGaa                                                                  >  2:937563/1‑141 (MQ=255)                                                                 ctctACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAgg                                                                >  2:682370/1‑141 (MQ=255)                                                                  tctACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCGATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGt                                                               <  1:870956/141‑1 (MQ=255)                                                                    tACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTa                                                             <  1:692556/141‑1 (MQ=255)                                                                      cTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCGATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATg                                                           >  2:871346/1‑141 (MQ=255)                                                                       tggaggTGTAACAATAGCAAAAGATAAAGAAGACGAAGGGAATCCTCAACCGATTAAATCTATGAAACACTCTGAATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATgg                                                          <  2:378186/141‑1 (MQ=255)                                                                        ggtggtCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCAGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGTAGGTTATGGa                                                         <  2:895075/141‑1 (MQ=255)                                                                          tggtCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGAtt                                                       <  1:160803/141‑1 (MQ=255)                                                                           ggtCTAACAATTGGTAAAGATAACGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTACAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTCTGGAttt                                                      >  1:765685/1‑141 (MQ=255)                                                                           ggtCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGAttt                                                      >  2:139021/1‑141 (MQ=255)                                                                           ggtCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTAAGAAACACTCGGGACCATATAGGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGAttc                                                      >  1:647773/1‑140 (MQ=255)                                                                             tCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGa                                                    <  1:581345/141‑1 (MQ=255)                                                                                aaCAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAAtt                                                 >  2:822839/1‑141 (MQ=255)                                                                                aaCAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAAtt                                                 >  1:705723/1‑141 (MQ=255)                                                                                  cAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAAtttt                                               >  1:218742/1‑141 (MQ=255)                                                                                  cAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCGATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGCTTTGAAtttt                                               >  2:119942/1‑141 (MQ=255)                                                                                      tGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTcaca                                           >  1:147080/1‑141 (MQ=255)                                                                                        gTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACACTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAg                                         >  2:436124/1‑141 (MQ=255)                                                                                            aGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAgatga                                     <  1:100944/141‑1 (MQ=255)                                                                                            aGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAgatga                                     <  2:808396/141‑1 (MQ=255)                                                                                            aGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAgatga                                     <  1:1435/141‑1 (MQ=255)                                                                                            aGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAgatga                                     <  1:713554/141‑1 (MQ=255)                                                                                             gATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAgatgat                                    >  1:226416/1‑141 (MQ=255)                                                                                                 aaGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAGATGATATGa                                >  1:115619/1‑141 (MQ=255)                                                                                                   ggATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGTACTCGCGGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACGAGATGATATGaaa                              >  1:19797/1‑141 (MQ=255)                                                                                                    gATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAGATGATATGaaaa                             >  1:582889/1‑141 (MQ=255)                                                                                                    gATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAGATGATATGaaaa                             >  2:150802/1‑141 (MQ=255)                                                                                                      tACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAGATGATATGAAAAAg                           >  1:532939/1‑141 (MQ=255)                                                                                                        cGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAGATGATATGAAAAAGGt                         <  2:80760/141‑1 (MQ=255)                                                                                                        cGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAGATGATATGAAAAAGGt                         >  2:1030650/1‑141 (MQ=255)                                                                                                         gAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAGATGATATGAAAAAGGtt                        >  1:164992/1‑141 (MQ=255)                                                                                                             ggAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAGATGATATGAAAAAGGTTaaaa                    >  1:696970/1‑141 (MQ=255)                                                                                                                aTCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCAAAAGGTGAAGACCCTGAAAAAGATATTAATGAAGGTTATGGATTTGAATTTTCACAAGATGATATGAAAAAGGTTAAAAAag                 >  2:607996/1‑141 (MQ=255)                                                                                                                aTCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATCACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAGATGATATGAAAAAGGTTAAAAAag                 >  1:572806/1‑141 (MQ=255)                                                                                                                  ccTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAGATGATATGAAAAAGGTTAAAAAagag               <  1:678967/141‑1 (MQ=255)                                                                                                                       aCCGATTAAATCTATGAAACACTCTGGATAATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAGATGATATGAAAAAGGTTAAAAAAGAGAATaa          <  2:824889/141‑1 (MQ=255)                                                                                                                            ttAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAGATGATATGAAAAAGGTTAAAAAAGAGAATAAAAAGc     <  1:57041/141‑1 (MQ=255)                                                                                                                            ttAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAGATGATATGAAAAAGGTTAAAAAAGAGAATAAAAAGc     >  2:802091/1‑141 (MQ=255)                                                                                                                              aaaTCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAGATGATATGAAAAAGGTTAAAAAAGAGAATAAAAAGCTa   >  1:994129/1‑141 (MQ=255)                                                                                                                               aaTCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAGATGATATGAAAAAGGTTAAAAAAGAGAATAAAAAGCTaa  <  2:889253/141‑1 (MQ=255)                                                                                                                               aaTCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAGATGATATGAAAAAGGTTAAAAAAGAGAATAAAAAGCTaa  <  2:387531/141‑1 (MQ=255)                                                                                                                               aaTCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAGATGATATGAAAAAGGTTAAAAAAGAGAATAAAAAGCTaa  >  1:492218/1‑141 (MQ=255)                                                                                                                               aaTCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAACGACACTTATTTGAATTTTCACAAGATGATATGAAAAAGGTTAAAAAAGAGAATAAAAAGCTaa  >  1:565585/1‑141 (MQ=255)                                                                                                                                   |                                                                                                                                         TCGATGACGCCATTCAATATGCTAAGAAACATAAAGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAGATGATATGAAAAAGGTTAAAAAAGAGAATAAAAAGCTAA  >  CP000730/2772124‑2772390

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