BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	867,645	(C)_8→7	coding (870/1839 nt)	b0829 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	867,638	0	C	.	96.5%	483.7 / 9.4	117	coding (863/1839 nt)	b0829	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (3/1); new base . (66/46); total (70/47)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.70e-01

Only 100 of 144 total aligned reads displayed.

GTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGG  >  NC_000913/867546‑867725

gTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGc                                                                                         >  2:269606/1‑93 (MQ=255)

gTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGc                                                                                         >  2:2194808/1‑93 (MQ=255)

 ttAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGcc                                                                                        >  1:3460264/1‑93 (MQ=255)

 ttAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGcc                                                                                        <  2:1787431/93‑1 (MQ=255)

  tAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGATTCCCGTTGATATCGCTTGAAGATCCAGCGCAACAGGccc                                                                                       >  2:4262426/1‑93 (MQ=255)

   aGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGcccc                                                                                      >  1:2679873/1‑93 (MQ=255)

     ctgctgATCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAAGcccccc                                                                                    >  2:4746538/1‑93 (MQ=255)

      tgctgTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGccccccc                                                                                   <  2:3450914/93‑1 (MQ=255)

      tgctgTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGccccccc                                                                                   >  1:2588611/1‑93 (MQ=255)

      tgctgTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGccccccc                                                                                   >  1:1545671/1‑93 (MQ=255)

       ggtgTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGccccccca                                                                                  <  2:1978833/91‑2 (MQ=255)

       gctgTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGccccccca                                                                                  <  2:3551968/93‑2 (MQ=255)

       gctgTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGccccccat                                                                                  <  2:3361513/93‑3 (MQ=255)

         tgTTCCGCAACTTGGTGCGATGCAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGcccccccatc                                                                                >  2:2992440/1‑90 (MQ=255)

         tgTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGcccccccatc                                                                                <  1:2097380/93‑4 (MQ=255)

          gTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGcccccccatcg                                                                               >  2:3864907/1‑89 (MQ=255)

          gTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGcccccccatcg                                                                               >  2:4374907/1‑89 (MQ=255)

           ttCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGa                                                                             <  1:1172689/93‑1 (MQ=255)

              cGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCa                                                                          >  2:367176/1‑93 (MQ=255)

               gCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAg                                                                         >  1:2749471/1‑93 (MQ=255)

                  aCTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGaaa                                                                      >  1:1196597/1‑93 (MQ=255)

                  aCTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGaaa                                                                      >  1:3054591/1‑93 (MQ=255)

                  aCTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGaaa                                                                      >  2:4527452/1‑93 (MQ=255)

                  aCTTGGGGGGATGAAAGGGTTAGATTAGCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGaaa                                                                      <  1:931319/93‑1 (MQ=255)

                   cTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGaaaa                                                                     >  2:4090690/1‑93 (MQ=255)

                    ttGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAAc                                                                    <  1:579299/93‑1 (MQ=255)

                    ttGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAAc                                                                    >  2:263060/1‑93 (MQ=255)

                      ggTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACgg                                                                  <  1:579895/93‑1 (MQ=255)

                      ggTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACgg                                                                  >  1:4740909/1‑93 (MQ=255)

                       gTGCGAGTAAAGGCTTAGATTATCCCCGGCGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACggt                                                                 <  1:339843/93‑1 (MQ=255)

                         gCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACggtgg                                                               >  1:1876571/1‑93 (MQ=255)

                           gATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGgtggtg                                                             >  1:123080/1‑93 (MQ=255)

                           gATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGgtggtg                                                             >  1:271929/1‑93 (MQ=255)

                           gATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGgtggtg                                                             >  2:1214273/1‑93 (MQ=255)

                           gATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGgtggtg                                                             >  1:1345725/1‑93 (MQ=255)

                             tGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGCGCAGAAAACGGTGGTGGa                                                           >  2:4143420/1‑93 (MQ=255)

                             tGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGa                                                           >  2:1442553/1‑93 (MQ=255)

                             tGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGCGGTGGa                                                           >  1:2382695/1‑93 (MQ=255)

                             tGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑ACCCCCCATCGAGCAGAAAACGGGGGTGGa                                                           >  2:2011276/1‑93 (MQ=255)

                               aaaGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAATATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATg                                                         >  1:3693826/1‑93 (MQ=255)

                               aaaGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATg                                                         <  2:3976886/93‑1 (MQ=255)

                               aaaGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATg                                                         <  2:3112083/93‑1 (MQ=255)

                                aaGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATgg                                                        <  1:2563032/93‑1 (MQ=255)

                                 aGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCCGAAAACGGTGGTGGATggg                                                       >  2:4376404/1‑92 (MQ=255)

                                 aGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGc                                                       >  1:4578834/1‑93 (MQ=255)

                                 aGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGc                                                       <  2:547655/93‑1 (MQ=255)

                                      tAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAAcc                                                  <  2:3686948/93‑1 (MQ=255)

                                      tAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAAcc                                                  <  1:3880569/93‑1 (MQ=255)

                                       aGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGTAGAAAACGGTGGTGGATGGCGAACCt                                                 >  1:4714343/1‑93 (MQ=255)

                                        gATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTg                                                >  2:814305/1‑93 (MQ=255)

                                         aTTATCCCCGACGTTTCCCGTTTATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGt                                               <  1:538779/93‑1 (MQ=255)

                                         aTTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGt                                               <  2:1369726/93‑1 (MQ=255)

                                          ttATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGtt                                              <  1:1418902/93‑1 (MQ=255)

                                          ttATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAAACTGtt                                              >  2:3755904/1‑93 (MQ=255)

                                           tATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGAGGTGGATGGCGAACCTGttt                                             >  1:3927023/1‑93 (MQ=255)

                                               cccGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGAGGATGGCGAACCCGttttgag                                         >  2:254690/1‑90 (MQ=255)

                                                ccGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGa                                        >  2:3956943/1‑93 (MQ=255)

                                                  gACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAgg                                      >  1:3937914/1‑92 (MQ=255)

                                                   aCGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACggagg                                     >  2:4070268/1‑89 (MQ=255)

                                                   aCGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTg                                     >  1:1607149/1‑93 (MQ=255)

                                                   aCGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACA‑AACCCCCCCATCGAGCAGCAAACGGTGGTGGACGGCGAACCTGTTTTACGAgac                                     >  1:2348100/1‑91 (MQ=255)

                                                    cGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGc                                    <  2:1769654/93‑1 (MQ=255)

                                                     gTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCg                                   >  1:3344114/1‑93 (MQ=255)

                                                      tttCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGt                                  >  1:3327415/1‑93 (MQ=255)

                                                        tCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGTGTaa                                >  2:3917451/1‑93 (MQ=255)

                                                         cccGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAAt                               >  2:777429/1‑93 (MQ=255)

                                                          ccGTTGATATCGCTTGAACATCCAGCGAACCAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATc                              <  2:3189615/93‑1 (MQ=255)

                                                          ccGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATc                              >  1:2672083/1‑93 (MQ=255)

                                                             ttGATATCGCTTGAAGATCCAGCGAAACAGG‑CCCCTCCTTCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTGACGAGTGCGTAATCTTg                           <  2:4410489/93‑1 (MQ=255)

                                                              tGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGt                          >  1:4136689/1‑93 (MQ=255)

                                                                atatCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCa                        >  2:190716/1‑93 (MQ=255)

                                                                atatCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCa                        <  1:2909379/93‑1 (MQ=255)

                                                                  atCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCAcc                      <  2:711023/93‑1 (MQ=255)

                                                                  atCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATCGCGAACCTGTTTTACGAGTGCGTAATCTTGTCAcc                      <  2:4695420/93‑1 (MQ=255)

                                                                   tCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCAccc                     >  2:1716639/1‑93 (MQ=255)

                                                                   tCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCAccc                     >  2:110033/1‑93 (MQ=255)

                                                                     gCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGt                   <  1:257272/93‑1 (MQ=255)

                                                                      cTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGtt                  >  1:1735404/1‑93 (MQ=255)

                                                                      cTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGGGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCAGCCGtt                  >  2:1311755/1‑93 (MQ=255)

                                                                       ttGAACATCCAGCGAAACAGGCCCCCCCCATCGACCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGtt                  <  1:4034668/93‑1 (MQ=255)

                                                                       ttGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGttt                 >  2:3200362/1‑93 (MQ=255)

                                                                       ttGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGttt                 >  1:693467/1‑93 (MQ=255)

                                                                       ttGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGGTTTACGAGTGCGTAATCTTGTCACCCGttt                 >  2:154405/1‑93 (MQ=255)

                                                                        tGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTc                <  1:2390511/93‑1 (MQ=255)

                                                                        tGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCAGTTTc                >  2:4260431/1‑93 (MQ=255)

                                                                            cATCCAGCGAATCAGGCCCCCCCCACCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCTCCCGTTTCCCt             <  2:4708215/93‑1 (MQ=255)

                                                                            cATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCtt            <  1:2419946/93‑1 (MQ=255)

                                                                             aTCCAGCGAAACAGG‑CCCTCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCttt           >  2:1680450/1‑93 (MQ=255)

                                                                             aTCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCttt           <  1:460694/93‑1 (MQ=255)

                                                                             aTCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCttt           <  2:4412157/93‑1 (MQ=255)

                                                                              tCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTTCCTTTg          >  2:3042847/1‑93 (MQ=255)

                                                                               ccAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTgc         <  2:1937386/93‑1 (MQ=255)

                                                                                cAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTgcg        >  2:1395092/1‑93 (MQ=255)

                                                                                 aGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTgcgc       >  2:901403/1‑93 (MQ=255)

                                                                                  gCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCa      <  1:3441526/93‑1 (MQ=255)

                                                                                     aaaCAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCg   >  2:3381881/1‑93 (MQ=255)

                                                                                     aaaCAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCg   >  1:1155957/1‑93 (MQ=255)

                                                                                      aaCAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCgg  >  2:260375/1‑93 (MQ=255)

                                                                                      aaCAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCgg  <  1:369071/93‑1 (MQ=255)

                                                                                      aaCAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCgg  >  1:1331862/1‑93 (MQ=255)

GTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGG  >  NC_000913/867546‑867725

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