BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	1,330,955	G→T	63.3%	D290Y (GAT→TAT)	acnA →	aconitate hydratase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	1,330,955	0	G	T	63.3%	172.4 / inf	278	D290Y (GAT→TAT)	acnA	aconitate hydratase 1
Reads supporting (aligned to +/- strand): ref base G (57/45); new base T (81/95); total (138/140)
Fisher's exact test for biased strand distribution p-value = 1.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01

Only 100 of 293 total aligned reads displayed.

AGCCGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTG  >  CP009273/1330806‑1331104

aGCCGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTg                                                                                                                                                       >  2:4451806/1‑150 (MQ=255)

    ggTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATgg                                                                                                                                                   >  2:3671541/1‑150 (MQ=255)

     gTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTtatggttatggt                                                                                                                                                  <  1:1416077/150‑1 (MQ=255)

      tttCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTc                                                                                                                                                 <  1:3338421/150‑1 (MQ=255)

          cATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGa                                                                                                                                             >  2:3299399/1‑150 (MQ=255)

            tGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGAtt                                                                                                                                           <  2:281853/150‑1 (MQ=255)

              cTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCa                                                                                                                                         >  2:2764983/1‑150 (MQ=255)

               ttATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCAc                                                                                                                                        >  2:2179412/1‑150 (MQ=255)

               ttATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCAc                                                                                                                                        <  2:3406625/150‑1 (MQ=255)

               ttATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCAc                                                                                                                                        >  2:924874/1‑150 (MQ=255)

                 aTCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTa                                                                                                                                      <  2:2874697/150‑1 (MQ=255)

                 aTCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTa                                                                                                                                      <  2:661457/150‑1 (MQ=255)

                  tCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTAc                                                                                                                                     <  2:1860635/150‑1 (MQ=255)

                  tCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTAc                                                                                                                                     <  2:1854277/150‑1 (MQ=255)

                        aTGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTgg                                                                                                                               <  2:245161/150‑1 (MQ=255)

                          gTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCg                                                                                                                             >  2:2592483/1‑150 (MQ=255)

                            aGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGa                                                                                                                           >  1:153046/1‑150 (MQ=255)

                            aGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGa                                                                                                                           <  2:944179/150‑1 (MQ=255)

                              tGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATc                                                                                                                         >  1:604762/1‑150 (MQ=255)

                                ggCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATcgc                                                                                                                       >  1:4256172/1‑150 (MQ=255)

                                 gCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATcgcg                                                                                                                      <  1:673311/150‑1 (MQ=255)

                                 gCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATcgcg                                                                                                                      >  1:2728233/1‑150 (MQ=255)

                                   ttCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGcc                                                                                                                    >  1:1493157/1‑150 (MQ=255)

                                       aaCTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGccacca                                                                                                                >  2:3884484/1‑150 (MQ=255)

                                         cTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAATCCTGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCAtt                                                                                                              <  1:1725024/150‑1 (MQ=255)

                                         cTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCAtt                                                                                                              <  2:1239800/150‑1 (MQ=255)

                                          ttACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAACTGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTg                                                                                                             <  2:588346/150‑1 (MQ=255)

                                            aCCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGcc                                                                                                           >  1:4390964/1‑150 (MQ=255)

                                             ccGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCa                                                                                                          <  2:3648357/150‑1 (MQ=255)

                                               ggAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAat                                                                                                        >  2:2224570/1‑150 (MQ=255)

                                                    aTTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTc                                                                                                   <  1:3444011/150‑1 (MQ=255)

                                                     ttACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCg                                                                                                  <  1:1182231/150‑1 (MQ=255)

                                                      tACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGc                                                                                                 >  2:2158634/1‑150 (MQ=255)

                                                        cGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCa                                                                                               <  1:4522107/150‑1 (MQ=255)

                                                         gTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAg                                                                                              <  1:1787838/150‑1 (MQ=255)

                                                          tGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGa                                                                                             <  2:814300/150‑1 (MQ=255)

                                                            aaGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAat                                                                                           <  1:3309708/150‑1 (MQ=255)

                                                             aGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAata                                                                                          <  1:104267/150‑1 (MQ=255)

                                                              ggTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAatat                                                                                         >  1:2461382/1‑150 (MQ=255)

                                                                 aTTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGt                                                                                      >  2:1744750/1‑150 (MQ=255)

                                                                     ccGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCa                                                                                  <  1:386860/150‑1 (MQ=255)

                                                                     ccGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCa                                                                                  >  2:1401365/1‑150 (MQ=255)

                                                                      cGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCAc                                                                                 >  1:3879149/1‑150 (MQ=255)

                                                                        ccACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCt                                                                               <  1:1887801/150‑1 (MQ=255)

                                                                         cacaGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGGTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCtg                                                                              >  1:3067717/1‑150 (MQ=255)

                                                                          acaGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCtgt                                                                             >  2:3998725/1‑150 (MQ=255)

                                                                            aGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTgg                                                                           <  2:200351/150‑1 (MQ=255)

                                                                            aGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTgg                                                                           >  1:424883/1‑150 (MQ=255)

                                                                              aCCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATCTGTCGCCAGAATATGGTGCCACCTGTGGct                                                                         >  2:2727403/1‑150 (MQ=255)

                                                                                cTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGcttc                                                                       >  1:4007103/1‑150 (MQ=255)

                                                                                 tGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGcttct                                                                      <  2:4482630/150‑1 (MQ=255)

                                                                                   gTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCttcttc                                                                    >  2:2240586/1‑150 (MQ=255)

                                                                                   gTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCttcttc                                                                    >  1:3700633/1‑150 (MQ=255)

                                                                                    ttCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTcc                                                                   <  2:885763/150‑1 (MQ=255)

                                                                                    ttCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTcc                                                                   <  1:153131/150‑1 (MQ=255)

                                                                                       tcACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCaa                                                                >  2:2414941/1‑150 (MQ=255)

                                                                                        cACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAAt                                                               >  2:1405523/1‑150 (MQ=255)

                                                                                          cTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCg                                                             >  1:3834091/1‑150 (MQ=255)

                                                                                          cTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCg                                                             <  1:2805884/150‑1 (MQ=255)

                                                                                          cTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCg                                                             <  1:2593203/150‑1 (MQ=255)

                                                                                          cTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCg                                                             <  2:4092577/150‑1 (MQ=255)

                                                                                            gTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGAt                                                           >  2:4541550/1‑150 (MQ=255)

                                                                                              tACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGc                                                         <  2:1644577/150‑1 (MQ=255)

                                                                                              gaCCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGc                                                         <  2:1643953/149‑1 (MQ=255)

                                                                                                 ccAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGt                                                      >  1:4139704/1‑150 (MQ=255)

                                                                                                  cAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTa                                                     >  2:1065669/1‑150 (MQ=255)

                                                                                                      tgctgcGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAAccc                                                 <  2:3034146/150‑1 (MQ=255)

                                                                                                           cgcAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGAt                                            <  1:2586679/150‑1 (MQ=255)

                                                                                                           cgcAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGAt                                            >  2:254563/1‑150 (MQ=255)

                                                                                                           cgcAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGAt                                            >  2:3608661/1‑150 (MQ=255)

                                                                                                             cAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTa                                          >  2:4395838/1‑150 (MQ=255)

                                                                                                              aaaCATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTAc                                         >  1:4531978/1‑150 (MQ=255)

                                                                                                              aaaCATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTAc                                         >  1:2294231/1‑150 (MQ=255)

                                                                                                              aaaCATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTAc                                         >  2:1381027/1‑150 (MQ=255)

                                                                                                                  aTGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGc                                     <  2:3948755/150‑1 (MQ=255)

                                                                                                                  aTGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGc                                     <  2:1056964/150‑1 (MQ=255)

                                                                                                                   tGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCg                                    >  1:2899759/1‑150 (MQ=255)

                                                                                                                   tGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCg                                    >  1:3558631/1‑150 (MQ=255)

                                                                                                                      cggggtggGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGttt                                 >  1:2300872/4‑150 (MQ=255)

                                                                                                                         ggtggGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAg                              <  2:3088192/150‑1 (MQ=255)

                                                                                                                          gtggGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGc                             <  1:4456329/150‑1 (MQ=255)

                                                                                                                            ggggAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCgg                           <  1:2991707/150‑1 (MQ=255)

                                                                                                                               gAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCATCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGgcg                        >  1:1858689/1‑150 (MQ=255)

                                                                                                                                aaaTTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGgcgc                       >  2:111672/1‑150 (MQ=255)

                                                                                                                                 aaTTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCa                      <  1:107158/150‑1 (MQ=255)

                                                                                                                                      gtcgAATTTTATGGTTATGGTCTGGATTCACAACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGaa                 >  2:2663816/1‑150 (MQ=255)

                                                                                                                                        cgAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGa               >  2:2092237/1‑150 (MQ=255)

                                                                                                                                         gAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGAt              <  1:2042218/150‑1 (MQ=255)

                                                                                                                                         gAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGAt              <  2:1856587/150‑1 (MQ=255)

                                                                                                                                          aaTTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATc             <  2:1423503/150‑1 (MQ=255)

                                                                                                                                           aTTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCa            <  2:3878869/150‑1 (MQ=255)

                                                                                                                                           aTTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCa            >  1:2556331/1‑150 (MQ=255)

                                                                                                                                             tttATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAgg          <  2:1255749/150‑1 (MQ=255)

                                                                                                                                               tatggttatggtCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTc        >  1:2949458/1‑150 (MQ=255)

                                                                                                                                               tatggttatggtCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTc        <  2:151657/150‑1 (MQ=255)

                                                                                                                                                aTGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCg       <  1:2224570/150‑1 (MQ=255)

                                                                                                                                                aTGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCg       >  1:12272/1‑150 (MQ=255)

                                                                                                                                                aTGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGAGGGCGCAGCGAAGATCAGGTCg       >  1:3296908/1‑150 (MQ=255)

                                                                                                                                                 tGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGa      <  1:924885/150‑1 (MQ=255)

                                                                                                                                                     gATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTg  <  1:3643950/150‑1 (MQ=255)

AGCCGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTG  >  CP009273/1330806‑1331104

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