Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,913,866	T→C	6.7%	Y213C (TAT→TGT)	ytbE ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,913,866	0	T	C	6.7%	181.8 / 3.3	60	Y213C (TAT→TGT)	ytbE	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (30/26);  new base C (2/2);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

TATTTTGATTATAATTTGAAATCGTTACATCGCCTTCATTTTTAAGATTAAAGTATAAAATTGAAGTTGGTGTATATTTGGCACCTAATTCTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGAATAATTTACAAATGATAATGAACGTTTGTCATCATTTTGATCCATCACAATAGTTTGCGGTCTAGATTTATCTAAATCCAATGTATCAAATACTTGTTCCATTGGTGGTAAATCTTTAAATTGACCGCCACTAATACC  >  CP000730/1913726‑1914002                                                                                                                                             |                                                                                                                                         tATTTTGATTATAATTTGAAATCGTTACATCGCCTTCATTTTTAAGATTAAAGTATAAAATTGAAGTTGGTGTATATTTGGCACCTAATTCTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGaat                                                                                                                                          <  2:241091/141‑1 (MQ=255)   ttttGATTATAATTTGAAATCGTTACATCGCCTTCATTTTTAAGATTAAAGTATAAAATTGAAGTTGGTGTATATTTGGCACCTAATTCTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGaataa                                                                                                                                        <  2:106617/141‑1 (MQ=255)   ttttGATTATAATTTGAAATCGTTACATCGCCTTCATTTTTAAGATTAAAGTATAAAATTGAAGTTGGTGTATATTTGGCACCTAATTCTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGaataa                                                                                                                                        >  1:418867/1‑141 (MQ=255)   ttttCATTATAATTTGAAACCGTTACATCGCCTTCATTTTTAAGACTAAAGAATAACATTGCAGTTGGTGTATATTTGGCACCTAATTCTTTTTGTAAGTGTTCAGCCAATTGTTTAATCGCCTCACTTTCATCTGaataa                                                                                                                                        >  2:491441/1‑141 (MQ=255)     ttGATTATAATTTGAAATCGTTACATCGCCTTCATTTTTAAGATTAAAGTATAAAATTGAAGTTGGTGTATATTTGGCACCTAATTCTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGAATAAtt                                                                                                                                      <  1:424973/141‑1 (MQ=255)       gATTATAATTTGAAATCGTTACATCGCCTTCATTTTTAAGATTAAAGTATAAAATTGAAGTTGGTGTATATTTGGCACCTAATTCTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGAATAATTTa                                                                                                                                    >  1:629953/1‑141 (MQ=255)       gATTATAATTTGAAATCGTTACATCGCCTTCATTTTTAAGATTAAAGTATAAAATTGAAGTTGGTGTATATTTGGCACCTAATTCTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGAATAATTTa                                                                                                                                    >  1:272802/1‑141 (MQ=255)       gATTATAATTTGAAATCGTTACATCGCCTTCATTTTTAAGATTAAAGTATAAAATTGAAGTTGGTGTATATTTGGCACCTAATTCTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGAACAATTTa                                                                                                                                    >  1:26496/1‑141 (MQ=255)       gATTATAATTTGAAATCGTTACATCGCCTTCATTTTTAAGATTAAAGTATAAAATTGAAGTTGGTGTATATTTGGCACCTAATTCTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGAACAATTTa                                                                                                                                    >  1:26493/1‑141 (MQ=255)        aTTATAATTTCAAATCACTACATAGACTTCATTTTTAAGATTAAAGTATAATATTGAAGTTGGTGTATCTTTGGCGCTTAATTCTTTTTCTAAGTTTTCAGCCAATTTTTTAATCGCCTCAATTTGATCTGCATAATTTAc                                                                                                                                   <  2:377905/141‑1 (MQ=255)          tataATTTGAAATCGTTACATCGCCTTCATTTTTAAGATTAAAGTATAAAATTGAAGTTGGTGTATATTTGGCACCTAATTCTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGAATAATTTACaa                                                                                                                                 <  1:214881/141‑1 (MQ=255)             aaTTTGAAATCGTTACATCGCCTTCATTTTTAAGATTAAAGTATAAAATTGAAGTTGGTGTATATTTGGCACCTAATTCTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGAATAATTTACAAATg                                                                                                                              <  2:696356/141‑1 (MQ=255)              aTTTGAAATCGTTACATCGCCTTCATTTTTAAGATTAAAGTATAAAATTGAAGTTGGTGTATATTTGGCACCTAATTCTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGAATAATTTACAAATGa                                                                                                                             >  1:643207/1‑141 (MQ=255)                  gAAATCGTTACATCGCCTTCATTTTTAAGATTAAAGTATAAAATTGAAGTTGGTGTATATTTGGCACCTAATTCTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGAATAATTTACAAATGATAAt                                                                                                                         <  2:562744/141‑1 (MQ=255)                   aaaTCGTTACATCGCCTTCATTTTTAAGATTAAAGTATAAAATTGAAGTTGGTGTATATTTGGCACCTAATTCTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGAATAATTTACAAATGATAATg                                                                                                                        <  2:622173/141‑1 (MQ=255)                    aaTCGTTACATCGCCTTCATTTTTAAGATTAAAGTATAAAATTGAAGTTGGTGTATATTTGGCACCTAATTCTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGAATAATTTACAAATGATAATGa                                                                                                                       >  2:115554/1‑141 (MQ=255)                     aTCGTTACATCGCCTTCATTTTTAAGATTAAAGTATAAAATTGAAGTTGGTGTATATTTGGCACCTAATTCTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGAATAATTTACAAATGATAATGaa                                                                                                                      >  2:429007/1‑141 (MQ=255)                      tCGTTACATCGCCTTCATTTTTAAGATTAAAGTATAAAATTGAAGTTGGTGTATATTTGGCACCTAATTCTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGAATAATTTACAAATGATAATGAAc                                                                                                                     >  1:113410/1‑141 (MQ=255)                        gTTACATCGCCTTCATTTTTAAGATTAAAGTATAAAATTGAAGTTGGTGTATATTTGGCACCTAATTCTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGAATAATTTACAAATGATAATGAACGt                                                                                                                   >  2:61758/1‑141 (MQ=255)                          tACATCGCCTTCATTTTTAAGATTAAAGTATAAAATTGAAGTTGGTGTATATTTGGCACCTAATTCTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCCATTTGATCTGAATAATTTACAAATGATAATGAACGttt                                                                                                                 <  2:351024/141‑1 (MQ=255)                            cATCGCCTTCATTTTTAAGATTAAAGTATAAAATTGAAGTTGGTGTATATTTGGCACCTAATTCTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGAATAATTTACAAATGATAATGAACGTTTGt                                                                                                               >  1:290220/1‑141 (MQ=255)                            cATCGCCTTCATTTTTAAGATTAAAGTATAAAATTGAAGTTGGTGTATATTTGGCACCTAATTCTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGAATAATTTACAAATGATAATGAACGTTTGt                                                                                                               >  2:161130/1‑141 (MQ=255)                                 ccTTCATTTTTAAGATTACACTATACCATCAACATGGGTGTATATCTGGCACCTAATTCTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGAATAATTTACAAATGATAATGAACGTTTGtcatca                                                                                                          <  1:670839/141‑1 (MQ=255)                                 ccTTCATTTTTAAGATTAAAGTATAAAATTGAAGTTGGTGTATATTTGGCACCTAATTCTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGAATAATTTACAAATGATAATGAACGTTTGtcatca                                                                                                          <  1:114465/141‑1 (MQ=255)                                  cTTCATTTTTAAGATTAAAGTATAAAATTGAAGTTGGTGTATATTTGGCACCTAATTCTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGAATAATTTACAAATGATAATGAACGTTTGTcatcat                                                                                                         >  1:546223/1‑141 (MQ=255)                                  cTTCATTTTTAAGATTAAAGTATAAAATTGAAGTTGGTGTATATTTGGCACCTAATTCTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGAATAATTTACAAATGATAATGAACGTTTGTcatcat                                                                                                         >  2:468978/1‑141 (MQ=255)                                     cATTTTTAAGATTAAAGTATAAAATTGAAGTTGGTGTATATTTGGCACCTAATTCTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGAATAATTTACAAATGATAATGAACGTTTGTCATCAtttt                                                                                                      >  1:185831/1‑141 (MQ=255)                                     cATTTTTAAGATTAAAGTATAAAATTGAAGTTGGTGTATATTTGGCACCTAATTCTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGAATAATTTACAAATGATAATGAACGTTTGTCATCAtttt                                                                                                      >  1:127261/1‑141 (MQ=255)                                     cATTTTTAAGATTAAAGTATAAAATTGAAGTTGGTGTATATTTGGCACCTAATTCTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGAATAATTTACAAATGATAATGAACGTTTGTCATCAtttt                                                                                                      <  1:94127/141‑1 (MQ=255)                                      aTTTTTAAGATTAAAGTATAAAATTGAAGTTGGTGTATATTTGGCACCTAATTCTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGAATAATTTACAAATGATAATGAACGTTTGTCATCATTTTg                                                                                                     >  2:466725/1‑141 (MQ=255)                                      aTTTTTAAGATTAAAGTATAAAATTGAAGTTGGTGTATATTTGGCACCTAATTCTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGAATAATTTACAAATGATAATGAACGTTTGTCATCATTTTg                                                                                                     >  2:433846/1‑141 (MQ=255)                                       tttttAAGATTAAAGTATAAAATTGAAGTTGGTGTATATTTGGCACCTAATTCTTTTTGTAAGTCTTCAGCCAATTGTTTACTCGCCTCACTTTGATCTGATTATTTTCACCATGAGAATGAACGTTTGTCATCATTTTgg                                                                                                    >  1:490868/1‑140 (MQ=255)                                       tttttAAGATTAAAGTATAAAATTGAAGTTGGTGTATATTTGGCACCTAATTCTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGaata                                                                                                                                         <  2:251913/104‑1 (MQ=255)                                       tttttAAGATTAAAGTATAAAATTGAAGTTGGTGTATATTTGGCACCTAATTCTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGaata                                                                                                                                         >  1:251913/1‑104 (MQ=255)                                                         aaaaTTGAAGTTGGTGTATATTTGGCACCTAATTCTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGAATAATTTACAAATGATAATGAACGTTTGTCATCATTTTGATCCATCACAATAGTTTGc                                                                                  <  1:17123/141‑1 (MQ=255)                                                         aaaaTTGAAGTTGGTGTATATTTGGCACCTAATTCTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGAATAATTTACAAATGATAATGAACGTTTGTCATCATTTTGATCCATCACAATAGTTTGc                                                                                  <  1:97793/141‑1 (MQ=255)                                                          aaaTTGAAGTTGGTGTATATTTGGCACCTAATTCTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGAATAATTTACAAATGATAATGAACGTTTGTCATCATTTTGATCCATCACAATAGTTTGCg                                                                                 >  1:263747/1‑141 (MQ=255)                                                            aTTGAAGTTGGTGTATATTTGGCACCTAATTCTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGAATAATTTACAAATGATAATGAACGTTTGTCATCATTTTGATCCATCACAATAGTTTGCGGt                                                                               <  1:171120/141‑1 (MQ=255)                                                             ttGAAGTTGGTGTATATTTGGCACCTAATTCTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGAATAATTTACAAATGATAATGAACGTTTGTCATCATTTTGATCCATCACAATAGTTTGCGGTc                                                                              <  1:450927/141‑1 (MQ=255)                                                                  gTTGGTGTATATTTGGCACCTAATTCTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGAATAATTTACAAATGATAATGAACGTTTGTCATCATTTTGATCCATCACAATAGTTTGCGGTCTAGAt                                                                         >  2:97810/1‑141 (MQ=255)                                                                    tGGTGTATATTTGGCACCTAATTCTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGAATAATTTACAAATGATAATGAACGTTTGTCATCATTTTGATCCATCACAATAGTTTGCGGTCTAGAttt                                                                       >  2:700033/1‑141 (MQ=255)                                                                      gtgtATATTTGGCACCTAATTCTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGAATAATTTACAAATGATAATGAACGTTTGTCATCATTTTGATCCATCACAATAGTTTGCGGTCTAGATTTAt                                                                     <  1:39698/141‑1 (MQ=255)                                                                                 gCACCTAATTCTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGAATAATTTACAAATGATAATGAACGTTTGTCATCATTTTGATCCATCACAATAGTTTGCGGTCTAGATTTATCTAAATCCAAt                                                          >  1:142239/1‑141 (MQ=255)                                                                                    ccTAATTCTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGAATAATTTACAAATGATAATGAACGTTTGTCATCATTTTGATCCATCACAATAGTTTGCGGTCTAGATTTATCTAAATCCAATGTa                                                       <  1:50109/141‑1 (MQ=255)                                                                                           cTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGAACAATTTACAAATGATAATGAACGTTTGTCATCATTTTGATCCATCACAATAGTTTGCGGTCTAGATTTATCTAAATCCAATGTATCAAATa                                                <  2:26496/141‑1 (MQ=255)                                                                                           cTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGAACAATTTACAAATGATAATGAACGTTTGTCATCATTTTGATCCATCACAATAGTTTGCGGTCTAGATTTATCTAAATCCAATGTATCAAATa                                                <  2:26493/141‑1 (MQ=255)                                                                                            tttttGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGAATAATTTACAAATGATAATGAACGTTTGTCATCATTTTGATCCATCACAATAGTTTGCGGTCTAGATTTATCTAAATCCAATGTATCAAATAc                                               >  2:288112/1‑141 (MQ=255)                                                                                            tttttGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGAATAATTTACAAATGATAATGAACGTTTGTCATCATTTTGATCCATCACAATAGTTTGCGGTCTAGATTTATCTAAATCCAATGTATCAAATAc                                               >  1:450836/1‑141 (MQ=255)                                                                                              tttGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGAATAATTTACAAATGATAATGAACGTTTGTCATCATTTTGATCCATCACAATAGTTTGCGGTCTAGATTTATCTAAATCCAATGTATCAAATACtt                                             >  1:612983/1‑141 (MQ=255)                                                                                                tGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGAATAATTTACAAATGATAATGAACGTTTGTCATCATTTTGATCCATCACAATAGTTTGCGGTCTAGATTTATCTAAATCCAATGTATCAAATACTTGt                                           <  1:487594/141‑1 (MQ=255)                                                                                                tGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGAATAATTTACAAATGATAATGAACGTTTGTCATCATTTTGATCCATCACAATAGTTTGCGGTCTAGATTTATCTAAATCCAATGTATCAAATACTTGt                                           <  1:419190/141‑1 (MQ=255)                                                                                                       ttttaGCCAATTGTTTAATCGCCTCAATTTGATCTGAATAATTTACAAATGATAATGAGCGTTTTTCATCATTTTTATCCATCACAATAGTTTGAGGTCTAGATATATCTAAATCCAATGTATCAAATACTTGTTCCATtg                                    >  2:279711/5‑141 (MQ=255)                                                                                                             ccAATTGTTTAATCGCCTCAATTTGATCTGAATAATTTACAAATGATAATGAACGTTTGTCATCATTTTGATCCATCACAATAGTTTGCGGTCTAGATTTATCTAAATCCAATGTATCAAATACTTGTTCCATTGGTGGTa                              <  2:629953/141‑1 (MQ=255)                                                                                                             ccAATTGTTTAATCGCCTCAATTTGATCTGAATAATTTACAAATGATAATGAACGTTTGTCATCATTTTGATCCATCACAATAGTTTGCGGTCTAGATTTATCTAAATCCAATGTATCAAATACTTGTTCCATTGGTGGTa                              >  1:145661/1‑141 (MQ=255)                                                                                                              cAATTGTTTAATCGCCTCAATTTGATCTGAATAATTTACAAATGATAATGAACGTTTGTCATCATTTTGATCCATCACAATAGTTTGCGGTCTAGATTTATCTAAATCCAATGTATCAAATACTTGTTCCATTGGTGGTaa                             <  2:113410/141‑1 (MQ=255)                                                                                                              cAATTGTTTAATCGCCTCAATTTGATCTGAATAATTTACAAATGATAATGAACGTTTGTCATCATTTTGATCCATCACAATAGTTTGCGGTCTAGATTTATCTAAATCCAATGTATCAAATACTTGTTCCATTGGTGGTaa                             >  1:248297/1‑141 (MQ=255)                                                                                                                   gTTTAATCGCCTCAATTTGATCTGAATAATTTACAAATGATAATGAACGTTTGTCATCATTTTGATCCATCACAATAGTTTGCGGTCTAGATTTATCTAAATCCAATGTATCAAATACTTGTTCCATTGGTGGTAAATCtt                        <  2:341718/141‑1 (MQ=255)                                                                                                                   gTTTAATCGCCTCAATTTGATCTGAATAATTTACAAATGATAATGAACGTTTGTCATCATTTTGATCCATCACAATAGTTTGCGGTCTAGATTTATCTAAATCCAATGTATCAAATACTTGTTCCATTGGTGGTAAATCtt                        <  1:498864/141‑1 (MQ=255)                                                                                                                      tAATCGCCTCAATTTGATCTGAATAATTTACAAATGATAATGAACGTTTGTCATCATTTTGATCCATCACAATAGTTTGCGGTCTAGATTTATCTAAATCCAATGTATCAAATACTTGTTCCATTGGTGGTAAATCTTTaa                     >  2:123554/1‑141 (MQ=255)                                                                                                                            ccTCAATTTGATCTGAATAATTTACAAATGATAATGAACGTTTGTCATCATTTTGATCCATCACAATAGTTTGCGGTCTAGATTTATCTAAATCCAATGTATCAAATACTTGTTCCATTGGTGGTAAATCTTTAAATTGAc               <  2:106648/141‑1 (MQ=255)                                                                                                                              tCAATTTGATCTGAATAATTTACAAATGATAATGAACGTTTGTCATCATTTTGATCCATCACAATAGTTTGCGGTCTAGATTTATCTAAATCCAATGTATCAAATACTTGTTCCATTGGTGGTAAATCTTTAAATTGACCg             <  1:492828/141‑1 (MQ=255)                                                                                                                                aaTTTGATCTGAATAATTTACAAATGATAATGAACGTTTGTCATCATTTTGATCCATCACAATAGTTTGCGGTCTAGATTTATCTAAATCCAATGTATCAAATACTTGTTCCATTGGTGGTAAATCTTTAAATTGACCGcc           >  2:305159/1‑141 (MQ=255)                                                                                                                                 aTTTGATCTGAATAATTTACAAATGATAATGAACGTTTGTCATCATTCTGATCCATCACAATAGTTTGCGGTCTAGATTTATCTAAATCCAATGTATCAAATACTTGTTCCATTGGTGGTAAATCTTTAAATTGACCGCCa          >  1:458098/1‑141 (MQ=255)                                                                                                                                   ttGATCTGAATAATTTACAAATGATAATGAACGTTTGTCATCATTTTGATCCATCACAATAGTTTGCGGTCTAGATTTATCTAAATCCAATGTATCAAATACTTGTTCCATTGGTGGTAAATCTTTAAATTGACCGCCACt        <  2:234214/141‑1 (MQ=255)                                                                                                                                   ttGATCTGAATAATTTACAAATGATAATGAACGTTTGTCATCATTTTGATCCATCACAATAGTTTGCGGTCTAGATTTATCTAAATCCAATGTATCAAATACTTGTTCCATTGGTGGTAAATCTTTAAATTGACCGCCACt        >  1:144419/1‑141 (MQ=255)                                                                                                                                         tGAATAATTTACAAATGATAATGAACGTTTGTCATCATTTTGATCCATCACAATAGTTTGCGGTCTAGATTTATCTAAATCCAATGTATCAAATACTTGTTCCATTGGTGGTAAATCTTTAAATTGACCGCCACTAATAcc  <  2:546223/141‑1 (MQ=255)                                                                                                                                             |                                                                                                                                         TATTTTGATTATAATTTGAAATCGTTACATCGCCTTCATTTTTAAGATTAAAGTATAAAATTGAAGTTGGTGTATATTTGGCACCTAATTCTTTTTGTAAGTCTTCAGCCAATTGTTTAATCGCCTCAATTTGATCTGAATAATTTACAAATGATAATGAACGTTTGTCATCATTTTGATCCATCACAATAGTTTGCGGTCTAGATTTATCTAAATCCAATGTATCAAATACTTGTTCCATTGGTGGTAAATCTTTAAATTGACCGCCACTAATACC  >  CP000730/1913726‑1914002

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