BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	986,409:1	+A	79.6%	coding (574/1089 nt)	ompF ←	outer membrane porin F

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	986,409	1	.	A	79.6%	319.4 / 71.9	93	coding (574/1089 nt)	ompF	outer membrane porin F
Reads supporting (aligned to +/- strand): ref base . (5/14); new base A (38/36); total (43/50)
Fisher's exact test for biased strand distribution p-value = 7.09e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.44e-01

AGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTT‑AGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGG  >  NC_000913/986311‑986504

aGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGtt                                                                                                  >  1:115942/1‑99 (MQ=255)

  cttcttGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTT‑aa                                                                                               <  2:939264/99‑2 (MQ=255)

    tcttGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTT‑AGaa                                                                                             <  2:2065859/99‑1 (MQ=255)

     cttGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTT‑aagaa                                                                                            <  1:202281/99‑5 (MQ=255)

     cttGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTT‑aagaa                                                                                            <  1:302985/99‑5 (MQ=255)

        gCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACgg                                                                                          >  1:1749241/1‑99 (MQ=255)

         cAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGc                                                                                         <  1:1717402/99‑1 (MQ=255)

          aGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCg                                                                                        >  2:1126737/1‑99 (MQ=255)

          aGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCg                                                                                        <  1:1083396/99‑1 (MQ=255)

           ggTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCCACACCGTCGCCGTTAAGAACGGCGt                                                                                       >  2:302712/1‑99 (MQ=255)

            gTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTg                                                                                      >  2:1507120/1‑99 (MQ=255)

             ttGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGc                                                                                     >  1:800870/1‑99 (MQ=255)

             ttGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGc                                                                                     >  2:1441575/1‑99 (MQ=255)

               ggTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAg                                                                                   <  2:2084032/99‑1 (MQ=255)

                gTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAgt                                                                                  <  2:1357193/99‑1 (MQ=255)

                gTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAgt                                                                                  >  1:377572/1‑99 (MQ=255)

                gTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAgt                                                                                  >  2:1256809/1‑99 (MQ=255)

                gTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTT‑AGAACGGCGTGCAgtg                                                                                 <  2:2347776/99‑1 (MQ=255)

                 tACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAgtg                                                                                 >  1:1671155/1‑99 (MQ=255)

                   cGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTc                                                                               >  1:271703/1‑99 (MQ=255)

                    ggTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCa                                                                              <  2:428888/99‑1 (MQ=255)

                    ggTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCa                                                                              <  2:765065/99‑1 (MQ=255)

                     gTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCAc                                                                             >  1:1258611/1‑99 (MQ=255)

                     gTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCAc                                                                             <  1:1921172/99‑1 (MQ=255)

                           tGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGt                                                                       <  2:980491/99‑1 (MQ=255)

                           tGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGt                                                                       >  2:825685/1‑99 (MQ=255)

                            gCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGtt                                                                      <  1:1337898/99‑1 (MQ=255)

                             cACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGttt                                                                     <  1:387956/99‑1 (MQ=255)

                              aCCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGtttt                                                                    >  1:2235438/1‑99 (MQ=255)

                               ccATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTT‑AGAACGGCGTGCAGTGTCACGCTCGTTTTTa                                                                  <  2:2200252/99‑1 (MQ=255)

                                cATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTa                                                                  <  2:40736/99‑1 (MQ=255)

                                cATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTT‑AGAACGGCGTGCAGTGTCACGCTCGTTTTTAc                                                                 <  2:386365/99‑1 (MQ=255)

                                 aTAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTAc                                                                 <  1:611186/99‑1 (MQ=255)

                                 aTAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTAc                                                                 >  2:895297/1‑99 (MQ=255)

                                 aTAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTT‑AGAACGGCGTGCAGTGTCACGCTCGTTTTTAcc                                                                <  1:1079802/99‑1 (MQ=255)

                                    aGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCa                                                              >  1:796895/1‑99 (MQ=255)

                                    aGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCa                                                              >  2:2311797/1‑99 (MQ=255)

                                     gCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTT‑AGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAgg                                                            <  2:2182243/99‑1 (MQ=255)

                                       aCCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTT‑AGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTa                                                          <  2:2317694/99‑1 (MQ=255)

                                           aCGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTg                                                       <  1:1532705/99‑1 (MQ=255)

                                           aCGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTT‑AGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGa                                                      >  1:127798/1‑99 (MQ=255)

                                             gATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGaa                                                     <  2:1493370/99‑1 (MQ=255)

                                              aTACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAAc                                                    <  1:37499/99‑1 (MQ=255)

                                              aTACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAAc                                                    <  1:890812/99‑1 (MQ=255)

                                              aTACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAAc                                                    <  2:2130826/99‑1 (MQ=255)

                                               tACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACa                                                   >  1:1032128/1‑99 (MQ=255)

                                               tACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACa                                                   >  1:1951726/1‑99 (MQ=255)

                                                aCCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAg                                                  >  2:686457/1‑99 (MQ=255)

                                                  cAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCg                                                >  2:438863/1‑99 (MQ=255)

                                                   aaaGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGa                                               <  2:1161020/99‑1 (MQ=255)

                                                    aaGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTT‑AGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAg                                             >  1:1637091/1‑99 (MQ=255)

                                                     aGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAg                                             <  1:1340998/99‑1 (MQ=255)

                                                     aGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAg                                             <  2:1783759/99‑1 (MQ=255)

                                                      gCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGt                                            >  1:1126349/1‑99 (MQ=255)

                                                       ccTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGtt                                           >  2:322919/1‑99 (MQ=255)

                                                       ccTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGtt                                           <  1:876374/99‑1 (MQ=255)

                                                        cTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTc                                          <  1:1209487/99‑1 (MQ=255)

                                                         ttcgtattcgtaGCTGATAGAACCGCCAACACCGTCGCCGTT‑AGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAg                                        <  2:550840/99‑1 (MQ=255)

                                                          tcgtattcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAg                                        <  1:410941/99‑1 (MQ=255)

                                                            gtattcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGc                                      <  1:167442/99‑1 (MQ=255)

                                                            gtattcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGc                                      <  1:1301745/99‑1 (MQ=255)

                                                             tattcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGcc                                     >  1:1070656/1‑99 (MQ=255)

                                                             tattcgtaGCTGATAGAACCGCCAACACCGTCGCCGTT‑AGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCa                                    <  2:1919385/99‑1 (MQ=255)

                                                             tattcgtaGCTGATAGAACCGCCAACACCGTCGCCGTT‑AGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCa                                    <  2:537442/99‑1 (MQ=255)

                                                               ttcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCAt                                   >  1:2007565/1‑99 (MQ=255)

                                                               ttcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCAt                                   >  1:816909/1‑99 (MQ=255)

                                                               ttcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCAt                                   <  2:2182/99‑1 (MQ=255)

                                                               ttcgtaGCTGATAGAACCGCCAACACCGTCGCCGTT‑AGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATc                                  >  2:2414019/1‑99 (MQ=255)

                                                                  gtaGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCaa                                >  1:1932593/1‑99 (MQ=255)

                                                                  gtaGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCaa                                >  1:2215579/1‑99 (MQ=255)

                                                                   taGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAAc                               >  1:736970/1‑99 (MQ=255)

                                                                    aGCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAAcc                              <  1:1221409/99‑1 (MQ=255)

                                                                     gCTGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCa                             >  1:2103146/1‑99 (MQ=255)

                                                                       tGATAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGa                           <  2:1973969/99‑1 (MQ=255)

                                                                       tGATAGAACCGCCAACACCGTCGCCGTT‑AGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGAc                          <  1:1397020/99‑1 (MQ=255)

                                                                         aTAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGAcc                         <  1:1777744/99‑1 (MQ=255)

                                                                         aTAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGAcc                         >  2:154909/1‑99 (MQ=255)

                                                                         aTAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGAcc                         <  2:468417/99‑1 (MQ=255)

                                                                         aTAGAACCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGAcc                         >  2:1158267/1‑99 (MQ=255)

                                                                          tAGAACCGCCAACACCGTCGCCGTT‑AGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTTCTGCACAGCGCATTACAGCCCAACCACCAGACTaa                       >  2:1947597/1‑99 (MQ=255)

                                                                              aCCGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAaaga                    <  1:1890854/99‑1 (MQ=255)

                                                                               ccGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAaagaa                   >  1:1017861/1‑99 (MQ=255)

                                                                               ccGCCAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAaagaa                   >  1:1277321/1‑99 (MQ=255)

                                                                                  ccAACACCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGtt                >  2:645461/1‑99 (MQ=255)

                                                                                    aaCACCGTCGCCGTT‑AGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGa             <  2:497959/99‑1 (MQ=255)

                                                                                     acacCGTCGCCGTT‑AGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAg            <  2:2096942/99‑1 (MQ=255)

                                                                                      cacCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAg            <  2:2215579/99‑1 (MQ=255)

                                                                                      cacCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAg            <  1:477634/99‑1 (MQ=255)

                                                                                      cacCGTCGCCGTT‑AGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGt           >  1:553421/1‑99 (MQ=255)

                                                                                       acCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGt           <  2:27627/99‑1 (MQ=255)

                                                                                       acCGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGt           >  1:735903/1‑99 (MQ=255)

                                                                                        ccGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGtt          >  2:4489/1‑99 (MQ=255)

                                                                                         cGTCGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTa         <  2:924398/99‑1 (MQ=255)

                                                                                            cGCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGa      <  1:367284/99‑1 (MQ=255)

                                                                                             gCCGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGAt     >  1:1237264/1‑99 (MQ=255)

                                                                                              ccGTTAAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATa    >  1:1402394/1‑99 (MQ=255)

                                                                                               cGTT‑AGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAgg  <  2:1246639/99‑1 (MQ=255)

AGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTT‑AGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGG  >  NC_000913/986311‑986504

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