Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	463,014	C→A	E47* (GAA→TAA)	USA300HOU_0437 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	463,014	0	C	A	94.6%	237.1 / 7.3	74	E47* (GAA→TAA)	USA300HOU_0437	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (1/3);  new base A (39/31);  total (40/34)
Fisher's exact test for biased strand distribution p-value = 3.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

GCCGTCTTGTTTAGCTTTGTTGATAATGTCTTTGATTGAATGTAGTCCTTTATCGGCGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTCTAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGCTAAAACACCTACCATTGCCGTAGCTACTAAAACCTTTTTAAAT  >  CP000730/462877‑463147                                                                                                                                          |                                                                                                                                      gCCGTCTTGTTTAGCTTTGTTGATAATGTCTTTGATTGAATGTAGTCCTTTATCGGCGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATaa                                                                                                                                    <  1:679040/141‑1 (MQ=255)  ccGTCTTGCTTAGCTTTGTTGATAATGTCTTTGATTGAATGTAGTCCTTTATCGGCGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCGATATAGCTTAATTCATGTTCTaag                                                                                                                                   <  2:76840/141‑1 (MQ=255)   cGTCTTGTTTAGCTTTGTTGATAATGTCTTTGATTGAATGTAGTCCTTTATCGGCGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATCCTTATCTAATAACCCATCAATATAACTTAATTCATGTTATaaga                                                                                                                                  >  2:35113/1‑141 (MQ=255)     tCTTGTTTAGCTTTGTTGATAATGTCTTTGATTGAATGTAGTCCTTTATCGGCGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATaagaag                                                                                                                                >  1:290888/1‑141 (MQ=255)       ttGTTTAGCTTTGTTGATAATGTCTTTGATTGAATGTAGTCCTTTATAGGCGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGtt                                                                                                                              <  1:450519/141‑1 (MQ=255)          tttAGCTTTGTTGATAATGTCTTTGATTGAATGTAGTCCTTTATCGGCGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGc                                                                                                                           <  1:235620/141‑1 (MQ=255)               cTTTGTTGATAATGTCTTTGATTGAATGTAGTCCTTTATCGGCGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGcatcat                                                                                                                      <  2:806765/141‑1 (MQ=255)                tttGTTGATAATGTCTTTGATTGAATGTAGTCCTTTATCGGCGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGcatcatc                                                                                                                     <  2:357631/141‑1 (MQ=255)                 ttgttgATAATGTCTTTGATTGAATGTAGTCCTTTATCGGCGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGcatcatca                                                                                                                    >  1:508109/1‑141 (MQ=255)                  tgttgATAATGTCTTTGATTGAATGTAGTCCTTTATCGGCGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGcatcatcat                                                                                                                   <  2:630812/141‑1 (MQ=255)                   gttgATAATGTCTTTGATTGAATGTAGTCCTTTATCGGCGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATg                                                                                                                  >  1:506598/1‑141 (MQ=255)                             tCTTTGATTGAATGTAGTCCTTTATCGGCGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGt                                                                                                        >  2:834614/1‑141 (MQ=255)                             tCTTTGATTGAATGTAGTCCTTTATCGGCGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGt                                                                                                        >  2:189551/1‑141 (MQ=255)                             tCTTTGATTGAATGTAGTCCTTTATCGGCGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGt                                                                                                        <  1:155747/141‑1 (MQ=255)                                ttgattgaATGTAGTCCTTTATCGGCGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCa                                                                                                     <  2:719658/141‑1 (MQ=255)                                   attgaATGTAGTCCTTTATCGGCGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCAttt                                                                                                  >  1:299066/1‑141 (MQ=255)                                     tgaATGTAGTCCTTTATCGGCGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTg                                                                                                <  1:555521/141‑1 (MQ=255)                                        aTGTAGTCCTTTATCGGCGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTgaga                                                                                             <  2:377894/141‑1 (MQ=255)                                         tGTAGTCCTTTATCGGCGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGaa                                                                                            <  1:526623/141‑1 (MQ=255)                                           tAGTCCATTATCGGCGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCGTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATa                                                                                          >  2:566062/1‑141 (MQ=255)                                               ccTTTATCGGCGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACt                                                                                      >  1:247748/1‑141 (MQ=255)                                                 tttATCGGCGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATATCTTAATTCATGTTATAAGAAGTTTGCATCGTCATGTGGTACGGGTCCATTTTGAGAATAGACttt                                                                                    <  2:75962/141‑1 (MQ=255)                                                 tttATCGGCGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACttt                                                                                    >  2:205378/1‑141 (MQ=255)                                                  ttATCGGCGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTc                                                                                   >  1:195133/1‑141 (MQ=255)                                                     tCGGCGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCa                                                                                >  2:514956/1‑141 (MQ=255)                                                       ggCGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATc                                                                              >  1:55922/1‑141 (MQ=255)                                                        gCGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCt                                                                             <  2:23206/141‑1 (MQ=255)                                                         cGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTCTAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTg                                                                            <  1:822376/141‑1 (MQ=255)                                                         cGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTg                                                                            >  2:786941/1‑141 (MQ=255)                                                         cGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTg                                                                            <  1:179736/141‑1 (MQ=255)                                                         cGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTg                                                                            >  2:355748/1‑141 (MQ=255)                                                           aaGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCt                                                                          >  1:817190/1‑141 (MQ=255)                                                              tATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGa                                                                       >  1:348588/1‑141 (MQ=255)                                                                tGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGAtt                                                                     >  2:202823/1‑141 (MQ=255)                                                                 gATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTa                                                                    >  2:590917/1‑141 (MQ=255)                                                                  aTCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTAc                                                                   <  1:279152/141‑1 (MQ=255)                                                                  aTCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTAc                                                                   <  1:221640/141‑1 (MQ=255)                                                                       aaGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTCTGAGAATAAACTTTCGCATCTGCTTGATTACCatat                                                              >  1:627669/1‑141 (MQ=255)                                                                        aGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATc                                                             >  1:408231/1‑141 (MQ=255)                                                                        aGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATc                                                             <  2:671125/141‑1 (MQ=255)                                                                          ttGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCca                                                           <  2:392971/141‑1 (MQ=255)                                                                           tGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCCTGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCcaa                                                          <  1:690492/141‑1 (MQ=255)                                                                            gTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCcaac                                                         >  2:388906/1‑141 (MQ=255)                                                                             tCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCcaaca                                                        <  1:653742/141‑1 (MQ=255)                                                                             tCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCcaaca                                                        <  2:813090/141‑1 (MQ=255)                                                                               ttttGTAGCTTGGTCAGCATTCTTATCTAATAACACATCCATATAACTTAATTCATGTTCTAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCaacaac                                                      <  1:448383/141‑1 (MQ=255)                                                                                 ttGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGc                                                    >  2:624377/1‑141 (MQ=255)                                                                                  tGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGcc                                                   <  2:198003/141‑1 (MQ=255)                                                                                    tAGCTTGGTCAGCATTCTTATCTAATACCGCGTCACTTTAACTTTATTCATGTTATAAGAAGTTTGCATCATCTTGTAGTACGAGTCCATTTTGAGAAGAAACTTTCGTATCTGCTTGATTACCATATCCAACAACGCCAg                                                 <  2:127387/141‑1 (MQ=255)                                                                                        ttGGTCAGCATTCTTATCTAATAACACATCAATCTAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGc                                             >  2:139379/1‑141 (MQ=255)                                                                                          ggTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGCTa                                           <  2:741201/141‑1 (MQ=255)                                                                                           gTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGCTaa                                          <  1:482495/141‑1 (MQ=255)                                                                                                cATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGCTAAAacac                                     >  1:234281/1‑141 (MQ=255)                                                                                                    cTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGCTAAAACACCTAc                                 >  1:158906/1‑141 (MQ=255)                                                                                                    cTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGCTAAAACACCTAc                                 >  1:650176/1‑141 (MQ=255)                                                                                                     ttATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGCTAAAACACCTAcc                                >  2:113603/1‑141 (MQ=255)                                                                                                     ttATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGCTAAAACACCTAcc                                <  1:233372/141‑1 (MQ=255)                                                                                                     ttATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCAACATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGCTAAAACACCTAcc                                >  2:565455/1‑141 (MQ=255)                                                                                                       aTCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGCTAAAACACCTACCAt                              >  1:44961/1‑141 (MQ=255)                                                                                                       aTCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGCTAAAACACCTACCAt                              >  1:405695/1‑141 (MQ=255)                                                                                                           aataaCACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGCTAAAACACCTACCATTGcc                          <  2:508109/141‑1 (MQ=255)                                                                                                            ataaCACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGCTAAAACACCTACCATTGCCg                         >  1:480177/1‑141 (MQ=255)                                                                                                             taaCACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGCTAAAACACCTACCATTGCCGt                        <  1:697468/141‑1 (MQ=255)                                                                                                              aaCACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGCTAAAACACCTACCATTGCCGTa                       >  1:338437/1‑141 (MQ=255)                                                                                                                      aaTATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGCTAAAACACCTACCATTGCCGTAGCTACTaa               >  1:569421/1‑141 (MQ=255)                                                                                                                       atatAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGCTAAAACACCTACCATTGCCGTAGCTACTaaa              >  2:317058/1‑141 (MQ=255)                                                                                                                            aCTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGCTAAAACACCTACCATTGCCGTAGCTACTAAAACCtt         >  2:171411/1‑141 (MQ=255)                                                                                                                             cTTAATTCATGTTCTAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGCTAAAACACCTACCATTGCCGTAGCTACTAAAACCttt        >  2:111894/1‑141 (MQ=255)                                                                                                                             cTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGCTAAAACACCTACCATTGCCGTAGCTACTAAAACCttt        >  1:618032/1‑141 (MQ=255)                                                                                                                              ttgaTTCATGTTATAAGAAGTTTGCATCATCATGTGGTACGAGTCCATTTAGAGAATAAACTTTCGCATCTGCGTGATTACCATATCCAACAACGCCAGTAGCTAAAGCGCCTACCATTGCCGGAGCGACGAAAACGtttt       <  2:309098/138‑1 (MQ=255)                                                                                                                               tAATTCATGTTATAAGAAGTTTGGATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGCTAAAACACCTACCATTGCCGTAGCTACTAAAACCttttt      <  2:255198/141‑1 (MQ=255)                                                                                                                                aaTTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGCTAAAACACCTACCATTGCCGTAGCTACTAAAACCTTTTTa     >  2:330963/1‑141 (MQ=255)                                                                                                                                  ttCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGCTAAAACACCTACCATTGCCGTAGCTACTAAAACCTTTTTaaa   >  1:187846/1‑141 (MQ=255)                                                                                                                                   tCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGCTAAAACACCTACCATTGCCGTAGCTACTAAAACCTTTTTAAAt  <  2:65455/141‑1 (MQ=255)                                                                                                                                          |                                                                                                                                      GCCGTCTTGTTTAGCTTTGTTGATAATGTCTTTGATTGAATGTAGTCCTTTATCGGCGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTCTAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGCTAAAACACCTACCATTGCCGTAGCTACTAAAACCTTTTTAAAT  >  CP000730/462877‑463147

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