BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,184,654	A→T	16.8%	N222Y (AAT→TAT)	yqiI →	fimbrial protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,184,654	0	A	T	16.8%	770.9 / 135.6	285	N222Y (AAT→TAT)	yqiI	fimbrial protein
Reads supporting (aligned to +/- strand): ref base A (124/113); new base T (33/15); total (157/128)
Fisher's exact test for biased strand distribution p-value = 3.96e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 291 total aligned reads displayed.

ATAAAGTAAATAAAAAACTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATATTTTAGGG  >  CP009273/3184505‑3184803

aTACAGTAAATAAAAAACTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTa                                                                                                                                                       <  1:2408847/150‑1 (MQ=255)

 tAAAGTAAATAAAAAACTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTaa                                                                                                                                                      >  1:4194883/1‑150 (MQ=255)

  aaaGTAAATAAAAAACTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAAt                                                                                                                                                     >  1:4389713/1‑150 (MQ=255)

  aaaGTAAATAAAAAACTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAAt                                                                                                                                                     >  1:1661179/1‑150 (MQ=255)

  aaaGTAAATAAAAAACTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGAGTAAATAACATTCATGTCCTTGACTGCTTTGTTAAt                                                                                                                                                     >  2:3082146/1‑150 (MQ=255)

   aaGTAAATAAAAAACTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATg                                                                                                                                                    >  1:390769/1‑150 (MQ=255)

    aGTAAATAAAAAACTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTTATGt                                                                                                                                                   >  1:441941/1‑150 (MQ=255)

    aGTAAATAAAAAACTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTTATGt                                                                                                                                                   >  1:2901397/1‑150 (MQ=255)

    aGTAAATAAAAAACTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGt                                                                                                                                                   >  1:2117989/1‑150 (MQ=255)

     gTAAATAAAAAACTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTTATGtt                                                                                                                                                  >  2:2421243/1‑150 (MQ=255)

          taaaAAACTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCt                                                                                                                                             >  2:1558161/1‑150 (MQ=255)

          taaaAAACTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCt                                                                                                                                             <  2:2505724/150‑1 (MQ=255)

           aaaaaaCTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTg                                                                                                                                            >  2:1030675/1‑150 (MQ=255)

             aaaaCTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGa                                                                                                                                          <  1:3958707/150‑1 (MQ=255)

             aaaaCTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGa                                                                                                                                          >  1:4045571/1‑150 (MQ=255)

               aaCTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAAc                                                                                                                                        <  2:1102163/150‑1 (MQ=255)

               aaCTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAAc                                                                                                                                        >  2:952281/1‑150 (MQ=255)

                 cTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCa                                                                                                                                      <  1:4626090/150‑1 (MQ=255)

                  ttGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAg                                                                                                                                     <  1:301960/150‑1 (MQ=255)

                   tGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGc                                                                                                                                    >  1:4163382/1‑150 (MQ=255)

                     tATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTg                                                                                                                                  <  1:1482012/150‑1 (MQ=255)

                        tACCTGACGAATATATAGCACTTGAAATTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATg                                                                                                                               <  1:1268604/150‑1 (MQ=255)

                          ccTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGc                                                                                                                             >  1:1884080/1‑150 (MQ=255)

                           cTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCg                                                                                                                            <  1:3365344/150‑1 (MQ=255)

                             gACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGtt                                                                                                                          >  1:4569136/1‑150 (MQ=255)

                              aCGAATATATAGCACTTGAATATAAAGGTATGAGCGGAGAAGGTGCTATAGAAGTAGATAAAAATCTGACATTAAGGATTCGCGGATTAAATAACATTCTTTTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTg                                                                                                                         >  2:350339/1‑150 (MQ=255)

                                 aaTATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCg                                                                                                                      >  2:294071/1‑150 (MQ=255)

                                    atataGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACt                                                                                                                   <  1:180853/150‑1 (MQ=255)

                                    atataGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACt                                                                                                                   >  1:3893298/1‑150 (MQ=255)

                                     tataGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACtt                                                                                                                  <  1:2973879/150‑1 (MQ=255)

                                      ataGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTTATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACttt                                                                                                                 >  2:733865/1‑150 (MQ=255)

                                      ataGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTTATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACttt                                                                                                                 >  2:3153748/1‑150 (MQ=255)

                                       taGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTg                                                                                                                <  1:1084823/150‑1 (MQ=255)

                                        aGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGTGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTgg                                                                                                               <  2:3636497/150‑1 (MQ=255)

                                        aGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTTATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTgg                                                                                                               >  1:2171663/1‑150 (MQ=255)

                                         gCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTTATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGt                                                                                                              >  2:4044510/1‑150 (MQ=255)

                                          cACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTa                                                                                                             >  2:279673/1‑150 (MQ=255)

                                               gAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAaaata                                                                                                        >  1:129291/1‑150 (MQ=255)

                                                aaTTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTTATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAaaataa                                                                                                       <  2:4086829/150‑1 (MQ=255)

                                                aaTTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAaaataa                                                                                                       <  1:930364/150‑1 (MQ=255)

                                                 aTTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAaaataaa                                                                                                      >  1:2668090/1‑150 (MQ=255)

                                                  tttAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAaaataaat                                                                                                     <  2:2994307/150‑1 (MQ=255)

                                                   ttAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAAtt                                                                                                    >  1:712072/1‑150 (MQ=255)

                                                   ttAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAAtt                                                                                                    >  1:4241285/1‑150 (MQ=255)

                                                     aaaGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTcc                                                                                                  >  2:2710456/1‑150 (MQ=255)

                                                      aaGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTccc                                                                                                 <  2:4542005/150‑1 (MQ=255)

                                                      aaGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTccc                                                                                                 <  1:3900569/150‑1 (MQ=255)

                                                      aaGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAAGTccc                                                                                                 <  2:3972321/150‑1 (MQ=255)

                                                          tATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAAc                                                                                             >  1:1523065/1‑150 (MQ=255)

                                                           aTGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACa                                                                                            >  2:1648392/1‑150 (MQ=255)

                                                               gCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTa                                                                                        <  2:2547675/150‑1 (MQ=255)

                                                                 ggAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTaaa                                                                                      <  2:1783532/150‑1 (MQ=255)

                                                                 ggAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTaaa                                                                                      >  1:3055301/1‑150 (MQ=255)

                                                                 ggAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTGTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTaaa                                                                                      >  1:3642524/1‑150 (MQ=255)

                                                                  gAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTTATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTaaaa                                                                                     >  1:2354515/1‑150 (MQ=255)

                                                                  gAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTaaaa                                                                                     <  2:1257579/150‑1 (MQ=255)

                                                                    gCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTTATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAAt                                                                                   >  1:323617/1‑150 (MQ=255)

                                                                     cAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATa                                                                                  >  1:3821139/1‑150 (MQ=255)

                                                                     cAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATa                                                                                  <  1:1952121/150‑1 (MQ=255)

                                                                     cAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATa                                                                                  >  1:3318483/1‑150 (MQ=255)

                                                                         tGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTTATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAg                                                                              >  1:3707267/1‑150 (MQ=255)

                                                                          gCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGc                                                                             >  2:2029020/1‑150 (MQ=255)

                                                                          gCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGc                                                                             <  1:489521/150‑1 (MQ=255)

                                                                           cTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTTATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTACAGAACAATTAAAAATACCAGCg                                                                            >  2:4427026/1‑150 (MQ=255)

                                                                               aGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCgtgag                                                                        <  1:3768685/150‑1 (MQ=255)

                                                                                gAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTTATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCgtgagt                                                                       >  2:385319/1‑150 (MQ=255)

                                                                                  aGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCgtgagtga                                                                     <  2:1354402/150‑1 (MQ=255)

                                                                                   gTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGtgagtgag                                                                    >  1:2723006/1‑150 (MQ=255)

                                                                                   gTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGtgagtgag                                                                    <  1:3479059/150‑1 (MQ=255)

                                                                                   gTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGtgagtgag                                                                    <  1:464275/150‑1 (MQ=255)

                                                                                     aGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGAc                                                                  <  2:4364373/150‑1 (MQ=255)

                                                                                                    cTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCAt                                                   >  1:885382/1‑150 (MQ=255)

                                                                                                     ttCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATg                                                  >  2:838943/1‑150 (MQ=255)

                                                                                                        aGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGAcc                                               >  1:2229733/1‑150 (MQ=255)

                                                                                                           aTTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCaaa                                            >  2:2965926/1‑150 (MQ=255)

                                                                                                                cgGATTAAATAACATTCATGTCCTTGACTGCTTTGTTTATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATcc                                       <  2:4235508/150‑1 (MQ=255)

                                                                                                                 ggATTAAATAACATTCATGTCCTTGACTGCTTTGTTTATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCg                                      <  2:3517344/150‑1 (MQ=255)

                                                                                                                       aaTAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCg                                <  1:823366/150‑1 (MQ=255)

                                                                                                                        aTAACATTCATGTCCTTGACTGCTTTGTTTATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCgg                               >  1:2716134/1‑150 (MQ=255)

                                                                                                                        aTAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCgg                               >  2:566666/1‑150 (MQ=255)

                                                                                                                          aaCATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTCGTCATGACCAAAGATCCGGGTGCGGcc                             >  2:2863648/1‑150 (MQ=255)

                                                                                                                             aTTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGt                          >  1:396040/1‑150 (MQ=255)

                                                                                                                                  tGTCCTTGACTGCTTTGTTTATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGa                     <  1:2421243/150‑1 (MQ=255)

                                                                                                                                      cTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAg                 <  2:2145388/150‑1 (MQ=255)

                                                                                                                                       ttGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAACGATCCGGGTGCGGCCTGTACTGAGCAGt                <  2:2612650/150‑1 (MQ=255)

                                                                                                                                        tGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGtt               >  2:716150/1‑150 (MQ=255)

                                                                                                                                         gACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGttt              <  2:326846/150‑1 (MQ=255)

                                                                                                                                           cTGCTTTGTTTATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTaa            >  2:977552/1‑150 (MQ=255)

                                                                                                                                            tGCTTTGTTTATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAat           <  1:3153546/150‑1 (MQ=255)

                                                                                                                                            tGCTTTGTTTATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAat           <  1:861913/150‑1 (MQ=255)

                                                                                                                                            tGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAat           >  1:550905/1‑150 (MQ=255)

                                                                                                                                            tGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAat           >  2:2023207/1‑150 (MQ=255)

                                                                                                                                              cTTTGTTAATGTTGATCTGGAACCAGCTTATGGCGTTGTAGACTTTGGTAAAATAAATTCCCGATCAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGATTAatat         <  2:2268467/150‑1 (MQ=255)

                                                                                                                                              cTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAatat         >  1:670969/1‑150 (MQ=255)

                                                                                                                                               tttGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATAtt        <  2:1861735/150‑1 (MQ=255)

                                                                                                                                               tttGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATAtt        <  2:159394/150‑1 (MQ=255)

                                                                                                                                               tttGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATAtt        >  1:2546791/1‑150 (MQ=255)

                                                                                                                                               tttGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATAtt        >  1:114739/1‑150 (MQ=255)

                                                                                                                                                 tGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATAtttt      >  1:601988/1‑150 (MQ=255)

                                                                                                                                                   ttAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATATTTTAg    <  1:979312/150‑1 (MQ=255)

ATAAAGTAAATAAAAAACTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATATTTTAGGG  >  CP009273/3184505‑3184803

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