BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	35,913	C→G	27.2%	S213S (TCG→TCC)	USA300HOU_0026 ←	IS431mec transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	35,913	0	C	G	27.2%	567.3 / inf	475	S213S (TCG→TCC)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base C (153/189); new base G (69/59); total (225/250)
Fisher's exact test for biased strand distribution p-value = 7.88e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

Only 100 of 484 total aligned reads displayed.

TTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTC  >  CP000730/35774‑36053

       agatccaaaataagatTGTGATCTAATAAACATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTGACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATcc                                                                                                                                      <  2:466239/135‑1 (MQ=17)

                    aagagtattttaataaAAAT‑GATGGTTCTGGTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAgaga                                                                                                                       <  1:309406/131‑1 (MQ=11)

                     agagtattttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                                                                                                      <  2:246904/132‑1 (MQ=11)

                     agagtattttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                                                                                                      <  2:723495/132‑1 (MQ=11)

                     agagtattttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                                                                                                      <  1:467858/132‑1 (MQ=11)

                      gagtattttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAcc                                                                                                                     <  2:927575/133‑1 (MQ=11)

                      gagtattttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAcc                                                                                                                     <  1:926845/133‑1 (MQ=11)

                         tgATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCt                                                                                                                    <  2:553618/141‑1 (MQ=255)

                         tgATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCt                                                                                                                    <  2:1099238/141‑1 (MQ=255)

                           ttttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑GAAAAATAGAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCGGAAGAGACCTGCgg                                                                                                                >  1:790543/4‑141 (MQ=11)

                           ttttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCCCTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg                                                                                                                <  2:67520/138‑1 (MQ=11)

                           ttttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg                                                                                                                <  1:452393/138‑1 (MQ=11)

                           ttttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg                                                                                                                >  1:152439/4‑141 (MQ=11)

                           ttttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg                                                                                                                <  1:871186/138‑1 (MQ=11)

                           ttttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTAAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGAGAAAATCAGTAGATCTGAAGAGAGCTGCgg                                                                                                                >  2:793457/4‑141 (MQ=11)

                           ttttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCGAATAACTAATATATCATGTAAGTGTTCGTTTAACTTGCTAGCATGATGCTAAATTCGTGGCACGGCGAAAAGCCGTAGATCTGAAGTGAcggtcgg                                                                                                                >  1:379696/4‑135 (MQ=11)

                            tttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt                                                                                                               <  1:741772/139‑1 (MQ=11)

                            tttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt                                                                                                               <  2:126752/139‑1 (MQ=11)

                             ttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGtt                                                                                                              <  1:437755/140‑1 (MQ=11)

                              taataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAACATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGttt                                                                                                             >  2:699047/1‑140 (MQ=11)

                                ataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAGAATCCGTAGATCTGAAGAGACCTGCGGTTCtt                                                                                                           >  1:341435/1‑141 (MQ=11)

                                ataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCtt                                                                                                           >  2:143822/1‑141 (MQ=11)

                                 taaAAATAGAAGGGGCTGTTGCGAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGGCCTGCGGTTCt                                                                                                            >  2:264607/1‑141 (MQ=255)

                                        aGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTATTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTtata                                                                                                     >  1:258374/1‑141 (MQ=17)

                                                                                 catctATCATGTCTGTGTTCGTTTAACCTGATACCATGATGCTAATTTCGTGGCATGGGGACAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTATTTg                                                            <  1:1059479/137‑1 (MQ=1)

                                                                                              aGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACTTGCTGTTCTTTTTGGATATAGCTTAAATACACTCAATACCTTTTATAGTATTCTTTGCTGTATTGATACt                                               <  1:808555/141‑1 (MQ=11)

                                                                                              aGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACt                                               >  1:588055/1‑141 (MQ=14)

                                                                                              aGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAACTATTCTTTGCTGTATTGATACt                                               <  2:807008/141‑1 (MQ=14)

                                                                                               gtgtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACtt                                              >  2:424321/1‑141 (MQ=255)

                                                                                               gtgtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACtt                                              <  1:394469/141‑1 (MQ=255)

                                                                                               gtgtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACtt                                              <  2:60095/141‑1 (MQ=255)

                                                                                               gtgtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACtt                                              <  1:219355/141‑1 (MQ=255)

                                                                                               gtgtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACtt                                              >  2:76628/1‑141 (MQ=255)

                                                                                               gtgtTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACtt                                              >  1:237782/1‑141 (MQ=32)

                                                                                                tgtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACttt                                             >  1:406257/1‑141 (MQ=255)

                                                                                                tgtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACttt                                             >  1:412649/1‑141 (MQ=17)

                                                                                                tgtTCGTTTAACTTGCTAGCATGATGCAAATGTCGTGGCATGGGGAAAATCCGTAGATCTGAGGAGACCTACGGTTCTTTGAATATAGAGCGTAAATACATTCTATACCTTTTAAAGTATTCTTTGCCGTATTGATACttt                                             >  2:1579/1‑141 (MQ=11)

                                                                                                tgtTCGTTTAACCTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACttt                                             <  1:125077/141‑1 (MQ=11)

                                                                                                  ttCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGa                                           <  2:78153/141‑1 (MQ=17)

                                                                                                   tCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACGCAAGACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGAt                                          >  1:84210/1‑141 (MQ=12)

                                                                                                    cgtttAAATTGCTAGCATGATGCTAATTTAGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTGTTTTTATATAGAGCGTAAATACACTCAATACATTTTCACGTTTCCGGGGCTTTATTGATGCTTTGATa                                         >  2:1055243/4‑141 (MQ=11)

                                                                                                    cGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATa                                         >  1:944883/1‑141 (MQ=17)

                                                                                                     gtttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATat                                        >  2:586127/3‑140 (MQ=255)

                                                                                                     gtttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGCCCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATAc                                        >  2:990396/3‑141 (MQ=12)

                                                                                                      tttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATAcc                                       <  1:542492/140‑1 (MQ=17)

                                                                                                      tttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATAcc                                       <  1:1025951/140‑1 (MQ=17)

                                                                                                      tttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTATTAAAGTATTCTTTGCTGTATTGATACTTTGATAcc                                       <  1:766612/140‑1 (MQ=12)

                                                                                                       ttAACTTGCCAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCt                                      <  2:190861/141‑1 (MQ=17)

                                                                                                        tAACTTGCTAGCATGATGCTGATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCTATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCtg                                     >  1:776448/1‑140 (MQ=11)

                                                                                                         aatctcgtgtcATGTTGCTAGTGTTGTGTCATTGCGAAAATCCGTAGAACTGAAGAGTCCTGCGTTTCTTTTTATATATAGCTTAAATACACTCAATCCCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTg                                    <  2:856790/131‑1 (MQ=11)

                                                                                                          aCTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGt                                   >  1:192090/1‑141 (MQ=17)

                                                                                                          aCTGGCAAGCATGATGTTATCTTCGTCGCATGGCGAAGATCCGGAGATCTGAAGAGCTCTGCGGGTCTTTTGATCTAGCGCGTAAATACACTCAATACCTTTTAAAGGCTTCTTTGCTGTATTGATACTTTGATACCTTGt                                   <  2:1038529/141‑1 (MQ=255)

                                                                                                           cTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtc                                  >  1:617482/1‑141 (MQ=17)

                                                                                                            ttGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtct                                 <  1:751064/141‑1 (MQ=32)

                                                                                                            ttGCTAGCATGATGCTAATGTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAGAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtct                                 >  1:719374/1‑141 (MQ=12)

                                                                                                              gCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATGTGTTGtcttt                               >  2:621508/1‑141 (MQ=11)

                                                                                                              gCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtcttt                               >  1:971919/1‑141 (MQ=17)

                                                                                                              gCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtcttt                               >  2:1078770/1‑141 (MQ=17)

                                                                                                               cTAGCATGATTCTAATTTCGAGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctttc                              <  2:474601/141‑1 (MQ=12)

                                                                                                               cTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctttc                              >  2:353521/1‑141 (MQ=17)

                                                                                                               cTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctttc                              >  2:964081/1‑141 (MQ=17)

                                                                                                                tAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctttct                             >  1:1005369/1‑141 (MQ=17)

                                                                                                                tAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctttct                             <  1:901455/141‑1 (MQ=17)

                                                                                                                  gCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTa                           >  2:650516/1‑141 (MQ=17)

                                                                                                                   cATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTAc                          >  1:727123/1‑141 (MQ=255)

                                                                                                                   cATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTAc                          <  2:224032/141‑1 (MQ=33)

                                                                                                                    atgatgCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATAGTTTGATCACTTGTATTTCTTACt                         >  2:209653/1‑141 (MQ=9)

                                                                                                                    atgatgCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATAGCTTGTCTTTCTTACt                         <  1:819991/141‑1 (MQ=21)

                                                                                                                    atgatgCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACt                         >  2:743369/1‑141 (MQ=18)

                                                                                                                     tgatgCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACtt                        >  2:1019293/1‑141 (MQ=255)

                                                                                                                     tgatgCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACtt                        <  1:283980/141‑1 (MQ=21)

                                                                                                                     agagtcTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACtt                        <  1:703362/136‑1 (MQ=11)

                                                                                                                         gCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAat                    >  1:609332/1‑141 (MQ=21)

                                                                                                                          cTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGATCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAata                   >  1:642263/1‑141 (MQ=11)

                                                                                                                            aaTTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATg                 >  2:310102/1‑141 (MQ=34)

                                                                                                                            aaTTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATg                 >  2:451909/1‑141 (MQ=21)

                                                                                                                            aaTTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATg                 <  1:425772/141‑1 (MQ=25)

                                                                                                                            aaTTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATg                 >  1:672067/1‑141 (MQ=21)

                                                                                                                             aTTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGa                <  1:363207/141‑1 (MQ=25)

                                                                                                                              tttCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGAc               >  1:1083581/1‑141 (MQ=25)

                                                                                                                                tCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTGTTTCTTACTTTAATATGACgg             >  2:511151/1‑141 (MQ=11)

                                                                                                                                 cGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGt            <  1:293894/141‑1 (MQ=255)

                                                                                                                                  ggggCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTg           <  2:493986/139‑1 (MQ=39)

                                                                                                                                  gTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATATAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGActgtt           <  2:393538/141‑5 (MQ=11)

                                                                                                                                  gTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTTTCTTTCTTACTTTAATATGACGGTg           >  1:172562/1‑141 (MQ=11)

                                                                                                                                  gTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTg           <  2:321992/141‑1 (MQ=255)

                                                                                                                                  gTGGCATGGCGAAACTCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATTACGGtt           >  2:376450/1‑140 (MQ=255)

                                                                                                                                  gTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGAggggg           >  1:500941/1‑139 (MQ=19)

                                                                                                                                  gTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTg           >  1:995463/1‑141 (MQ=38)

                                                                                                                                   tGGCATGGCGAAAATCCGTAGATCTGAAGGGACCTGGGGTTCTTTTTATATAGAGCGGATATACATTCAGTACCTTTTAAAGTATTCTGAGGTGTATTGAGACTGTGATGCGTGATCTGTCTTCATTTAATGTGACGGTGa          >  2:525258/1‑141 (MQ=2)

                                                                                                                                   ggggaTGGGGGAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATACAGAGCGTAAATACATTCAGTACATTTAAAAGTATTCTTTGCTGTATTGATACTTTGATGGCTAGTCTTTCTTACTTTACTATGACTGTGa          <  2:438877/137‑1 (MQ=11)

                                                                                                                                    ggCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTggt         >  1:1009167/1‑139 (MQ=37)

                                                                                                                                    ggCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGAt         >  1:959655/1‑141 (MQ=37)

                                                                                                                                     gCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATc        <  2:820109/141‑1 (MQ=255)

                                                                                                                                      cATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCt       <  1:79232/141‑1 (MQ=39)

                                                                                                                                        tGGCGAAAATCCTTAGATCTGAAGATACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTg     <  1:627011/141‑1 (MQ=255)

                                                                                                                                        agccGAAACTGCGTGGATATGAAGAGACAAGCGGTTCTTACTCCAGAGAGCATAAAGACATTCAATACCTATTAAAGTATTAATTGCTGTATTGATACTTTGATACCTTATCTTTCTTACTTTAATATGACGGTGAGCTTg     >  2:89310/4‑141 (MQ=255)

                                                                                                                                         ggCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGc    <  1:562917/141‑1 (MQ=255)

                                                                                                                                           cGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTc  >  2:182555/1‑141 (MQ=37)

                                                                                                                                           cGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTc  >  1:517354/1‑141 (MQ=38)

TTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTC  >  CP000730/35774‑36053

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