BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	36,471	A→G	22.0%	S27S (AGT→AGC)	USA300HOU_0026 ←	IS431mec transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	36,471	0	A	G	22.0%	224.3 / 61.7	139	S27S (AGT→AGC)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base A (51/57); new base G (13/18); total (64/75)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 142 total aligned reads displayed.

GCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTAT  >  CP000730/36332‑36611

gCCATTTGTAATAATCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCA‑TATTCTTGACCCCAACGGTAGACCG‑TTGAATAATGAACGTTTACACCACGTTCCCTTAGTATTCCAGATATATCACGATAAc                                                                                                                                             <  1:128837/141‑1 (MQ=11)

gCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAAc                                                                                                                                             >  1:124922/1‑141 (MQ=21)

gCCATTTGTAATAAGCTTTTTTATGCGTTTTCTTCCAAATTTGATATAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAAc                                                                                                                                             <  1:158398/141‑1 (MQ=18)

 ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGATCG‑TTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCt                                                                                                                                            >  1:105576/1‑141 (MQ=255)

 ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGATCG‑TTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCt                                                                                                                                            >  1:268128/1‑141 (MQ=255)

 ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGCTAAcc                                                                                                                                            >  2:647836/1‑140 (MQ=19)

 ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACt                                                                                                                                            >  2:595669/1‑141 (MQ=21)

 ccATTTGTAATAAGCTTTTTTATGATTTTTCTTCCAAAGTTGATATAAAATTGGGGCA‑TATTCTTGAACCCAACGGAAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCGGATATATCACGATAACt                                                                                                                                            >  1:426839/1‑141 (MQ=11)

  cATTTGTAATAAGCTTTTTTATGCTTTTTCTTTCAAATTTGATATAAAAGTTGGGCA‑TATTCTTGAACCCAACGGTAGATCG‑TTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTc                                                                                                                                           <  2:47083/141‑1 (MQ=11)

  cATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGATCG‑TTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTc                                                                                                                                           <  2:66244/141‑1 (MQ=11)

  cATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGATCG‑TTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTc                                                                                                                                           <  2:31454/141‑1 (MQ=255)

  cATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTc                                                                                                                                           <  1:228831/141‑1 (MQ=18)

  cATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTc                                                                                                                                           >  2:355511/1‑141 (MQ=25)

  cATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATAACGATAACTc                                                                                                                                           >  1:323577/1‑141 (MQ=21)

   aTTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAACTTTGATATAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCa                                                                                                                                          >  2:206630/1‑141 (MQ=25)

    tttGTAATCAGCTTTTTTATGCTTTTTCTTCCCAATTTGATCTACTATTATGGCA‑TACTCCTGAGCCCAACGGTTGACCG‑TTGAATGATGAACGGTTACACCACGCTCCCTTAAAATTTCAGATATATCACGATAACTcca                                                                                                                                         >  2:304558/1‑139 (MQ=11)

    tttGTAATAATCTTTTCTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCA‑TATCCTTGAACCCAACGGTAGATCG‑TTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCaa                                                                                                                                         <  1:49771/141‑1 (MQ=11)

     ttttaataaGCTTTTTTATGCTTTTTCTTCCACATTTGATATAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGATCG‑TTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAAt                                                                                                                                        <  2:349696/138‑1 (MQ=11)

     ttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTTATATAAAATTGGGGCA‑TATTATTGAACCCAACGGTAGATCG‑TTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAAt                                                                                                                                        <  2:200333/141‑1 (MQ=11)

     ttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAAt                                                                                                                                        <  1:445465/141‑1 (MQ=21)

     ttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAAt                                                                                                                                        >  1:134731/1‑141 (MQ=21)

     ttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAAt                                                                                                                                        <  1:396555/141‑1 (MQ=21)

     ttGTAAAAAGCTTTTTTATGCTTTTTCTACCAAATTTGATATAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAAt                                                                                                                                        <  2:503362/141‑1 (MQ=255)

      tttaataaGTTTTTTTATGCTTTTTCTTCCAAAGTGGATATAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATg                                                                                                                                       <  2:240638/139‑1 (MQ=255)

      tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGATCG‑TTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAAtc                                                                                                                                       >  2:28471/1‑140 (MQ=255)

      tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGATCG‑TTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATg                                                                                                                                       <  1:659163/141‑1 (MQ=255)

      tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGATCG‑TTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATg                                                                                                                                       <  1:323834/141‑1 (MQ=11)

      tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGATCG‑TTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATg                                                                                                                                       <  2:535239/141‑1 (MQ=255)

      tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATg                                                                                                                                       <  2:32369/141‑1 (MQ=21)

      tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATg                                                                                                                                       <  2:35781/141‑1 (MQ=21)

      tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCA‑TATTCTTGACCCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATg                                                                                                                                       <  2:148808/141‑1 (MQ=255)

       ggaataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCT‑TTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAACTCAATgc                                                                                                                                      <  1:227807/139‑2 (MQ=255)

       gTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCCCGTTCCCTTAATATTTCAGATATATCACGATAACTCAATgc                                                                                                                                      >  2:121550/1‑140 (MQ=16)

       gTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTACTATTTCAGATATATCACGATAACTCAATgc                                                                                                                                      >  2:86735/1‑140 (MQ=12)

       gTAATAAGCTTTTTTATGCCTTTTCTTCCAAATTTGATATAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGATCG‑TTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgc                                                                                                                                      >  1:306230/1‑140 (MQ=11)

        taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGATCG‑TTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgca                                                                                                                                     <  2:299070/141‑3 (MQ=11)

        taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGATCG‑TTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgca                                                                                                                                     <  2:282562/141‑3 (MQ=11)

        taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATgca                                                                                                                                     <  1:350882/141‑3 (MQ=17)

        taataaGCTTTTTTATGCTGTTTCTTCCAAATTTGATATAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGATCG‑TTGAATGATGAACGTTTACAGCACGGTGCCTTAATATTTCAGATATATCACGATAGCTCAATgca                                                                                                                                     >  2:138969/1‑139 (MQ=11)

         aataaGCTTTTTTATGCTTTTTCTTCCACATTTGATATAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGATCG‑TTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGAAATATCACGATAGCCCAAtccct                                                                                                                                    >  2:187439/1‑137 (MQ=11)

         aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGATCG‑TTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgcaa                                                                                                                                    >  2:591847/1‑138 (MQ=11)

          ataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCATTGtata                                                                                                                                   >  2:54104/1‑141 (MQ=255)

          ataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGtcta                                                                                                                                   >  2:378104/1‑141 (MQ=255)

                 tttttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCA‑TATTCTTGAACAAAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGa                                                                                                                            >  1:625013/1‑141 (MQ=25)

                    ttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtag                                                                                                                         <  1:100280/141‑1 (MQ=35)

                    ttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATAACTTAGAtag                                                                                                                         >  1:597979/1‑141 (MQ=21)

                            tttATTGCAAATTTGATACAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTACTATGTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAAc                                                                                                                 <  1:157664/141‑1 (MQ=255)

                              tcttcCAAATTTGATACAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACgg                                                                                                               >  2:547771/1‑141 (MQ=35)

                              tcttcCAAATTTGATACAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACgg                                                                                                               <  2:77742/141‑1 (MQ=35)

                                ttcCAAATTTGATACAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCt                                                                                                             <  1:349945/141‑1 (MQ=35)

                                  ccAAATTTGATACAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAAATAGTAGCCAACggatgc                                                                                                           >  2:478196/1‑137 (MQ=25)

                                  ccAAATTTGATACAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACCCAATGTATATCTTAGATAGTAGCCAACGGCTAc                                                                                                           >  1:197381/1‑141 (MQ=25)

                                  ccAAATTTGATACAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTAAAATATATCACGATAACTCAATGTATATCTTAGATGGTAGCCAACGGCTAc                                                                                                           >  1:286692/1‑141 (MQ=255)

                                   cAAATTTGATACAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACa                                                                                                          >  2:519093/1‑141 (MQ=35)

                                        ttGATACAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTATAGTGAt                                                                                                     <  2:368021/141‑1 (MQ=25)

                                        ttGATACAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGAt                                                                                                     <  1:265276/141‑1 (MQ=35)

                                        ttGATACAAAATTGGGGAATTGTTCTT‑AACCCAACGGTAGCTAGACCGAACG‑GGAACGTTCACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGAt                                                                                                     <  1:651257/141‑1 (MQ=11)

                                         tGATACAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCATAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATa                                                                                                    >  1:465783/1‑141 (MQ=25)

                                            tACAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACa                                                                                                 >  2:563279/1‑141 (MQ=35)

                                            tACAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATAATTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCCACGGCTACAGTGCTAACa                                                                                                 >  2:482998/1‑141 (MQ=255)

                                            tACAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACAG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAAcc                                                                                                 >  2:110381/1‑140 (MQ=21)

                                                 aaTTGGTGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGAAGACCGTTTTCACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCtt                                                                                            <  2:270466/141‑1 (MQ=255)

                                                    tGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTGACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGtt                                                                                         <  1:123703/141‑1 (MQ=17)

                                                        gCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAAt                                                                                     <  1:207029/141‑1 (MQ=32)

                                                          a‑tatTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTg                                                                                   >  1:263384/1‑141 (MQ=32)

                                                            tatTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGCTAACATCCTTGTTCAATTGt                                                                                  >  1:599675/1‑141 (MQ=17)

                                                               tCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTta                                                                               <  2:577180/141‑1 (MQ=32)

                                                                    aaCCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTACTATTTCAGATATATCACGATAACTCAATGTATCTCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTCAATTGTTTATATCt                                                                          >  2:24713/1‑141 (MQ=12)

                                                                    aaCCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCt                                                                          <  1:519093/141‑1 (MQ=32)

                                                                    aaCCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCt                                                                          >  1:551794/1‑141 (MQ=32)

                                                                      cccAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGa                                                                        <  1:37194/141‑1 (MQ=32)

                                                                            ggTAGAGCG‑TTGAATGATGAACGTTTACAGCCCGTGCCCTAAATATTTCAGATATATCACGATAACACAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTCGTTAAATTGTTTATATCTGAAATAGt                                                                  >  1:489431/1‑141 (MQ=11)

                                                                                gACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATAAAGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAAAATCCTTGTTAAATTGTTTATATCTGAACTAGttact                                                              >  1:136728/1‑138 (MQ=11)

                                                                                  ccG‑TTGAATGATGAACGTTTACTCCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATAc                                                            <  1:167627/141‑1 (MQ=17)

                                                                                    g‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACag                                                          >  2:591781/1‑141 (MQ=32)

                                                                                    g‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACag                                                          >  1:667046/1‑141 (MQ=32)

                                                                                    g‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTCTTTATATCTCAAATAGTTCATACag                                                          <  1:334069/141‑1 (MQ=17)

                                                                                         aaTGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACagaaga                                                      <  1:355511/141‑1 (MQ=32)

                                                                                               gaACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCtt                                                >  1:29860/1‑141 (MQ=32)

                                                                                                aaCGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCttt                                               <  2:258131/141‑1 (MQ=32)

                                                                                                aaCGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCttt                                               <  2:275691/141‑1 (MQ=32)

                                                                                                     ttACACCACGTTCCCTTAATATTTCCGCTATATTACGATCACTCAATGTAGATCTTCGATAGTAGCCAACGGCTACAGGGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACGGAAGACTCCTTTTTGtt                                          <  1:137284/141‑1 (MQ=11)

                                                                                                       acacCACGTTCCCTTAATATTTCAGATATATGCCGATAACTCAGTGTATGTCTTAGATAGTAGCCAACGGCTACAGTGATGACAGCCTTGTGAAATTGTTTATATCTGCAATAGTTCATACAGAATCAGCATTTTTGTTaa                                        >  2:282174/1‑141 (MQ=11)

                                                                                                         accacGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTaaaa                                      >  2:530222/1‑141 (MQ=32)

                                                                                                                    ttAATATTTCAGATATATCACGCTAACTCAATGTATAACCTAGATAGTAGCCAACGGCTACAGTGTTAACATCCTTGTCAAATTGTTTATATCTAAAATAGTTCATACAGAAGTGTCGTTTTTGTTAAAATAGTACTATaa                           >  2:48732/1‑141 (MQ=11)

                                                                                                                       atctttCAGCCATATCACGATAGCTCAATGCATATCTTAGATAGTAGCCAACTGCTACAGTGATAACATCCTTGTTAAATTGTTCATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAAtt                        <  1:554685/138‑1 (MQ=1)

                                                                                                                        tatTTCAGATATATCACGATAGCTCAACGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTc                       <  1:385854/141‑1 (MQ=255)

                                                                                                                         atTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCa                      <  1:201487/141‑1 (MQ=35)

                                                                                                                          tttCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCaa                     <  2:496702/141‑1 (MQ=35)

                                                                                                                          ttacAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCaa                     <  1:478306/138‑1 (MQ=25)

                                                                                                                            tCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACt                   >  1:174276/1‑141 (MQ=35)

                                                                                                                            tCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACt                   >  1:461650/1‑141 (MQ=35)

                                                                                                                                   tagcACGATAACTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAAc            <  2:183814/138‑1 (MQ=2)

                                                                                                                                     tCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAg          >  1:288255/1‑141 (MQ=35)

                                                                                                                                       aCGATAACTCAATGTATATCTTAGATAGTAGACAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTTTTTTTTTTCTATAAATTCGACTTTTGATCAGaa        >  2:140991/1‑141 (MQ=11)

                                                                                                                                         gATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAAcc      >  2:550425/1‑141 (MQ=34)

                                                                                                                                          aTAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCg     >  1:573116/1‑141 (MQ=35)

                                                                                                                                           tAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGt    <  1:638325/141‑1 (MQ=35)

                                                                                                                                           tAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTCTGCAACAGAACCGt    <  1:356719/141‑1 (MQ=37)

                                                                                                                                             aCTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTCAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAAACGtat  >  1:486273/1‑141 (MQ=255)

GCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTAT  >  CP000730/36332‑36611

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