BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,722,918	T→A	N247Y (AAC→TAC)	b1647 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,722,918	0	T	A	99.6%	1751.5 / NA	468	N247Y (AAC→TAC)	b1647	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/0); major base A (226/239); minor base C (0/2); total (227/241)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 505 total aligned reads displayed.

GTGTTTCTGCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTTAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACAGAC  >  NC_000913/1722825‑1723010

 tgtTTCTGCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTCAac                                                                                           >  1:9562630/1‑96 (MQ=255)

 tgtTTCTGCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAac                                                                                           >  2:3418228/1‑96 (MQ=255)

 tatttCTGCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAac                                                                                           >  1:14300402/3‑96 (MQ=255)

   tttCTGCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGGTGGGAAACGGAATACCCAGGCGTAAacca                                                                                         >  2:5546864/1‑96 (MQ=255)

   tgtCTGCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAacca                                                                                         <  1:10017990/94‑1 (MQ=255)

    ttCTGCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAaccac                                                                                        <  1:12280486/96‑1 (MQ=255)

     tCTGCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGTGATGGTAAACGTAATACCCAGGCGTAAaccacc                                                                                       <  2:11613305/96‑1 (MQ=255)

     tCTGCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAaccacc                                                                                       >  2:8501854/1‑96 (MQ=255)

      cTGCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCt                                                                                      <  1:10461146/96‑1 (MQ=255)

      cTGCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCt                                                                                      >  2:4009893/1‑96 (MQ=255)

       tGCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTg                                                                                     >  1:428155/1‑96 (MQ=255)

       tGCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTg                                                                                     <  1:5859358/96‑1 (MQ=255)

       tGCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTg                                                                                     >  1:1562503/1‑96 (MQ=255)

         cTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGtt                                                                                   <  2:12384456/96‑1 (MQ=255)

           tCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGGAAACGTAATACCCAGGCGTAAACCACCTGTTCa                                                                                 >  2:6737277/1‑96 (MQ=255)

            cGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCaa                                                                                <  1:4608468/96‑1 (MQ=255)

            cGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCaa                                                                                <  1:9079354/96‑1 (MQ=255)

            cGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCaa                                                                                >  1:13995230/1‑96 (MQ=255)

             gACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAAt                                                                               >  2:11469109/1‑96 (MQ=255)

             gACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAAt                                                                               >  1:8085463/1‑96 (MQ=255)

             gACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGAGGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCaag                                                                               >  1:6142370/1‑95 (MQ=255)

               cAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCg                                                                             <  1:8855891/96‑1 (MQ=255)

                aGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGa                                                                            <  2:2074614/96‑1 (MQ=255)

                aGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGa                                                                            <  1:13539572/96‑1 (MQ=255)

                 gCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAg                                                                           >  2:4504472/1‑96 (MQ=255)

                 gCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAg                                                                           <  1:36800/96‑1 (MQ=255)

                   tGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGcc                                                                         >  2:11315016/1‑96 (MQ=255)

                    gCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCTGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGccc                                                                        <  2:13150684/96‑1 (MQ=255)

                     cccGTACGCGCTCAATTTTTCCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCcg                                                                       >  2:12111968/1‑96 (MQ=255)

                     cccGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCcg                                                                       >  2:13243107/1‑96 (MQ=255)

                      ccGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCcgc                                                                      <  2:9139832/96‑1 (MQ=255)

                       cGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCcgcg                                                                     >  2:5148005/1‑96 (MQ=255)

                        gTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGt                                                                    >  2:5287539/1‑96 (MQ=255)

                         tACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGtt                                                                   >  1:3243440/1‑96 (MQ=255)

                         tACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGtt                                                                   >  1:1206470/1‑96 (MQ=255)

                         tACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGtg                                                                   >  1:5574112/1‑95 (MQ=255)

                          aCGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGttt                                                                  >  1:11552325/1‑96 (MQ=255)

                           cgcgCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTa                                                                 >  1:7821491/1‑96 (MQ=255)

                           cgcgCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTa                                                                 <  2:4228909/96‑1 (MQ=255)

                           cgcgCTCAACTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTa                                                                 <  2:1516458/96‑1 (MQ=255)

                              gcTCAATTTTACATGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCCCCTGTTCAATCGAGCCCGCGTTTAAct                                                              <  1:4081115/96‑1 (MQ=255)

                               cTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAActc                                                             <  2:2348255/96‑1 (MQ=255)

                                tCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAActcc                                                            >  1:7401799/1‑96 (MQ=255)

                                     tttACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGc                                                       <  2:1801682/96‑1 (MQ=255)

                                      ttACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGcc                                                      <  2:6435962/96‑1 (MQ=255)

                                       tACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGcca                                                     <  1:8520920/96‑1 (MQ=255)

                                         ccTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGccacc                                                   >  2:2955763/1‑96 (MQ=255)

                                          cTGAACCGATAGTTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGccacca                                                  <  1:6585461/96‑1 (MQ=255)

                                          cTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGccacca                                                  >  1:9360430/1‑96 (MQ=255)

                                            gAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACcaca                                                >  2:5034689/1‑96 (MQ=255)

                                             aaCCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAg                                               >  1:6278321/1‑96 (MQ=255)

                                               ccGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGcc                                             <  2:8719899/96‑1 (MQ=255)

                                               ccGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACAACAGcc                                             >  2:12947573/1‑96 (MQ=255)

                                                 gATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAg                                           <  1:12966283/96‑1 (MQ=255)

                                                  aTAATTGGCAGCGGCTGCGATAGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGt                                          >  2:5378265/1‑96 (MQ=255)

                                                   tAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGATTAACTCCTCTGCCACGACAGCCAggt                                         >  1:3531741/1‑94 (MQ=255)

                                                    aaTTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTtc                                        >  2:7711451/1‑96 (MQ=255)

                                                     aTTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTtca                                       <  2:789043/96‑1 (MQ=255)

                                                      ttGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTtcat                                      >  2:3325089/1‑96 (MQ=255)

                                                       tGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTtcatc                                     >  1:2352753/1‑96 (MQ=255)

                                                       tGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTtcatc                                     >  1:13235120/1‑96 (MQ=255)

                                                       tGGCAGCGGCTGCGATGGTAAACGTAATACCCACGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTtcatc                                     >  2:14052363/1‑96 (MQ=255)

                                                        ggCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTtcatca                                    <  1:9938580/96‑1 (MQ=255)

                                                         gCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCAc                                   <  1:1073304/96‑1 (MQ=255)

                                                           aGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGGTTAACTCCTCTGCCACCACAGCCAGTTCATCACGc                                 >  1:1129108/1‑96 (MQ=255)

                                                            gCGGCTGCGATGGTAAACGTAATGCCCAGGCGTAAACCACCTGTTCAACCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCa                                <  2:1027224/96‑1 (MQ=255)

                                                            gCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACgga                                <  1:11809334/96‑3 (MQ=255)

                                                             cGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAg                               <  1:2906901/96‑1 (MQ=255)

                                                               gCTGCGATGGTAAACGTAATACCCAGGCGTAACCCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCg                             <  1:2730301/96‑1 (MQ=255)

                                                               gCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCg                             <  2:12512223/96‑1 (MQ=255)

                                                                cTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCgg                            >  1:5480199/1‑96 (MQ=255)

                                                                 tGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGc                           >  1:5207100/1‑96 (MQ=255)

                                                                  gCGATGGTAACCGTAATACCCAGGCGTAACCCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCt                          <  2:7408833/96‑1 (MQ=255)

                                                                  gCGATGGTAACCGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCt                          <  1:4981042/96‑1 (MQ=255)

                                                                  gCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCt                          >  2:11194609/1‑96 (MQ=255)

                                                                  gCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCt                          <  2:71919/96‑1 (MQ=255)

                                                                   cGATGGTAAACGTAATGCCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTg                         <  1:12790416/96‑1 (MQ=255)

                                                                    gATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTgg                        >  2:10755623/1‑96 (MQ=255)

                                                                    gATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTgg                        >  2:10063193/1‑96 (MQ=255)

                                                                    gATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTgg                        >  1:1874549/1‑96 (MQ=255)

                                                                     aTGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGa                       <  1:10966431/96‑1 (MQ=255)

                                                                      tGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGaa                      <  1:2823757/96‑1 (MQ=255)

                                                                        atAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAaat                    >  1:13788596/2‑96 (MQ=255)

                                                                         tAAACGTAATCCCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAaata                   <  2:5961863/96‑1 (MQ=255)

                                                                          aaaCGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAaataa                  >  2:13896538/1‑96 (MQ=255)

                                                                          aaaCGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTACATCACGCAGCGGCTGGAaataa                  >  1:12531694/1‑96 (MQ=255)

                                                                              gTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAatcat              <  2:3819200/96‑1 (MQ=255)

                                                                                 aTACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATcatcat           >  1:3877770/1‑96 (MQ=255)

                                                                                 aTACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATcatcat           <  2:3224639/96‑1 (MQ=255)

                                                                                   aCCCAGGCGTCAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTa         <  1:4866789/96‑1 (MQ=255)

                                                                                     ccAGGCGTAACCCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCCCAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTaaa       <  2:1048891/96‑1 (MQ=255)

                                                                                     ccAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTaaa       <  2:9376904/96‑1 (MQ=255)

                                                                                      cAGGCGTAAACCCCCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAAc      <  2:9562630/96‑1 (MQ=255)

                                                                                      cAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAAc      <  2:3715807/96‑1 (MQ=255)

                                                                                       aGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACa     <  2:7106486/96‑1 (MQ=255)

                                                                                       aGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACa     >  2:264648/1‑96 (MQ=255)

                                                                                        ggCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACAg    >  2:8176232/1‑96 (MQ=255)

                                                                                         gCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACAga   <  2:5835914/96‑1 (MQ=255)

                                                                                         gCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACAga   <  1:7022501/96‑1 (MQ=255)

                                                                                         gCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACAga   <  2:14300402/96‑1 (MQ=255)

GTGTTTCTGCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTTAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACAGAC  >  NC_000913/1722825‑1723010

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