Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	2,001,669	T→C	6.0%	noncoding (14/76 nt)	USA300HOU_t0051 ←	tRNA‑Ala

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	2,001,669	0	T	C	6.0%	201.7 / 3.8	67	noncoding (14/76 nt)	USA300HOU_t0051	tRNA‑Ala
Reads supporting (aligned to +/- strand):  ref base T (18/45);  new base C (2/2);  total (20/47)
Fisher's exact test for biased strand distribution p-value = 5.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.41e-01

ACTCTTTATTCACTCGGTTTTGCTTGGTAAAATCTATATTTTACTTACTTATCTAGTTTTCAATGTACAAATAATGGTGGAGACTAGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCTAAGTGGACTCGAACCACCGACCTCACGCTTATCAGGCGTGCGCTCTAACCAGCTGAGCTATAGGCCCATTAATTTGAATGAACAAACATTCAAAACTGA  >  CP000730/2001533‑2001809                                                                                                                                         |                                                                                                                                             aCTCTTTATTCACTCGGTTTTGCTTGGTAAGATCTATATTTTACTTACTTATCTAGTTTTCAATGTACAAATAATGGTGGAGACTAGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGc                                                                                                                                          <  2:52695/141‑1 (MQ=255) aCTCTTTATTCACTCGGTTTTGCTTGGTAAAATCTATATTTTACTTACTTATCTAGTTTTCAATGTACAAATAATGGTGGAGACTAGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGc                                                                                                                                          >  2:312471/1‑141 (MQ=255)   tctTTATTCACTCGGTTTTGCTTGGTAAAATCTATATTTTACTTACTTATCTAGTTTTCAATGTACAAATAATGGTGGAGACTAGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTa                                                                                                                                        >  2:483118/1‑141 (MQ=255)    ctTTATTCACTCGGTTTTGCTTGGTAAAATCTATATTTTACTTACTTATCTAGTTTTCAATGTACAAATAATGGTGGAGACTAGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTaa                                                                                                                                       >  1:131372/1‑141 (MQ=255)    ctTTATTCACTCGGTTTTGCTTGGTAAAATCTATATTTTACTTACTTATCTAGTTTTCAATGTACAAATAATGGTGGAGACTAGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTaa                                                                                                                                       >  2:613710/1‑141 (MQ=255)       tATTCACTCGGTTTTGCTTGGTAAAATCTATATTTTACTTACTTATCTAGTTTTCAATGTACAAATAATGGTGGAGACTAGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCCGAGCTAAGcc                                                                                                                                    >  1:708736/1‑141 (MQ=255)       tATTCACTCGGTTTTGCTTGGTAAAATCTATATTTTACTTACTTATCTAGTTTTCAATGTACAAATAATGGTGGAGACTAGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCCGAGCTAAGcc                                                                                                                                    >  1:708737/1‑141 (MQ=255)         ttCACTCGGTTTTGCTTGGTAAAATCTATATTTTACTTACTTATCTAGTTTTCAATGTACAAATAATGGTGGAGACTAGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGcccc                                                                                                                                  >  2:644603/1‑141 (MQ=255)                ggTTTTGCTTGGTAAAATTGCGAGTTTACTTACTTATCTAGTTTTCAATGTACAAATAATGGTGGAGACTAGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAg                                                                                                                           <  1:598893/141‑1 (MQ=21)                 gTTTTGCTTGGTAAAATCTATATTTTACTTACTTATCTAGTTTTCAATGTACAAATCAGGGTGGAGACTAGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAgg                                                                                                                          <  1:173509/141‑1 (MQ=255)                    ttGCTTGGTAAAATCTATATTTTACTTACTTATCTAGTTTTCAATGTACAAATAATGGTGGAGACTAGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTat                                                                                                                       <  2:631978/141‑1 (MQ=255)                     tGCTTGGTAAAATCTATATTTTACTTACTTATCTAGTTTTCAATGTACAAATAATGGTGGAGACTAGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTatt                                                                                                                      >  2:755628/1‑141 (MQ=255)                     tGCTTGGTAAAATCTATATTTTACTTACTTATCTAGTTTTCAATGTACAAATAATGGTGGAGACTAGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTatt                                                                                                                      <  2:230580/141‑1 (MQ=255)                      gCTTGGTAAAATCTATATTTTACTTACTTATCTAGTTTTCAATGTACAAATAATGGTGGAGACTAGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTatta                                                                                                                     >  2:784461/1‑141 (MQ=255)                        ttGGTAAAATCTATATTTTACTTACTTATCTAGTTTTCAATGTACAAATAATGGTGGAGACTAGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTattaaa                                                                                                                   <  2:192543/141‑1 (MQ=255)                         tGGTAAAATCTATATTTTACTTACTTATCTAGTTTTCAATGTACAAATAATGGTGGAGACTAGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTattaaat                                                                                                                  >  2:4923/1‑141 (MQ=255)                         tGGTAAAATCTATATTTTACTTACTTATCTAGTTTTCAATGTACAAATAATGGTGGAGACTAGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTattaaat                                                                                                                  <  1:96602/141‑1 (MQ=255)                               aaTCTATATTTTACTTACTTATCTAGTTTTCAATGTACAAATAATGGTGGAGACTAGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAAtgg                                                                                                            <  2:98597/141‑1 (MQ=255)                                 tCTATATTTTACTTACTTATCTAGTTTTCAATGTACTAATAATGGTGGAGACTAGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAAtggtg                                                                                                          <  2:332315/141‑1 (MQ=255)                                  cTATATTTTACTTACTTATCTAGTTTTCAATGTACAAATAATGGTGGAGACTAGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAAtggtgg                                                                                                         <  1:85411/141‑1 (MQ=255)                                    atatTTTACTTACTTATCTAGTTTTCAATGTACAAATAATGGTGGAGACTAGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGc                                                                                                       <  2:221899/141‑1 (MQ=255)                                         ttacttacTTATCTAGTTTTCAATGTACAAATAATGGTGGAGACTAGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCTAAg                                                                                                  >  2:129339/1‑141 (MQ=255)                                          tacttactTATCTAGTTTTCAATGTACAAATAATGGTGGAGACTAGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCTAAGt                                                                                                 <  2:669879/141‑1 (MQ=255)                                          tacttactTATCTAGTTTTCAATGTACAAATAATGGTGGAGACTAGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCTAAGt                                                                                                 <  2:23895/141‑1 (MQ=255)                                          tacttactTATCTAGTTTTCAATGTACAAATAATGGTGGAGACTAGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCTAAGt                                                                                                 <  1:445925/141‑1 (MQ=255)                                              tacttaTCTAGTTTTCAATGTACAAATAATGGTGGAGACTAGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCTAAGTGGAc                                                                                             <  2:206563/141‑1 (MQ=255)                                               acttaTCTAGTTTTCAATGTACAAATAATGGTGGAGACTAGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAAtggtggttcgagtccactt                                                                                            <  1:634001/141‑14 (MQ=255)                                               acttaTCTAGTTTTCAATGTACAAATAATGGTGGAGACTAGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCTTAAGTGGAcg                                                                                            <  2:717573/141‑2 (MQ=255)                                                 ttaTCTAGTTTTCAATGTACAAATAATGGTGGAGACTAGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCTAAGTGGACTCg                                                                                          <  2:268639/141‑1 (MQ=255)                                                                  tACAAATAATGGTGGAGACTAGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCTAAGTGGACTCGAACCACCGACCTCACGc                                                                         <  1:197430/141‑1 (MQ=255)                                                                  tACAAATAATGGTGGAGACTAGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCAAAGTGGACTCGAACCACCGACCTCACGc                                                                         <  2:778116/141‑1 (MQ=255)                                                                    cAAATAATGGTGGAGACTAGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCTAAGTGGACTCGAACCACCGACCTCACGCtt                                                                       <  1:710368/141‑1 (MQ=255)                                                                       ataatGGTGGAGACTAGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCTAAGTGGACTCGAACCACCGACCTCACGcctgat                                                                    <  2:34780/141‑6 (MQ=255)                                                                       ataatGGTGGAGACTAGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCTAAGTGGACTCGAACCACCGACCTCACGCTTATc                                                                    <  1:242470/141‑1 (MQ=255)                                                                       ataatGGTGGAGACTAGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCTAAGTGGACTCGAACCACCGACCTCACGCTTATc                                                                    >  1:188668/1‑141 (MQ=255)                                                                           tggtggAGACTAGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCTAAGTGGACTCGAACCACCGACCTCACGCTTATCAGGc                                                                <  1:706381/141‑1 (MQ=255)                                                                            ggtggAGACTAGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCTAAGTGGACTCGAACCACCGACCTCACGCTTATCAGGCg                                                               <  2:734131/141‑1 (MQ=255)                                                                             gtggAGACTAGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCTAAGTGGACTCGAACCACCGACCTCACGCTTATCAGGCGt                                                              <  2:455570/141‑1 (MQ=255)                                                                                 agaCTAGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCTAAGTGGACTCGAACCACCGACCTCACGCTTATCAGGCGTgcgc                                                          <  2:479533/141‑1 (MQ=255)                                                                                    cTAGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCTAAGTGGACTCGAACCACCGACCTCACGCTTATCCGGCGTGCGctct                                                       >  2:10555/1‑141 (MQ=255)                                                                                    cTAGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCTAAGTGGACTCGAACCACCGACCTCACGCTTATCAGGCGTGCGctct                                                       <  1:504627/141‑1 (MQ=255)                                                                                     tAGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCTAAGTGGACTCGAACCACCGACCTCACGCTTATCAGGCGTGCGCTCTa                                                      <  2:698566/141‑1 (MQ=255)                                                                                      aGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCTAAGTGGACTCGAACCACCGACCTCACGCTTATCAGGCGTGCGCTCTaa                                                     <  1:375816/141‑1 (MQ=255)                                                                                      aGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCTAAGTGGACTCGAACCACCGACCTCACGCTTATCAGGCGTGCGCTCTaa                                                     <  1:464805/141‑1 (MQ=255)                                                                                       gCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCTAAGTGGACTCGAACCACCGACCTCACGCTTATCAGGCGTGCGCTCTAAc                                                    <  2:383364/141‑1 (MQ=255)                                                                                               gAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCTAAGTGGACTCGAACCACCGACCTCACGCTTATCAGGCGTGCGCTCTAACCAGCTGAg                                            <  1:61632/141‑1 (MQ=255)                                                                                                aaCCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCTAAGTGGACTCGAACCACCGACCTCACGCTTCTCAGGCGTGCGCTCTAACCCGCTGAGc                                           >  1:285175/1‑141 (MQ=255)                                                                                                  ccGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCTAAGTGGACTCGAACCACCGACCTCACGCTTATCAGGCGTGCGCTCTAACCAGCTGAGCta                                         <  2:541263/141‑1 (MQ=255)                                                                                                    gCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCTAAGTGGACTCGAACCACCGACCTCACGCTTATCAGGCGTGCGCTCTAACCAGCTGAGCtata                                       >  2:236409/1‑141 (MQ=255)                                                                                                       gACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCTAAGTGGACTCGAACCACCGACCTCACGCTTATCAGGCGTGCGCTCTAACCAGCTGAGCTATAGGc                                    <  1:1022/141‑1 (MQ=255)                                                                                                        aCCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCTAAGTGGACTCGAACCACCGACCTCACGCTTATCAGGCGTGCGCTCTAACCAGCTGAGCTATAGGcc                                   <  2:560130/141‑1 (MQ=255)                                                                                                         acccctCCGTGCAAAGCAGGCGCTCTCCCAGGTGAGCTAAGCCCCCAAATAGTTATTAAATTAATGGAGGGCCGAAGTGGCCTCGAACCACCGACCGCACGCTTATCAGCCGTACGCGCTATCCACCTGTGCTATAGGccc                                  >  2:110459/4‑141 (MQ=25)                                                                                                           tcctGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCTAAGTGGACTCGAACCACCGACCTCACGCTTATCAGGCGTGCGCTCTAACCAGCTGAGCTATAGGCCCAt                                <  1:333963/141‑1 (MQ=255)                                                                                                               gCGTGCAAAGCCGGCGCTCTCCCAGCTGCGCTAAGCCCCCAAATAGGTATTACATTAATGGTGGGCCTAAGTGGACTCGAACCACCGACCCCACGCTTCGCAGCCGTGCGCTCTAACCAGCTGAGCTACAGGCCCATTGAt                            >  2:762384/1‑141 (MQ=21)                                                                                                                cGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCTAAGTGGACTCGAACCACCGACCTCACGCTTATCAGGCGTGCGCTCTAACCAGCTGAGCTATAGGCCCATTAAtt                           <  1:278390/141‑1 (MQ=255)                                                                                                                cGTGCAAAGCAGGCGCTCTCCCAGCCGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCTAAGTGGACTCGAACCACCGACCTCACGCTTATCAGGCGTGCGCTCTAACCAGCTGAGCTATAGGCCCATTAAtt                           <  2:708737/141‑1 (MQ=255)                                                                                                                cGTGCAAAGCAGGCGCTCTCCCAGCCGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCTAAGTGGACTCGAACCACCGACCTCACGCTTATCAGGCGTGCGCTCTAACCAGCTGAGCTATAGGCCCATTAAtt                           <  2:708736/141‑1 (MQ=255)                                                                                                                 gTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCTAAGTGGACTCGAACCACCGACCTCACGCTTATCAGGCGTGCGCTCTAACCAGCTGAGCTATAGGCCCATTAAttt                          <  1:715839/141‑1 (MQ=255)                                                                                                                  tGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCTAAGTGGACTCGAACCACCGACCTCACGCTTATCAGGCGTGCGCTCTAACCAGCTGAGCTATAGGCCCATTAATTtg                         >  1:446554/1‑141 (MQ=255)                                                                                                                    cAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCTAAGTGGACTCGAACCACCGACCTCACGCTTATCAGGCGTGCGCTCTAACCAGCTGAGCTATAGGCCCATTAATTtgaa                       >  1:66359/1‑141 (MQ=255)                                                                                                                      aaGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCTAAGTGGACTCGAACCACCGACCTCACGCTTATCAGGCGTGCGCTCTAACCAGCTGAGCTATAGGCCCATTAATTtgaatg                     <  1:613710/141‑1 (MQ=255)                                                                                                                        gCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCTAAGTGGACTCGAACCACCGACCTCACGCTTATCAGGCGTGCGCTCTAACCAGCTGAGCTATAGGCCCATTAACTtgaatgaa                   <  1:692521/141‑1 (MQ=255)                                                                                                                            gcgcTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCTAAGTGGACTCGAACCACCGACCTCACGCTTATCAGGCGTGCGCTCTAACCAGCTGAGCTATAGGCCCATTAATTTGAATGaacaaa               <  2:25539/141‑1 (MQ=255)                                                                                                                            gcgcTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCTAAGTGGACTCGAACCACCGACCTCACGCTTATCAGGCGTGCGCTCTAACCAGCTGAGCTATAGGCCCATTAATTTGAATGaacaaa               <  2:25534/141‑1 (MQ=255)                                                                                                                              gcTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCTAAGTGGACTCGAACCACCGACCTCACGCTTATCAGGCGTGCGCTCTAACCAGCTGAGCTATAGGCCCATTAATTTGAATGaacaaaca             >  2:140404/1‑141 (MQ=255)                                                                                                                               ctctCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCTAAGTGGACTCGAACCACCGACCTCACGCTTATCAGGCGTGCGCTCTAACCAGCTGAGCTATAGGCCCATTAATTTGAATGAACAAACAt            <  1:20824/141‑1 (MQ=255)                                                                                                                                   cccAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCTAAGTGGACTCGAACCACCGACCTCACGCTTATCAGGCGTGCGCTCTAACCAGCTGAGCTATAGGCCCATTAATTTGAATGAACAAACATTCaa        <  1:4923/141‑1 (MQ=255)                                                                                                                                         tGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCTAAGTGGACTCGAACCACCGACCTCACGCTTATCAGGCGTGCGCTCTAACCAGCTGAGCTATAGGCCCATTAATTTGAATGAACAAACATTCAAAACTGa  <  2:131372/141‑1 (MQ=255)                                                                                                                                         |                                                                                                                                             ACTCTTTATTCACTCGGTTTTGCTTGGTAAAATCTATATTTTACTTACTTATCTAGTTTTCAATGTACAAATAATGGTGGAGACTAGCGGGATCGAACCGCTGACCTCCTGCGTGCAAAGCAGGCGCTCTCCCAGCTGAGCTAAGCCCCCAAATAGGTATTAAATTAATGGTGGGCCTAAGTGGACTCGAACCACCGACCTCACGCTTATCAGGCGTGCGCTCTAACCAGCTGAGCTATAGGCCCATTAATTTGAATGAACAAACATTCAAAACTGA  >  CP000730/2001533‑2001809

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