BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,055,546	C→A	10.3%	I446I (ATC→ATA)	scpA →	methylmalonyl‑CoA mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,055,546	0	C	A	10.3%	1145.7 / 77.2	320	I446I (ATC→ATA)	scpA	methylmalonyl‑CoA mutase
Reads supporting (aligned to +/- strand): ref base C (149/138); new base A (12/21); total (161/159)
Fisher's exact test for biased strand distribution p-value = 1.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 325 total aligned reads displayed.

TACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGC  >  CP009273/3055397‑3055695

tACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATc                                                                                                                                                       <  2:712056/150‑1 (MQ=255)

tACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATc                                                                                                                                                       <  2:1461078/150‑1 (MQ=255)

 aCATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCg                                                                                                                                                      <  1:1009607/150‑1 (MQ=255)

 aCATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCg                                                                                                                                                      >  2:4653659/1‑150 (MQ=255)

     tGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCa                                                                                                                                                  <  1:707105/150‑1 (MQ=255)

     tGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCa                                                                                                                                                  <  1:2533188/150‑1 (MQ=255)

        gTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAggg                                                                                                                                               >  1:507914/1‑150 (MQ=255)

            cTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAg                                                                                                                                           >  1:459907/1‑150 (MQ=255)

              gACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCg                                                                                                                                         <  2:3811938/150‑1 (MQ=255)

                 cGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCgtgt                                                                                                                                      <  1:3545739/150‑1 (MQ=255)

                  gATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTc                                                                                                                                     <  1:3112771/150‑1 (MQ=255)

                  gATCAAATCGTAAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTc                                                                                                                                     <  1:2436230/150‑1 (MQ=255)

                   aTCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCa                                                                                                                                    <  1:2056026/150‑1 (MQ=255)

                     cAAATCGTCAAAAAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGAGTCATc                                                                                                                                  >  2:4084269/1‑150 (MQ=255)

                      aaaTCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCg                                                                                                                                 >  1:4460994/1‑150 (MQ=255)

                       aaTCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGt                                                                                                                                >  1:1817085/1‑150 (MQ=255)

                        aTCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGtt                                                                                                                               <  1:2762858/150‑1 (MQ=255)

                           gTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGgt                                                                                                                            <  2:2080821/150‑1 (MQ=255)

                            tCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGgtg                                                                                                                           >  1:15851/1‑150 (MQ=255)

                              aaaCAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTc                                                                                                                         >  2:541333/1‑150 (MQ=255)

                              aaaCAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTc                                                                                                                         >  1:3504816/1‑150 (MQ=255)

                               aaCAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTca                                                                                                                        <  1:1883581/150‑1 (MQ=255)

                                aCAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTcaa                                                                                                                       >  2:1822955/1‑150 (MQ=255)

                                 cAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTcaac                                                                                                                      <  1:627268/150‑1 (MQ=255)

                                 cAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTcaac                                                                                                                      >  2:3253888/1‑150 (MQ=255)

                                 cAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTcaac                                                                                                                      >  2:3890695/1‑150 (MQ=255)

                                  aaGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTcaaca                                                                                                                     >  1:4697630/1‑150 (MQ=255)

                                     ccAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGt                                                                                                                  >  1:4072389/1‑150 (MQ=255)

                                       agagCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACACGTAc                                                                                                                >  2:4012891/1‑150 (MQ=255)

                                       agagCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTAc                                                                                                                >  1:1559288/1‑150 (MQ=255)

                                       agagCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTAc                                                                                                                >  2:2836905/1‑150 (MQ=255)

                                       agagCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTAc                                                                                                                >  1:4430437/1‑150 (MQ=255)

                                       agagCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTAc                                                                                                                <  1:3952156/150‑1 (MQ=255)

                                         agCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACaa                                                                                                              >  2:4535218/1‑150 (MQ=255)

                                           cTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAAc                                                                                                            <  1:1490601/150‑1 (MQ=255)

                                             attatCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTg                                                                                                          >  1:3157337/1‑150 (MQ=255)

                                             attatCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTg                                                                                                          <  1:958769/150‑1 (MQ=255)

                                              ttatCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTgg                                                                                                         >  2:3462120/1‑150 (MQ=255)

                                               tatCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGa                                                                                                        >  1:1923625/1‑150 (MQ=255)

                                                    aaCAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACg                                                                                                   >  2:2576681/1‑150 (MQ=255)

                                                       aGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAg                                                                                                <  1:3939870/150‑1 (MQ=255)

                                                         aTCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGAc                                                                                              >  1:92865/1‑150 (MQ=255)

                                                         aTCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGAc                                                                                              >  1:1726770/1‑150 (MQ=255)

                                                               gaAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAAcc                                                                                        >  1:4651499/1‑150 (MQ=255)

                                                                aaGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCg                                                                                       >  1:1995090/1‑150 (MQ=255)

                                                                  gCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGAt                                                                                     <  2:1274214/150‑1 (MQ=255)

                                                                   ccGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATg                                                                                    <  2:4482878/150‑1 (MQ=255)

                                                                       tGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACt                                                                                >  2:3686711/1‑150 (MQ=255)

                                                                       tGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACt                                                                                >  2:2300365/1‑150 (MQ=255)

                                                                       tGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACt                                                                                <  2:4351713/150‑1 (MQ=255)

                                                                       tGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATAGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACt                                                                                >  2:2927409/1‑150 (MQ=255)

                                                                          cATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATAGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTga                                                                             <  2:60637/150‑1 (MQ=255)

                                                                           aTGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGCTGTACTTgag                                                                            >  1:3688930/1‑150 (MQ=255)

                                                                              gCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATc                                                                         >  2:3515927/1‑150 (MQ=255)

                                                                               cGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATAGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCg                                                                        <  1:373304/150‑1 (MQ=255)

                                                                               cGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATAGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCg                                                                        <  1:374113/150‑1 (MQ=255)

                                                                                gAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATAGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGa                                                                       <  2:647873/150‑1 (MQ=255)

                                                                                      gATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGt                                                                 >  1:4447102/1‑150 (MQ=255)

                                                                                         cGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATAGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGAt                                                              <  1:4153745/150‑1 (MQ=255)

                                                                                            aGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGt                                                           <  2:1042063/150‑1 (MQ=255)

                                                                                            aGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGt                                                           <  2:1251157/150‑1 (MQ=255)

                                                                                             gCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTg                                                          <  1:214503/150‑1 (MQ=255)

                                                                                              cAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGc                                                         >  2:2172845/1‑150 (MQ=255)

                                                                                               aGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCg                                                        <  2:2325590/150‑1 (MQ=255)

                                                                                                ggTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGt                                                       >  2:92648/1‑150 (MQ=255)

                                                                                                 gTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTa                                                      <  1:1122263/150‑1 (MQ=255)

                                                                                                    tGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACg                                                   <  1:4300980/150‑1 (MQ=255)

                                                                                                       cAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGc                                                >  1:232161/1‑150 (MQ=255)

                                                                                                        aaaaCGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCa                                               >  1:3602528/1‑150 (MQ=255)

                                                                                                         aaaCGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCaa                                              >  2:1756645/1‑150 (MQ=255)

                                                                                                            cGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACCACGTGATGGTGCGTAACGAGCAAAtt                                           >  2:2491500/1‑150 (MQ=255)

                                                                                                            cGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAAtt                                           <  2:3152890/150‑1 (MQ=255)

                                                                                                              aaTGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGc                                         <  2:4121100/150‑1 (MQ=255)

                                                                                                                tGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCtt                                       <  2:1164763/150‑1 (MQ=255)

                                                                                                                 gATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATAGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTc                                      <  2:2650690/150‑1 (MQ=255)

                                                                                                                   tCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGc                                    >  2:2447552/1‑150 (MQ=255)

                                                                                                                     gAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTg                                  >  1:2950624/1‑150 (MQ=255)

                                                                                                                     gAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTg                                  >  2:320063/1‑150 (MQ=255)

                                                                                                                      aaGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTgg                                 >  1:1244884/1‑150 (MQ=255)

                                                                                                                      aaGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTgg                                 <  1:786508/150‑1 (MQ=255)

                                                                                                                       agagGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGa                                <  2:2023939/150‑1 (MQ=255)

                                                                                                                       agagGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGa                                >  1:3148335/1‑150 (MQ=255)

                                                                                                                            ccTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCa                           <  2:4627322/150‑1 (MQ=255)

                                                                                                                             cTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCAt                          <  1:817449/150‑1 (MQ=255)

                                                                                                                              tCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCAtt                         <  1:3840470/150‑1 (MQ=255)

                                                                                                                              tCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCAtt                         <  1:1449844/150‑1 (MQ=255)

                                                                                                                              tCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCAtt                         >  1:2293683/1‑150 (MQ=255)

                                                                                                                              tCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCAtt                         >  1:991085/1‑150 (MQ=255)

                                                                                                                              tCAGCGCGCGAACAGTCGCTGATAGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCAtt                         <  1:2881814/150‑1 (MQ=255)

                                                                                                                               cAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTc                        <  1:1225229/150‑1 (MQ=255)

                                                                                                                               cAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTc                        >  2:814451/1‑150 (MQ=255)

                                                                                                                               cAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTc                        <  1:1168540/150‑1 (MQ=255)

                                                                                                                               cAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTc                        <  1:3253888/150‑1 (MQ=255)

                                                                                                                               cAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTc                        <  2:2720135/150‑1 (MQ=255)

                                                                                                                               cAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACGGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTc                        <  2:3720074/150‑1 (MQ=255)

                                                                                                                                aGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTcg                       >  2:2909290/1‑150 (MQ=255)

                                                                                                                                aGCGCGCGAACAGTCGCTGATAGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTcg                       >  1:2657424/1‑150 (MQ=255)

                                                                                                                                              cGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGtgatg         <  1:3306806/150‑1 (MQ=255)

                                                                                                                                               gCTGATAGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGtgatga        <  2:2151417/150‑1 (MQ=255)

                                                                                                                                                 tGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGatgatg      >  2:331166/1‑150 (MQ=255)

TACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGC  >  CP009273/3055397‑3055695

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