BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,312,170	G→C	6.1%	L194V (CTG→GTG)	ompC ←	outer membrane porin C

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,312,170	0	G	C	6.1%	557.6 / 18.4	164	L194V (CTG→GTG)	ompC	outer membrane porin C
Reads supporting (aligned to +/- strand): ref base G (99/55); new base C (7/3); total (106/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

Only 100 of 167 total aligned reads displayed.

TGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCT  >  NC_000913/2312023‑2312319

  tAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTAGGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCg                                                                                                                                                   <  2:165488/150‑1 (MQ=255)

    ggTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTc                                                                                                                                                 >  2:2450143/1‑150 (MQ=255)

    ggTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTc                                                                                                                                                 >  2:2418525/1‑150 (MQ=255)

    ggTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTc                                                                                                                                                 >  2:1971263/1‑150 (MQ=255)

      tttCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCAc                                                                                                                                               <  1:399480/150‑1 (MQ=255)

         cAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGAc                                                                                                                                            >  2:1812943/1‑150 (MQ=255)

         cAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGAc                                                                                                                                            >  2:1393848/1‑150 (MQ=255)

         cAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGAc                                                                                                                                            >  1:730747/1‑150 (MQ=255)

          aGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCACTGCGTCaaaacg                                                                                                                                           <  1:1565349/150‑6 (MQ=255)

           gCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCg                                                                                                                                          >  2:2066785/1‑150 (MQ=255)

           gCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCACTGCGTCaaaacgg                                                                                                                                          <  1:2407061/150‑7 (MQ=255)

            cACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCgt                                                                                                                                         <  2:1375875/150‑1 (MQ=255)

              cGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCgttg                                                                                                                                       >  2:288129/1‑150 (MQ=255)

                gTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCgttgtt                                                                                                                                     >  2:1182902/1‑150 (MQ=255)

                    ccGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTgttagtt                                                                                                                                 >  1:2161763/1‑150 (MQ=255)

                     cGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTgttagtta                                                                                                                                >  2:522276/1‑150 (MQ=255)

                     cGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTgttagtta                                                                                                                                >  2:1232523/1‑150 (MQ=255)

                      gTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTAc                                                                                                                               >  1:2386762/1‑150 (MQ=255)

                          ccGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCa                                                                                                                           >  2:1313122/1‑150 (MQ=255)

                            gATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACt                                                                                                                         >  1:397009/1‑150 (MQ=255)

                                tAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTa                                                                                                                     <  1:1113280/150‑1 (MQ=255)

                                 aaGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTgttagt                                                                                                                                  <  1:945287/136‑1 (MQ=255)

                                 aaGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTAAGTTACGCCACTAGTaa                                                                                                                    >  1:141755/1‑150 (MQ=255)

                                   gcagcGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAg                                                                                                                  >  2:562257/1‑150 (MQ=255)

                                    cagcGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGc                                                                                                                 >  2:2141767/1‑150 (MQ=255)

                                    cagcGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGc                                                                                                                 >  1:578670/1‑150 (MQ=255)

                                         gtgtTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACGGTTGTTAGTTACGCCACTAGTAAAGGCTTCa                                                                                                            <  2:1440391/150‑1 (MQ=255)

                                         gtgtTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCa                                                                                                            >  2:244135/1‑150 (MQ=255)

                                          tgtTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCAc                                                                                                           <  1:1343112/150‑1 (MQ=255)

                                             tCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAg                                                                                                        >  2:1374027/1‑150 (MQ=255)

                                             tCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCCTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAg                                                                                                        >  1:2298463/1‑150 (MQ=255)

                                              cTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGa                                                                                                       >  1:1476065/1‑150 (MQ=255)

                                                tagaatgaGTACGTTTGGAGCTGTAGATCGCACCACTGATTCCGAAACCTTCGTAATCATAAGTGATATAACCGCCGACGCCGTAGGCGTTTTGACGCAGTGCGTCAAGAGCTTTGTTAGTTACGCCACTAGTAAAGCCTTCACAAGATg                                                                                                     <  2:2249805/143‑1 (MQ=255)

                                                   cATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGt                                                                                                  >  2:791766/1‑150 (MQ=255)

                                                    aTCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAg                                                                                                                                                       <  1:120520/96‑1 (MQ=255)

                                                    aTCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAg                                                                                                                                                       >  2:120520/1‑96 (MQ=255)

                                                           cGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGttt                                                                                          >  1:662763/1‑150 (MQ=255)

                                                            gTTTGGAGCTGGAGATCGCACCACCGATACCTAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGACGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACTCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGtttt                                                                                         <  2:565071/150‑1 (MQ=255)

                                                            gTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGtttt                                                                                         >  2:885162/1‑150 (MQ=255)

                                                               tGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTAc                                                                                      >  2:55240/1‑150 (MQ=255)

                                                               tGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTAc                                                                                      >  2:97117/1‑150 (MQ=255)

                                                               tGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTAc                                                                                      <  2:1097286/150‑1 (MQ=255)

                                                                 gAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTAccc                                                                                    >  1:509847/1‑150 (MQ=255)

                                                                   gCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTg                                                                                  >  1:1313878/1‑150 (MQ=255)

                                                                      ggAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTa                                                                               >  1:2232439/1‑150 (MQ=255)

                                                                         gaTCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTg                                                                            >  2:1512226/1‑150 (MQ=255)

                                                                            cGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAAc                                                                         >  2:2471039/1‑150 (MQ=255)

                                                                               acgaccGATACCTTAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCTGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATTGGTTGCCGTTTTTACCCTGGTACTGAACAGc                                                                      <  2:509847/147‑1 (MQ=255)

                                                                                 caccGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCaa                                                                    <  2:671426/150‑1 (MQ=255)

                                                                                 caccGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCaa                                                                    >  1:476768/1‑150 (MQ=255)

                                                                                  accGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCaaa                                                                   >  1:643410/1‑150 (MQ=255)

                                                                                    cGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGt                                                                 >  2:1849500/1‑150 (MQ=255)

                                                                                     gATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGtt                                                                >  2:670687/1‑150 (MQ=255)

                                                                                     gATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGtt                                                                >  2:1750768/1‑150 (MQ=255)

                                                                                         ccGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAgg                                                            >  2:440950/1‑150 (MQ=255)

                                                                                         ccGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAgg                                                            >  2:1702368/1‑150 (MQ=255)

                                                                                         ccGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCATGGTACTGAACAGCAAAGTTCAgg                                                            >  2:1867638/1‑150 (MQ=255)

                                                                                           gAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGcc                                                          <  1:791766/150‑1 (MQ=255)

                                                                                               ccTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCa                                                      >  2:2108837/1‑150 (MQ=255)

                                                                                               ccTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCa                                                      >  2:2220656/1‑150 (MQ=255)

                                                                                                  tCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAAcc                                                   <  1:2458549/150‑1 (MQ=255)

                                                                                                  tCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAAcc                                                   <  1:2450143/150‑1 (MQ=255)

                                                                                                  tCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAAcc                                                   <  2:1004641/150‑1 (MQ=255)

                                                                                                   cGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCa                                                  <  1:1374027/150‑1 (MQ=255)

                                                                                                    gTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAg                                                 <  2:43186/150‑1 (MQ=255)

                                                                                                      aaTCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGAc                                               <  2:731265/150‑1 (MQ=255)

                                                                                                           tAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCgaag                                          >  1:1971818/1‑150 (MQ=255)

                                                                                                            aaGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCgaaga                                         >  1:1164159/1‑150 (MQ=255)

                                                                                                             aGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCgaagaa                                        >  2:1559117/1‑150 (MQ=255)

                                                                                                               tGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGt                                      <  1:440950/150‑1 (MQ=255)

                                                                                                                gATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTc                                     <  1:1512226/150‑1 (MQ=255)

                                                                                                                 aTAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCa                                    >  1:296990/1‑150 (MQ=255)

                                                                                                                  tAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAg                                   <  1:679224/150‑1 (MQ=255)

                                                                                                                    gAACCGCCGACGCCGTAGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAgtg                                 >  2:464208/1‑150 (MQ=255)

                                                                                                                     aaCCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCTAAGAAGTCAgtgt                                <  1:2418525/150‑1 (MQ=255)

                                                                                                                      aCCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTCCGCCCCTAGTACAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGtt                               >  1:33074/1‑150 (MQ=255)

                                                                                                                       ccgccgACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGATACGCCACAAGTAAACCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTa                              >  2:2277660/1‑150 (MQ=255)

                                                                                                                        cgccgACGCCGTCGCCGTTTTGACGCAGTGCGTCACGCCCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTAc                             >  1:727071/1‑150 (MQ=255)

                                                                                                                        cgccgACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTAc                             >  1:2293357/1‑150 (MQ=255)

                                                                                                                          ccgACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                           <  1:1010203/150‑1 (MQ=255)

                                                                                                                            gACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTa                         <  1:1387529/150‑1 (MQ=255)

                                                                                                                                ccgtcgccgtTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTc                     >  1:238981/1‑150 (MQ=255)

                                                                                                                                  gtcgccgtTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTcgc                   >  2:1180459/1‑150 (MQ=255)

                                                                                                                                   tcgccgtTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTcgcg                  <  1:1393848/150‑1 (MQ=255)

                                                                                                                                    cgccgtTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGAc                                               >  1:1756930/1‑120 (MQ=255)

                                                                                                                                     acggtcgtGACGCACTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGaa                >  2:2322277/8‑150 (MQ=255)

                                                                                                                                      ccgtTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAg               <  2:643410/150‑1 (MQ=255)

                                                                                                                                        gtcgtGACGCACTGTGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGCTGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGcc             >  2:825013/5‑150 (MQ=255)

                                                                                                                                        gtcgtGACGCACTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGcc             >  2:392753/5‑150 (MQ=255)

                                                                                                                                        gtTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGcc             >  1:1239891/1‑150 (MQ=255)

                                                                                                                                        gtTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGcc             >  1:1854609/1‑150 (MQ=255)

                                                                                                                                           ttGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGtt          <  2:1547601/150‑1 (MQ=255)

                                                                                                                                           ttGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGtt          <  2:314600/150‑1 (MQ=255)

                                                                                                                                            tGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTa         <  2:1437200/150‑1 (MQ=255)

                                                                                                                                            tGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTa         <  1:2046204/150‑1 (MQ=255)

                                                                                                                                              aCGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTAcc       >  1:2125584/1‑150 (MQ=255)

                                                                                                                                              aCGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTAcc       <  2:233883/150‑1 (MQ=255)

                                                                                                                                               cGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCa      >  1:935290/1‑150 (MQ=255)

                                                                                                                                                gCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCAc     <  1:2471039/150‑1 (MQ=255)

                                                                                                                                                  aGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACgc   <  1:387044/150‑1 (MQ=255)

TGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCT  >  NC_000913/2312023‑2312319

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`