BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,615,947	G→T	91.4%	E73D (GAG→GAT)	radA →	DNA repair protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,615,947	0	G	T	91.4%	463.7 / 37.6	174	E73D (GAG→GAT)	radA	DNA repair protein
Reads supporting (aligned to +/- strand): ref base G (7/8); new base T (82/77); total (89/85)
Fisher's exact test for biased strand distribution p-value = 7.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

Only 100 of 174 total aligned reads displayed.

AGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCCA  >  CP009273/4615804‑4616090

agtgcCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGc                                                                                                                                             <  2:644886/148‑1 (MQ=255)

agtgcCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGc                                                                                                                                             <  2:535697/148‑1 (MQ=255)

   gcCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCcgc                                                                                                                                          >  2:1457378/1‑148 (MQ=255)

    ccTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCAATTGGTGGCGCGTAACTAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCcgcg                                                                                                                                         <  1:1232809/148‑1 (MQ=255)

     cTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGt                                                                                                                                        <  2:2322922/148‑1 (MQ=255)

     cTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGt                                                                                                                                        <  2:934462/148‑1 (MQ=255)

      tGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGtt                                                                                                                                       >  2:2374911/1‑148 (MQ=255)

      tGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGtt                                                                                                                                       <  1:2844691/148‑1 (MQ=255)

       gTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATCTAAGCCTAGAAGATCTGCCGCGttt                                                                                                                                      >  2:1216839/1‑148 (MQ=255)

        tCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGtttt                                                                                                                                     <  1:911225/148‑1 (MQ=255)

           tGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTcc                                                                                                                                  >  1:1056305/1‑148 (MQ=255)

           tGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTcc                                                                                                                                  <  2:3283885/148‑1 (MQ=255)

            gCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCa                                                                                                                                 <  2:3145868/148‑1 (MQ=255)

            gCCCGGCACCCCATCACCGAGGGGCGTCTTGCGGCGTCGCCAAGGGTGGCGCGTAACGAGCGTCTCAGCGGCAATGCCGGAAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCa                                                                                                                                 <  2:3145395/148‑1 (MQ=255)

               tGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCg                                                                                                                              >  1:990277/1‑148 (MQ=255)

                  aaCACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGAt                                                                                                                           <  2:1266969/148‑1 (MQ=255)

                   acacCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGAtt                                                                                                                          >  2:2809809/1‑148 (MQ=255)

                    cacCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGAttt                                                                                                                         <  1:1380779/148‑1 (MQ=255)

                       cATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTaaa                                                                                                                      >  1:1339230/1‑148 (MQ=255)

                       cATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTaaa                                                                                                                      >  1:423995/1‑148 (MQ=255)

                         tCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAaga                                                                                                                    >  1:3066821/1‑148 (MQ=255)

                         tCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAaga                                                                                                                    <  2:362463/148‑1 (MQ=255)

                          cACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAagag                                                                                                                   >  1:1575655/1‑148 (MQ=255)

                           aCCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGt                                                                                                                  <  2:1557557/148‑1 (MQ=255)

                              gAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCg                                                                                                               <  1:784069/148‑1 (MQ=255)

                               aGGTGCGTATTGCTGCGTCGACAATGGTGGCGCGTAACGAGCGTCTCAGAGTCTATGCCGGTAGCGCCGGGGTGGCAAAATTCCAGAAACTCTCCTATATCAGCCTTTAAGATCTGCCGCGTTTTTCCACCGTATTTAAAGAGTTCGa                                                                                                              <  2:583184/148‑1 (MQ=255)

                                ggTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGAc                                                                                                             >  1:371788/1‑148 (MQ=255)

                                     gTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGt                                                                                                        <  1:2179565/148‑1 (MQ=255)

                                     gTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGt                                                                                                        <  1:2359413/148‑1 (MQ=255)

                                      tCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTa                                                                                                       >  1:1503672/1‑148 (MQ=255)

                                      tCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTa                                                                                                       <  1:302503/148‑1 (MQ=255)

                                       cTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTAc                                                                                                      >  2:934630/1‑148 (MQ=255)

                                       cTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTAc                                                                                                      >  1:2535668/1‑148 (MQ=255)

                                       cTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTAc                                                                                                      >  1:1984811/1‑148 (MQ=255)

                                       cTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTAc                                                                                                      >  2:1024874/1‑148 (MQ=255)

                                                 cGCCAATGGTGGCGCGTAACGAGCGTCGCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAggcggcgg                                                                                            <  1:1856497/148‑1 (MQ=255)

                                                        ggtggCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCgtggtg                                                                                     <  1:2156782/148‑1 (MQ=255)

                                                         gtggCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGc                                                                                    >  2:1550745/1‑148 (MQ=255)

                                                                 tAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGt                                                                            >  1:2807816/1‑148 (MQ=255)

                                                                  aaCGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTg                                                                           >  2:1280749/1‑148 (MQ=255)

                                                                   aCGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTATGCGGCGGCGTGGTGCCAGGAAGTGc                                                                          >  1:785900/1‑148 (MQ=255)

                                                                   aCGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGc                                                                          <  2:2760692/148‑1 (MQ=255)

                                                                       gCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCAtt                                                                      <  2:1927829/148‑1 (MQ=255)

                                                                       gCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCAtt                                                                      >  1:1288824/1‑148 (MQ=255)

                                                                       gCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCAtt                                                                      >  1:2226831/1‑148 (MQ=255)

                                                                        cGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTc                                                                     <  2:984846/148‑1 (MQ=255)

                                                                           ctcAGCGGCTATGCGTGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGACATTCTGa                                                                  <  1:2146110/148‑1 (MQ=255)

                                                                           ctcAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGa                                                                  >  1:1278526/1‑148 (MQ=255)

                                                                           ctcAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGa                                                                  <  1:1303437/148‑1 (MQ=255)

                                                                           ctcAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGa                                                                  >  1:1347942/1‑148 (MQ=255)

                                                                           ctcAGCGGCTAGGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGa                                                                  <  2:2863556/148‑1 (MQ=255)

                                                                            tcAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGTCGTGGTGCCAGGAAGTGCCATTCTGAt                                                                 >  1:667364/1‑148 (MQ=255)

                                                                             cAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGAtt                                                                <  2:2364684/148‑1 (MQ=255)

                                                                               gCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTgg                                                              >  1:935723/1‑148 (MQ=255)

                                                                                cGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGc                                                             <  2:2658981/148‑1 (MQ=255)

                                                                                cGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGc                                                             >  2:2045130/1‑148 (MQ=255)

                                                                                cGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGc                                                             <  2:1626038/148‑1 (MQ=255)

                                                                                     aTGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTaa                                                        <  1:622436/148‑1 (MQ=255)

                                                                                     aTGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTaa                                                        <  1:622484/148‑1 (MQ=255)

                                                                                       gCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAAcc                                                      <  1:1809163/148‑1 (MQ=255)

                                                                                       gCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAAcc                                                      >  2:1800381/1‑148 (MQ=255)

                                                                                         cGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCt                                                    <  1:1532948/148‑1 (MQ=255)

                                                                                          ggTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCCGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTg                                                   <  1:793331/148‑1 (MQ=255)

                                                                                          ggTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTg                                                   <  1:793368/148‑1 (MQ=255)

                                                                                          ggTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTg                                                   <  1:2905515/148‑1 (MQ=255)

                                                                                          ggTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTg                                                   <  1:1861076/148‑1 (MQ=255)

                                                                                           gTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTgg                                                  <  1:2589531/148‑1 (MQ=255)

                                                                                           gTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTgg                                                  <  2:671504/148‑1 (MQ=255)

                                                                                            tAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGt                                                 <  2:1575655/148‑1 (MQ=255)

                                                                                             aGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTg                                                <  2:1751044/148‑1 (MQ=255)

                                                                                             aGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTg                                                <  1:112855/148‑1 (MQ=255)

                                                                                                  cGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCATACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCgggg                                           <  2:838755/148‑1 (MQ=255)

                                                                                                    gggTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAATGTCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGACAGTCTGATTGGCGGTAACCCTGGTGCGGGGaa                                         >  1:1696708/1‑148 (MQ=255)

                                                                                                      gTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAAt                                       <  1:122566/148‑1 (MQ=255)

                                                                                                        ggCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATcc                                     <  1:2454692/148‑1 (MQ=255)

                                                                                                         gCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCa                                    >  1:212538/1‑148 (MQ=255)

                                                                                                         gCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCa                                    >  2:1433188/1‑148 (MQ=255)

                                                                                                          cAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCAc                                   >  2:956980/1‑148 (MQ=255)

                                                                                                           aaaaGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACg                                  <  1:1940571/148‑1 (MQ=255)

                                                                                                            aaaGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACgc                                 >  1:183329/1‑148 (MQ=255)

                                                                                                             aaGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACgct                                >  1:1634067/1‑148 (MQ=255)

                                                                                                              aGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACgctg                               >  1:2159625/1‑148 (MQ=255)

                                                                                                              aGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACgctg                               >  2:3025535/1‑148 (MQ=255)

                                                                                                               gTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACgctgc                              >  2:598246/1‑148 (MQ=255)

                                                                                                                tCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACgctgct                             >  2:2494608/1‑148 (MQ=255)

                                                                                                                tCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACgctgat                             >  1:88123/1‑146 (MQ=255)

                                                                                                                 ccAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTa                            >  1:173832/1‑148 (MQ=255)

                                                                                                                 ccAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTa                            <  2:265259/148‑1 (MQ=255)

                                                                                                                 ccAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGCCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTa                            <  1:3088264/148‑1 (MQ=255)

                                                                                                                    gAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTg                         >  2:2064074/1‑148 (MQ=255)

                                                                                                                     aaaCTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGc                        >  2:1120661/1‑148 (MQ=255)

                                                                                                                       aCTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCaa                      >  2:2160763/1‑148 (MQ=255)

                                                                                                                         tctcCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAGATCCACGCTGCTACTGCAAAc                    >  2:3262472/1‑148 (MQ=255)

                                                                                                                          ctcCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACg                   >  2:2988837/1‑148 (MQ=255)

                                                                                                                              gATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCtgt               >  1:1029960/1‑148 (MQ=255)

                                                                                                                                  tCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCaa           <  2:1802349/148‑1 (MQ=255)

                                                                                                                                   cAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCaaa          >  2:2739169/1‑148 (MQ=255)

                                                                                                                                         tGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGcc    >  2:387915/1‑148 (MQ=255)

                                                                                                                                         tGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGcc    <  1:743731/148‑1 (MQ=255)

                                                                                                                                           aaGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCca  >  1:2400075/1‑148 (MQ=255)

AGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCCA  >  CP009273/4615804‑4616090

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