BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,557,815	T→A	57.1%	K709N (AAA→AAT)	glgP ←	glycogen phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,557,815	0	T	A	57.1%	179.0 / inf	298	K709N (AAA→AAT)	glgP	glycogen phosphorylase
Reads supporting (aligned to +/- strand): ref base T (46/81); new base A (74/96); total (121/177)
Fisher's exact test for biased strand distribution p-value = 2.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 302 total aligned reads displayed.

CCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCA  >  CP009273/3557667‑3557960

   ccAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGt                                                                                                                                               <  1:4578342/150‑1 (MQ=255)

    cAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTa                                                                                                                                              <  1:2912717/150‑1 (MQ=255)

     aGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAg                                                                                                                                             >  2:910062/1‑150 (MQ=255)

     aGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAg                                                                                                                                             <  1:570088/150‑1 (MQ=255)

     aGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAg                                                                                                                                             <  2:1965373/150‑1 (MQ=255)

     aGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAg                                                                                                                                             >  2:3055896/1‑150 (MQ=255)

        aCCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATGTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGccc                                                                                                                                          >  2:4505377/1‑150 (MQ=255)

           tGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTga                                                                                                                                       >  1:732801/1‑150 (MQ=255)

            ggTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTgac                                                                                                                                      >  1:3600172/1‑150 (MQ=255)

                gTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGacgacg                                                                                                                                  >  1:2047147/1‑150 (MQ=255)

                  gATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCa                                                                                                                                <  2:2541398/150‑1 (MQ=255)

                  gATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCa                                                                                                                                <  2:545379/150‑1 (MQ=255)

                  gATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCa                                                                                                                                <  2:3116406/150‑1 (MQ=255)

                    tCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGt                                                                                                                              >  1:4355466/1‑150 (MQ=255)

                    tCGCCGAAGTTGATCAGAGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGt                                                                                                                              <  1:3544915/150‑1 (MQ=255)

                     cGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGtt                                                                                                                             >  1:3385034/1‑150 (MQ=255)

                     cGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGtt                                                                                                                             >  1:2211921/1‑150 (MQ=255)

                         gAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGttcttc                                                                                                                         >  1:341661/1‑150 (MQ=255)

                           aGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCa                                                                                                                       <  1:3975078/150‑1 (MQ=255)

                           aGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCa                                                                                                                       <  1:1648676/150‑1 (MQ=255)

                                aTCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCActtct                                                                                                                  >  1:235565/1‑150 (MQ=255)

                                aTCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCActtct                                                                                                                  <  2:4046556/150‑1 (MQ=255)

                                   aGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTcc                                                                                                               <  1:2300905/150‑1 (MQ=255)

                                    gCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCg                                                                                                              >  2:3175254/1‑150 (MQ=255)

                                    gCGAATCAACCAGATCGCGATAGCGACCCGGATCTTACGGACTGAATACACCGCTGCCGATTTTCGACAGCACCTGAAGCAGCTCCTCAAAATTCTCGAAGTATTCACGCGGGATGTATCCCAGACGACGCAGTTCATCCAATACATCCg                                                                                                              >  2:1284169/1‑150 (MQ=255)

                                      gAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCt                                                                                                            <  2:748747/150‑1 (MQ=255)

                                        aTCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCtgt                                                                                                          <  2:3384430/150‑1 (MQ=255)

                                        aTCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCtgt                                                                                                          >  2:4318546/1‑150 (MQ=255)

                                          cAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTgtgt                                                                                                        <  2:1430640/150‑1 (MQ=255)

                                           aaCCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGtt                                                                                                       <  1:3167847/150‑1 (MQ=255)

                                           aaCCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGtt                                                                                                       <  1:170934/150‑1 (MQ=255)

                                           aaCCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGtt                                                                                                       <  1:1216895/150‑1 (MQ=255)

                                            aCCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTa                                                                                                      <  2:3502811/150‑1 (MQ=255)

                                              cAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTAcc                                                                                                    <  2:1318554/150‑1 (MQ=255)

                                              cAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTAcc                                                                                                    <  2:4578121/150‑1 (MQ=255)

                                                gATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCaa                                                                                                  >  1:1857346/1‑150 (MQ=255)

                                                 aTCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCaaa                                                                                                 <  2:605825/150‑1 (MQ=255)

                                                 aTCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCaaa                                                                                                 <  2:2600980/150‑1 (MQ=255)

                                                    gcgATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAAt                                                                                              <  1:1253897/150‑1 (MQ=255)

                                                    gcgATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAAt                                                                                              <  2:3562326/150‑1 (MQ=255)

                                                         aGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGa                                                                                         >  2:4032826/1‑150 (MQ=255)

                                                          gCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGat                                                                                        <  1:769616/150‑1 (MQ=255)

                                                          gCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGat                                                                                        <  1:1408251/150‑1 (MQ=255)

                                                           cGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGata                                                                                       >  2:485382/1‑150 (MQ=255)

                                                            gACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGatat                                                                                      <  1:4475797/150‑1 (MQ=255)

                                                            gACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGatat                                                                                      >  1:1938911/1‑150 (MQ=255)

                                                            gACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGatat                                                                                      >  1:662791/1‑150 (MQ=255)

                                                            gACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGatat                                                                                      >  1:4252525/1‑150 (MQ=255)

                                                              cccGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTg                                                                                    >  2:4036037/1‑150 (MQ=255)

                                                                 ggATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCa                                                                                 >  1:4094753/1‑150 (MQ=255)

                                                                   aTCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGc                                                                               >  2:1105293/1‑150 (MQ=255)

                                                                    tCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCa                                                                              >  2:3864756/1‑150 (MQ=255)

                                                                         cGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGAc                                                                         >  1:1103104/1‑150 (MQ=255)

                                                                           gACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACAt                                                                       >  2:1670108/1‑150 (MQ=255)

                                                                           gACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACAt                                                                       >  1:2256591/1‑150 (MQ=255)

                                                                           gACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACAt                                                                       <  1:2511712/150‑1 (MQ=255)

                                                                              tGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGAt                                                                    <  2:2311464/150‑1 (MQ=255)

                                                                               gAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATc                                                                   >  2:4164507/1‑150 (MQ=255)

                                                                                   acacCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGc                                                               >  2:2122737/1‑150 (MQ=255)

                                                                                     acCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCAt                                                             >  1:3020362/1‑150 (MQ=255)

                                                                                       cGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCAtct                                                           >  1:3851309/1‑150 (MQ=255)

                                                                                       cGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCAtct                                                           <  2:4440493/150‑1 (MQ=255)

                                                                                        gCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCCTCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCAtctc                                                          <  1:2694490/150‑1 (MQ=255)

                                                                                        gCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCAtctc                                                          <  2:235565/150‑1 (MQ=255)

                                                                                         cTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCg                                                         <  1:312312/150‑1 (MQ=255)

                                                                                         cTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCg                                                         <  1:2629663/150‑1 (MQ=255)

                                                                                         cTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCg                                                         >  1:1587498/1‑150 (MQ=255)

                                                                                          tGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGa                                                        <  2:328167/150‑1 (MQ=255)

                                                                                                    cGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCAcc                                              >  2:1641368/1‑150 (MQ=255)

                                                                                                     gTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCg                                             <  2:2480887/150‑1 (MQ=255)

                                                                                                         gCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCa                                         <  2:2210076/150‑1 (MQ=255)

                                                                                                           aCCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAAc                                       <  1:1658884/150‑1 (MQ=255)

                                                                                                            ccTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACg                                      <  1:930860/150‑1 (MQ=255)

                                                                                                                 tGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCg                                 >  1:4033730/1‑150 (MQ=255)

                                                                                                                  gCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGa                                <  2:2880903/150‑1 (MQ=255)

                                                                                                                   cAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGAt                               <  2:4605995/150‑1 (MQ=255)

                                                                                                                   cAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGAt                               <  1:4051088/150‑1 (MQ=255)

                                                                                                                   cAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGAt                               >  2:2336712/1‑150 (MQ=255)

                                                                                                                    aGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATa                              >  1:2464507/1‑150 (MQ=255)

                                                                                                                      ctcctcATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGt                            <  1:1096563/150‑1 (MQ=255)

                                                                                                                      ctcctcATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGt                            >  2:4088482/1‑150 (MQ=255)

                                                                                                                       tcctcATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTc                           >  2:2333529/1‑150 (MQ=255)

                                                                                                                        cctcATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCa                          <  2:362130/150‑1 (MQ=255)

                                                                                                                         ctcATCTTTCTTGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAg                         >  1:2014807/1‑150 (MQ=255)

                                                                                                                         ctcATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAg                         <  2:2741512/150‑1 (MQ=255)

                                                                                                                           cATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAgcg                       >  2:3388753/1‑150 (MQ=255)

                                                                                                                             tCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCa                     <  2:1726979/150‑1 (MQ=255)

                                                                                                                              cTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCAc                    >  2:709385/1‑150 (MQ=255)

                                                                                                                              cTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCAc                    >  2:1552668/1‑150 (MQ=255)

                                                                                                                              cTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCAc                    >  1:893278/1‑150 (MQ=255)

                                                                                                                                   tcGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTa               <  1:1502607/150‑1 (MQ=255)

                                                                                                                                    cGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTACCACTTCTTCCGCTGTGTTACCAAAAATAAAGAGATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTaa              >  2:632669/1‑150 (MQ=255)

                                                                                                                                      tagtaTTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAgc            >  2:2060014/1‑150 (MQ=255)

                                                                                                                                      tagtaTTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAgc            <  1:4505377/150‑1 (MQ=255)

                                                                                                                                        gtaTTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAgcgc          >  2:209125/1‑150 (MQ=255)

                                                                                                                                        gtaTTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAgcgc          >  2:1019751/1‑150 (MQ=255)

                                                                                                                                        gtaTTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAgcgc          >  2:2965180/1‑150 (MQ=255)

                                                                                                                                          aTTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCTTCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCaa        <  2:2710867/150‑1 (MQ=255)

                                                                                                                                              aCGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACtt    >  1:3483143/1‑150 (MQ=255)

                                                                                                                                                gcgGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCa  <  1:2122737/150‑1 (MQ=255)

CCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCA  >  CP009273/3557667‑3557960

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