BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,127,022	T→G	92.8%	E188D (GAA→GAC)	yijF ←	DUF1287 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,127,022	0	T	G	92.8%	1213.9 / 46.8	317	E188D (GAA→GAC)	yijF	DUF1287 family protein
Reads supporting (aligned to +/- strand): ref base T (17/6); new base G (154/140); total (171/146)
Fisher's exact test for biased strand distribution p-value = 5.20e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 323 total aligned reads displayed.

GATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTT  >  CP009273/4126875‑4127168

    gtgctgGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTcc                                                                                                                                              <  1:3244217/150‑1 (MQ=255)

    gtgctgGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTcc                                                                                                                                              <  1:2984321/150‑1 (MQ=255)

     tgctgGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCTTCTTCCt                                                                                                                                             <  2:1328132/150‑1 (MQ=255)

     tgctgGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCt                                                                                                                                             <  2:3518076/150‑1 (MQ=255)

             cTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCt                                                                                                                                     <  2:470192/150‑1 (MQ=255)

              tGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTg                                                                                                                                    >  1:2616009/1‑150 (MQ=255)

                         aGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTa                                                                                                                         <  2:3483763/150‑1 (MQ=255)

                            cGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGt                                                                                                                      >  2:4005875/1‑150 (MQ=255)

                            cGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGt                                                                                                                      >  2:4005458/1‑150 (MQ=255)

                            cGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGt                                                                                                                      >  1:2639723/1‑150 (MQ=255)

                             gTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTa                                                                                                                     >  2:3077110/1‑150 (MQ=255)

                               aaGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATc                                                                                                                   >  1:1063429/1‑150 (MQ=255)

                                aGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCa                                                                                                                  <  1:73895/150‑1 (MQ=255)

                                aGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCa                                                                                                                  >  2:2840886/1‑150 (MQ=255)

                                  aaaaaaGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCAcc                                                                                                                <  1:3104512/150‑1 (MQ=255)

                                   aaaaaGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCa                                                                                                               >  2:1588445/1‑150 (MQ=255)

                                   aaaaaGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCa                                                                                                               >  1:1294863/1‑150 (MQ=255)

                                   aaaaaGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCAACCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGACAATGTTATGTATCACCa                                                                                                               >  1:4404934/1‑150 (MQ=255)

                                    aaaaGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAg                                                                                                              <  1:1636125/150‑1 (MQ=255)

                                    aaaaGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAg                                                                                                              >  2:2742493/1‑150 (MQ=255)

                                    aaaaGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTGATGTATCACCAg                                                                                                              >  1:516935/1‑150 (MQ=255)

                                       agaagCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAgcgg                                                                                                           >  2:3168784/1‑150 (MQ=255)

                                        gaagCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCACTGTTATGTATCACCAgcggc                                                                                                          >  1:2483616/1‑150 (MQ=255)

                                         aagCCTTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAgcggcg                                                                                                         <  1:3131565/150‑1 (MQ=255)

                                           gCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTa                                                                                                       <  1:625835/150‑1 (MQ=255)

                                              gTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCg                                                                                                    >  1:702231/1‑150 (MQ=255)

                                               tCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGt                                                                                                   >  1:1538920/1‑150 (MQ=255)

                                                 aGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTcg                                                                                                 >  2:3183017/1‑150 (MQ=255)

                                                   tgtAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTcgcg                                                                                               >  2:1396051/1‑150 (MQ=255)

                                                   tgtAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTcgcg                                                                                               >  2:225013/1‑150 (MQ=255)

                                                      aaaaaGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGc                                                                                            >  2:2271125/1‑150 (MQ=255)

                                                      aaaaaGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGc                                                                                            >  1:3641883/1‑150 (MQ=255)

                                                         aaGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGaa                                                                                         >  2:4569689/1‑150 (MQ=255)

                                                         aaGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCGTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGaa                                                                                         >  2:4050498/1‑150 (MQ=255)

                                                            gATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGcc                                                                                      >  2:424665/1‑150 (MQ=255)

                                                             aTCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCa                                                                                     >  2:293961/1‑150 (MQ=255)

                                                              tCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCAt                                                                                    <  2:4788805/150‑1 (MQ=255)

                                                                tGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCt                                                                                  >  1:3709493/1‑150 (MQ=255)

                                                                  aaTTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGa                                                                                >  1:377047/1‑150 (MQ=255)

                                                                  aaTTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGa                                                                                >  2:260537/1‑150 (MQ=255)

                                                                    ttttttCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGaaa                                                                              >  1:3948386/1‑150 (MQ=255)

                                                                     tttttCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAac                                                                             >  2:1269281/1‑150 (MQ=255)

                                                                      ttttCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAacc                                                                            <  2:1145628/150‑1 (MQ=255)

                                                                       tttCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAacca                                                                           >  1:483594/1‑150 (MQ=255)

                                                                              aTGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAat                                                                    <  1:4445538/150‑1 (MQ=255)

                                                                               tgtgGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAata                                                                   >  2:3376396/1‑150 (MQ=255)

                                                                                gtgGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAatat                                                                  <  2:615416/150‑1 (MQ=255)

                                                                                   gggCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATAtgtg                                                               <  1:3077110/150‑1 (MQ=255)

                                                                                        aTCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGc                                                          >  2:1279228/1‑150 (MQ=255)

                                                                                         tctcTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGcc                                                         <  2:167238/150‑1 (MQ=255)

                                                                                           tctTATTTATCAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCg                                                       <  2:538187/150‑1 (MQ=255)

                                                                                             ttATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGtt                                                     <  2:3267155/150‑1 (MQ=255)

                                                                                             ttATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGtt                                                     <  1:1178874/150‑1 (MQ=255)

                                                                                             ttATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGtt                                                     <  1:4219564/150‑1 (MQ=255)

                                                                                              tATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTa                                                    >  2:863190/1‑150 (MQ=255)

                                                                                               aTTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTAt                                                   <  1:2034879/150‑1 (MQ=255)

                                                                                                tttAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATc                                                  <  1:3503027/150‑1 (MQ=255)

                                                                                                 ttAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATcc                                                 <  1:459345/150‑1 (MQ=255)

                                                                                                  tAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCa                                                <  2:3632437/150‑1 (MQ=255)

                                                                                                        aaTAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGTAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGcc                                          <  1:4655079/150‑1 (MQ=255)

                                                                                                        aaTAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGcc                                          <  1:2253918/150‑1 (MQ=255)

                                                                                                         aTAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCa                                         <  2:3872651/150‑1 (MQ=255)

                                                                                                          tAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTCTGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCCATATGTGCCAGCCCGTTCTCCAGTCGCCAg                                        >  1:2813635/1‑150 (MQ=255)

                                                                                                             cGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGgag                                     <  1:2844858/150‑1 (MQ=255)

                                                                                                              gATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTTTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGgaga                                    <  1:2156435/150‑1 (MQ=255)

                                                                                                              gATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGgaga                                    >  2:3841002/1‑150 (MQ=255)

                                                                                                              gATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGgaga                                    <  2:2543159/150‑1 (MQ=255)

                                                                                                              gATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGgaga                                    >  2:2664301/1‑150 (MQ=255)

                                                                                                               aTAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGAc                                   <  2:149648/150‑1 (MQ=255)

                                                                                                               aTAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGAc                                   >  1:851655/1‑150 (MQ=255)

                                                                                                                tAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACa                                  >  2:4520736/1‑150 (MQ=255)

                                                                                                                 aaTGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACaa                                 <  1:3839164/150‑1 (MQ=255)

                                                                                                                 aaTGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACaa                                 <  1:330547/150‑1 (MQ=255)

                                                                                                                  aTGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAat                                >  1:4013954/1‑150 (MQ=255)

                                                                                                                    gCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAatat                              <  1:777643/150‑1 (MQ=255)

                                                                                                                     ccccACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATc                             >  1:2724989/1‑150 (MQ=255)

                                                                                                                      cccACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCg                            >  2:286656/1‑150 (MQ=255)

                                                                                                                      cccACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCg                            >  2:1598969/1‑150 (MQ=255)

                                                                                                                        caccaTCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGcc                          <  1:756986/150‑1 (MQ=255)

                                                                                                                        caccaTCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGcc                          <  1:896529/150‑1 (MQ=255)

                                                                                                                         accaTCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGccc                         >  1:4063401/1‑150 (MQ=255)

                                                                                                                         accaTCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGccc                         >  1:3468241/1‑150 (MQ=255)

                                                                                                                          ccaTCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCg                        >  2:2460875/1‑150 (MQ=255)

                                                                                                                           caTCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGc                       >  2:2959312/1‑150 (MQ=255)

                                                                                                                           caTCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGc                       >  2:890367/1‑150 (MQ=255)

                                                                                                                           caTCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGc                       >  2:2155599/1‑150 (MQ=255)

                                                                                                                                gCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGAt                  <  2:3172435/150‑1 (MQ=255)

                                                                                                                                 ccAGTTAAACAGCACATCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATa                 >  1:157542/1‑150 (MQ=255)

                                                                                                                                  cAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATGTGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCGTGATAg                >  2:4350615/1‑150 (MQ=255)

                                                                                                                                   aGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGt               >  1:331703/1‑150 (MQ=255)

                                                                                                                                       aaaCAGCACATCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACt           >  2:2083544/1‑150 (MQ=255)

                                                                                                                                         aCAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTgg         <  1:98672/150‑1 (MQ=255)

                                                                                                                                          cAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTggg        >  1:2388501/1‑150 (MQ=255)

                                                                                                                                          cAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTggg        <  1:166023/150‑1 (MQ=255)

                                                                                                                                          cAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTggg        >  2:345819/1‑150 (MQ=255)

                                                                                                                                           aGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTgggg       >  2:1872972/1‑150 (MQ=255)

                                                                                                                                            gCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGt      <  1:1397594/150‑1 (MQ=255)

                                                                                                                                               caTCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCt   <  1:1035971/150‑1 (MQ=255)

                                                                                                                                               caTCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCt   <  2:2904890/150‑1 (MQ=255)

                                                                                                                                                aTCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCtt  <  2:1796805/150‑1 (MQ=255)

GATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTT  >  CP009273/4126875‑4127168

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