BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	2,237,098	G→T	94.2%	H12N (CAT→AAT)	folE ←	GTP cyclohydrolase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	2,237,098	0	G	T	94.2%	875.6 / 42.0	276	H12N (CAT→AAT)	folE	GTP cyclohydrolase I
Reads supporting (aligned to +/- strand): ref base G (6/10); new base T (143/117); total (149/127)
Fisher's exact test for biased strand distribution p-value = 2.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.51e-01

Only 100 of 276 total aligned reads displayed.

TCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACATAATCCAGTATTCACAG  >  CP009273/2236951‑2237246

tCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTaa                                                                                                                                                    >  2:4725569/1‑150 (MQ=255)

 cAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAAc                                                                                                                                                   <  1:1078889/150‑1 (MQ=255)

 cAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAAc                                                                                                                                                   >  1:4890851/1‑150 (MQ=255)

  aaaCTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAAcc                                                                                                                                                  >  2:731185/1‑150 (MQ=255)

  aaaCTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAAcc                                                                                                                                                  >  1:4069354/1‑150 (MQ=255)

  aaaCTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAAcc                                                                                                                                                  >  2:4428562/1‑150 (MQ=255)

    aCTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAg                                                                                                                                                <  1:4233650/150‑1 (MQ=255)

     cTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAgg                                                                                                                                               >  2:1346693/1‑150 (MQ=255)

      tGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAggg                                                                                                                                              <  1:3851455/150‑1 (MQ=255)

      tGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAggg                                                                                                                                              <  2:637282/150‑1 (MQ=255)

       gTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGc                                                                                                                                             <  2:431099/150‑1 (MQ=255)

        tcatcaGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGcc                                                                                                                                            <  2:1522469/150‑1 (MQ=255)

            caGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGctt                                                                                                                                        <  2:5055166/150‑1 (MQ=255)

              gCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGcttct                                                                                                                                      <  2:268764/150‑1 (MQ=255)

              gCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGcttct                                                                                                                                      >  1:3492849/1‑150 (MQ=255)

               ccAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCCGGGCCGcttctt                                                                                                                                     >  2:3460341/1‑150 (MQ=255)

               ccAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGcttctt                                                                                                                                     <  1:4442885/150‑1 (MQ=255)

               ccAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATATCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGcttctt                                                                                                                                     >  2:3317655/1‑150 (MQ=255)

                cAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCttt                                                                                                                                    >  1:3402938/1‑150 (MQ=255)

                 aGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTa                                                                                                                                   <  1:4702137/150‑1 (MQ=255)

                   gTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACt                                                                                                                                 <  1:3199202/150‑1 (MQ=255)

                    tCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTg                                                                                                                                >  2:4386721/1‑150 (MQ=255)

                     cGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGa                                                                                                                               <  1:5055008/150‑1 (MQ=255)

                     cGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGa                                                                                                                               >  1:4723684/1‑150 (MQ=255)

                     cGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGa                                                                                                                               >  1:3267902/1‑150 (MQ=255)

                      gagaTTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAg                                                                                                                              >  1:2301356/1‑150 (MQ=255)

                           tCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATg                                                                                                                         <  1:2463950/150‑1 (MQ=255)

                               cagcTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCACCTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCAt                                                                                                                     >  2:2571966/1‑150 (MQ=255)

                                agcTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCAtt                                                                                                                    >  1:4103235/1‑150 (MQ=255)

                                  cTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTa                                                                                                                  <  1:4671324/150‑1 (MQ=255)

                                   tGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTAt                                                                                                                 >  2:2258198/1‑150 (MQ=255)

                                    gCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATg                                                                                                                >  1:3483126/1‑150 (MQ=255)

                                       tGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGAtt                                                                                                             >  2:860934/1‑150 (MQ=255)

                                        gATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGAttt                                                                                                            >  1:2060003/1‑150 (MQ=255)

                                        gATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGAttt                                                                                                            <  1:2803070/150‑1 (MQ=255)

                                        gATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGAttt                                                                                                            <  1:4393949/150‑1 (MQ=255)

                                          tttCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTtct                                                                                                          <  1:3732866/150‑1 (MQ=255)

                                          tttCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTtct                                                                                                          >  2:2379848/1‑150 (MQ=255)

                                           ttCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTtctc                                                                                                         <  2:2151564/150‑1 (MQ=255)

                                           ttCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTtctc                                                                                                         <  1:129457/150‑1 (MQ=255)

                                              ggTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTg                                                                                                      >  1:153022/1‑150 (MQ=255)

                                                     tGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGgtgt                                                                                               >  2:672356/1‑150 (MQ=255)

                                                     tGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGgtgt                                                                                               <  2:402124/150‑1 (MQ=255)

                                                      gACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGtgtg                                                                                              <  1:3936924/150‑1 (MQ=255)

                                                       aCCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGa                                                                                             >  2:407640/1‑150 (MQ=255)

                                                        ccAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGAc                                                                                            <  1:1539270/150‑1 (MQ=255)

                                                             aaTAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTc                                                                                       >  2:3655131/1‑150 (MQ=255)

                                                              aTAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTcc                                                                                      >  1:572971/1‑150 (MQ=255)

                                                                   gCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCt                                                                                 <  1:2453512/150‑1 (MQ=255)

                                                                    cTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTg                                                                                >  2:4567025/1‑150 (MQ=255)

                                                                     ttttGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTGCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGc                                                                               <  1:3361182/150‑1 (MQ=255)

                                                                          cgcgTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTggg                                                                          >  2:4583589/1‑150 (MQ=255)

                                                                          cgcgTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTggg                                                                          >  1:1155857/1‑150 (MQ=255)

                                                                          cgcgTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTggg                                                                          >  1:3646633/1‑150 (MQ=255)

                                                                          cgcgTTTCGTTATCCATTTCATGAACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTTATGGCATGTATTATTTATACTGCAGGTGTGACGCATCCGCCATGCTTgtg                                                                          >  2:2064721/1‑148 (MQ=255)

                                                                           gcgTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTgggg                                                                         <  1:1762331/150‑1 (MQ=255)

                                                                            cgTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTggggg                                                                        <  1:2591858/150‑1 (MQ=255)

                                                                              tttCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCa                                                                      >  2:2045427/1‑150 (MQ=255)

                                                                                tCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCaaa                                                                    <  2:4828913/150‑1 (MQ=255)

                                                                                 cGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCGTTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCGAAg                                                                   >  1:3025044/1‑150 (MQ=255)

                                                                                  gTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGt                                                                  >  1:2980059/1‑150 (MQ=255)

                                                                                  gTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACAGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGt                                                                  >  1:4751358/1‑150 (MQ=255)

                                                                                      tCCATGTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGttatt                                                              <  2:4533518/150‑1 (MQ=255)

                                                                                        cATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTaa                                                            <  2:3460687/150‑1 (MQ=255)

                                                                                         aTTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAg                                                           >  2:3783190/1‑150 (MQ=255)

                                                                                          tttCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGc                                                          <  1:712001/150‑1 (MQ=255)

                                                                                           ttCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCt                                                         >  1:1183475/1‑150 (MQ=255)

                                                                                           ttCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCt                                                         <  1:1595179/150‑1 (MQ=255)

                                                                                             cATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGa                                                       >  2:3433820/1‑150 (MQ=255)

                                                                                                gCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGAttt                                                    >  2:2313438/1‑150 (MQ=255)

                                                                                                gCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGAttt                                                    <  1:82443/150‑1 (MQ=255)

                                                                                                  aCGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTAc                                                  >  2:1473251/1‑150 (MQ=255)

                                                                                                   cgggcgggCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACa                                                 <  1:1482859/150‑1 (MQ=255)

                                                                                                   cgggcgggCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACa                                                 >  2:3246352/1‑150 (MQ=255)

                                                                                                    gggcgggcGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAg                                                <  1:4005272/150‑1 (MQ=255)

                                                                                                       cgggcgCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCt                                             <  2:2843582/150‑1 (MQ=255)

                                                                                                        gggcgCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATTATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTg                                            >  1:481923/1‑150 (MQ=255)

                                                                                                              aGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTTCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGAGTTACAGCCTGattatt                                      >  1:4606524/1‑150 (MQ=255)

                                                                                                              aGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGattatt                                      >  2:1330858/1‑150 (MQ=255)

                                                                                                              aGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTGCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCACCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGattagt                                      >  1:2443920/1‑148 (MQ=255)

                                                                                                                   tgtTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCg                                 <  1:3618724/150‑1 (MQ=255)

                                                                                                                   tgtTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTATTTACTGAGTGATGGCATTTATGATTTATCCTGCAGGTGTGACGCCGCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCg                                 >  1:4501556/1‑150 (MQ=255)

                                                                                                                       tCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAg                             <  2:1076913/150‑1 (MQ=255)

                                                                                                                            tCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGGGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAggcggc                        >  1:5063123/1‑150 (MQ=255)

                                                                                                                             ccTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGgcggcg                       <  2:3066100/150‑1 (MQ=255)

                                                                                                                             ccTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGgcggcg                       <  1:3613462/150‑1 (MQ=255)

                                                                                                                               tCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGca                     <  2:178085/150‑1 (MQ=255)

                                                                                                                                 gcgcAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGcaca                   <  1:2908778/150‑1 (MQ=255)

                                                                                                                                  cgcAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACAt                  >  1:1992942/1‑150 (MQ=255)

                                                                                                                                  cgcAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACAt                  <  1:301137/150‑1 (MQ=255)

                                                                                                                                   gcAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACATa                 >  2:1168464/1‑150 (MQ=255)

                                                                                                                                          aCGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACATAATCCAGt          >  2:194697/1‑150 (MQ=255)

                                                                                                                                           cGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACATAATCCAGTa         <  2:3738656/150‑1 (MQ=255)

                                                                                                                                           cGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACATAATCCAGTa         >  1:2138393/1‑150 (MQ=255)

                                                                                                                                            gCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGTTGGCATTTATGATTTCAACTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACATAATACAGTAt        >  1:4310509/1‑150 (MQ=255)

                                                                                                                                            gCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCCGGGGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACATAATCCAGTAt        >  1:3276894/1‑150 (MQ=255)

                                                                                                                                            gCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACATAATCCAGTAt        <  1:1039223/150‑1 (MQ=255)

                                                                                                                                             cTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACATAATCCAGTAtt       <  1:134412/150‑1 (MQ=255)

                                                                                                                                                cATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACATAATCCAGTATTcac    <  2:2450858/150‑1 (MQ=255)

                                                                                                                                                cATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACATAATCCAGTATTcac    <  1:2676864/150‑1 (MQ=255)

TCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACATAATCCAGTATTCACAG  >  CP009273/2236951‑2237246

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