BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	158,300	0	A	G	65.4%	267.2 / 302.0	295	F265S (TTC→TCC)	pcnB	poly(A) polymerase I
Reads supporting (aligned to +/- strand): ref base A (50/52); new base G (96/97); total (146/149)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 306 total aligned reads displayed.

GGATACGCGTATCGGTATTCTTCAGCACCTGTTCAATGATCCGCTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGAACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTCGTAACCGTAGCCCGCTTGTAGCAGTTTAAGCGATTCTTCAAACAGG  >  NC_000913/158208‑158392

   tACGCGTATCGGTATTCTTCAGCACCTGTTCAATGATCCGCTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGGACa                                                                                           <  2:1006795/93‑1 (MQ=255)

    aCGCGTATCGGTATTCTTCAGCACCTGTTCAATGATCCGCTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGGACAg                                                                                          >  1:780640/1‑93 (MQ=255)

    aCGCGTATCGGTATTCTTCAGCACCTGTTCAATGATCCGCTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGGACAg                                                                                          >  2:7924688/1‑93 (MQ=255)

     cgcgTATCGGTATTCTTCAGCACCTGTTCAATGATCCGCTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGGACAGc                                                                                         <  1:6827865/93‑1 (MQ=255)

     cgcgTATCGGTATTCTTCAGCACCTGTTCAATGATCCGCTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGGACAGc                                                                                         <  1:1715287/93‑1 (MQ=255)

     cgcgTATCGGTATTCTTCAGCACCTGTTCAATGATCCGCTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGAACAGc                                                                                         <  1:1008307/93‑1 (MQ=255)

     cgcgTATCGGTATTCTTCAGCACCTGTTCAATGATCCGCTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGGAATGGTCGGGGACAGc                                                                                         >  2:7308832/1‑93 (MQ=255)

      gcgTATCGGTATTCTTCAGCACCTGTTCAATGATCCGCTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGAACAGCg                                                                                        <  1:5102769/93‑1 (MQ=255)

        gTATCGGTATTCTTCAGCACCTGTTCAATGATCCGCTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGGACAGCGGc                                                                                      <  2:8736991/93‑1 (MQ=255)

          aTCGGTATTCTTCAGCACCTGTTCAATGATCCGCTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGGACAGCGGCTg                                                                                    <  1:1164717/93‑1 (MQ=255)

          aTCGGTATTCTTCAGCACCTGTTCAATGATCCGCTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGAACAGCGGCTg                                                                                    >  1:6200647/1‑93 (MQ=255)

           tCGGTATTCTTCAGCACCTGTTCAATGATCCGCTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGGACAGCGGCTgg                                                                                   >  1:1545109/1‑93 (MQ=255)

            cGGTATTCTTCAGCACCTGTTCAATGATCCGCTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGGACAGCGGCTGGa                                                                                  >  1:822984/1‑93 (MQ=255)

            cGGTATTCTTCAGCACCTGTTCAATGATCCGCTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGAACAGCGGCTGGa                                                                                  >  2:1887601/1‑93 (MQ=255)

              gTATTCTTCAGCACCTGTTCAATGATCCGCTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGGACAGCGGCTGGAAc                                                                                >  2:3909867/1‑93 (MQ=255)

               tATTCTTCAGCACCTGTTCAATGATCCGCTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGGACAGCGGCTGGAACa                                                                               >  1:7472339/1‑93 (MQ=255)

                aTTCTTCAGCACCTGTTCAATGATCCGCTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGGACAGCGGCTGGAACAg                                                                              >  1:5949896/1‑93 (MQ=255)

                aTTCTTCAGCACCTGTTCAATGATCCGCTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGGACAGCGGCTGGAACAg                                                                              >  1:3062741/1‑93 (MQ=255)

                 ttcttcATCAACTGGTCAATGATCCGCTCCATCGGGCTGTTGCCATTTTCGGTGAAGCAGCGGGTAATGGTCGGGGACAGCGGCTGGGACatg                                                                             >  2:1468613/1‑91 (MQ=255)

                 ttcttcAGCACCTGTTCAATGATCCGCTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGGACAGCGGCTGGAACAga                                                                             <  1:8493052/93‑1 (MQ=255)

                    ctcAGCACCTGGTCAATGATCCGCTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGGACAGCGGCTGGAACAgatga                                                                          <  1:2408792/92‑1 (MQ=255)

                      cAGCACCTGTTCAATGATCCGCTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGGACAGCGGCTGGAACAGATGata                                                                        >  1:2640188/1‑93 (MQ=255)

                      cAGCACCTGTTCAATGATCCGCTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGGACAGCGGCTGGAACAGATGata                                                                        <  1:8807004/93‑1 (MQ=255)

                      cAGCACCTGTTCAATGATCCGCTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGAACAGCGGCTGGAACAgatgatt                                                                        >  1:5488582/1‑92 (MQ=255)

                       aGCACCTGTTCAATGATCCGCTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGGACAGCGGCTGGAACAGATGatat                                                                       <  1:1053459/93‑1 (MQ=255)

                        gCACCTGTTCAATGATCCGCTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGGACAGCGGCTGGAACAGATGATAtt                                                                      >  2:4114248/1‑93 (MQ=255)

                        gCACCTGTTCAATGATCCGCTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGGAATGGGCGGGAACAGCGGCGGGAACAGATGATAtt                                                                      >  1:1600036/1‑93 (MQ=255)

                         cACCTGTTCAATGATCCGCTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGGACAGCGGCTGGAACAGATGATATTc                                                                     <  1:2526556/93‑1 (MQ=255)

                         cACCTGTTCAATGATCCGCTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGAACAGCGGCTGGAACAGATGATATTc                                                                     >  2:3522318/1‑93 (MQ=255)

                          aCCTGTTCAATGATCCGCTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGGACAGCGGCTGGAACAGAAGATATTca                                                                    >  2:6349595/1‑93 (MQ=255)

                           ccTGTTCAATGATCCGCTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGAACAGCGGCTGGAACAGAAGATATTcat                                                                   >  2:1161097/1‑92 (MQ=255)

                            cTGTTCAATGATCCGCTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGAACAGCGGCTGGAACAGATGATATTcaca                                                                  >  1:585324/1‑93 (MQ=255)

                             tGTTCAATGATCCGCTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGAACAGCGGCTGGAACAGATGATATTcacac                                                                 >  1:4726275/1‑93 (MQ=255)

                               ttaaaTGATCCGCGCGATCGGGCTGTCGCCATTTCCCGTGAAGTAGGGGGTAATGGTCGGGAACAGCGGCT‑TAACAGATGATATTCACacaac                                                              <  2:4390953/90‑1 (MQ=255)

                               ttCAATGATCCGCTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGGACAGCGGCTGGAACAGATGATATTCACacaa                                                               >  2:1768936/1‑93 (MQ=255)

                                tCAATGATCCGCTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGGACAGCGGCTGGAACAGATGATATTCACacaac                                                              <  2:1624302/93‑1 (MQ=255)

                                 cAATGATCCGCTCCATCTGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGGACAGCGGCTGGAACAGATGATATTCACacaaca                                                             >  2:7815818/1‑93 (MQ=255)

                                 cAATGATCCGCTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGGACAGCGGCTGGAACAGATGATATTCACacaaca                                                             >  2:3639813/1‑93 (MQ=255)

                                 cAATGATCCGCTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGGACAGCGGCTGGAACAGATGATATTCACacaaca                                                             <  1:2258020/93‑1 (MQ=255)

                                  aaTGATCCGCTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGGACAGCGGCTGGAACAGATGATATTCACACAACAg                                                            >  2:7196805/1‑93 (MQ=255)

                                   attaTCCCCTCCATCGGGCTGTCCCCATTTTCCCTTAAGGAGCGGGGAATGGTTGGGGACAGCGGCTGGAACAGATGATATTCACACAACAGc                                                           <  1:3936547/90‑1 (MQ=255)

                                    tGATCCGCTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGGACAGCGGCTGGAACAGATGATATTCACCCACCAGCt                                                          >  2:2758071/1‑93 (MQ=255)

                                     gATCCGCTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGGACAGCGGCTGGAACAGATGATATTCACACAACAGCtt                                                         >  1:7657115/1‑93 (MQ=255)

                                       tCCGCTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGGACAGCGGCTGGAACAGATGATATTCACACAACAGCTtat                                                       <  1:1760330/93‑1 (MQ=255)

                                        ccGCTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGGACAGCGGCTGGAACAGATGATATTCACACAACAGCTtata                                                      <  1:8622391/93‑1 (MQ=255)

                                        ccGCTCCATCGGGCTGTCGCCATTGTCCGTGAAGTAGCGGGTAATGCTCGGGGACAGCGGCTGGAACAGATGATATTCACACAACAGCTtata                                                      <  1:7194560/93‑1 (MQ=255)

                                         cGCTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGGACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAg                                                     <  1:7938414/93‑1 (MQ=255)

                                           cTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGAACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGt                                                   <  2:8594295/93‑1 (MQ=255)

                                           cTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGGAATGGTCGGGGACAGCGGCTGGAACCGAAGATATTTCCAACACAGCTTATAGGt                                                   >  2:6888343/1‑93 (MQ=255)

                                             ccATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGGACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGttt                                                 <  2:5507409/93‑1 (MQ=255)

                                             ccATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGAACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGttt                                                 >  1:1952246/1‑93 (MQ=255)

                                              cATCGGGCTGTCGCTATTTTCCGTGAAGTAGCGGGTAATGGTCGGGAACAGCGTCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTc                                                <  2:4851794/93‑1 (MQ=255)

                                              cATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGGACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTc                                                >  2:2210714/1‑93 (MQ=255)

                                              cATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGAACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTc                                                <  1:8022838/93‑1 (MQ=255)

                                                  gggCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGGACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTCGTaa                                            >  2:5177017/1‑93 (MQ=255)

                                                  gggCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGAACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTCGTaa                                            >  1:7778535/1‑93 (MQ=255)

                                                    gCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCTGGGACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTCGTaaac                                          >  1:180394/1‑91 (MQ=255)

                                                     gtGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGGACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTCGTAACCg                                         <  2:6774197/92‑1 (MQ=255)

                                                     cTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGGACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTCGTAACCg                                         >  1:5372341/1‑93 (MQ=255)

                                                      tGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGGACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTCGTAACCGt                                        >  2:8148362/1‑93 (MQ=255)

                                                       gTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGGACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTCGTAACCGTa                                       >  2:7678736/1‑93 (MQ=255)

                                                       gTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGAACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTCGTAACCGTa                                       >  2:5977590/1‑93 (MQ=255)

                                                        tCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGGACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTCGTAACCGTAg                                      >  1:6983660/1‑93 (MQ=255)

                                                         cGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGGACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTCGTAACCGTAGc                                     <  2:62095/93‑1 (MQ=255)

                                                          gCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGGACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTCGTAACCGTCGcc                                    >  1:2181335/1‑93 (MQ=255)

                                                           ccATTTTCCGTGAAGTAGCGGGTAATGGTCGGGAACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTCGTAACCGTAGccc                                   >  2:1334319/1‑93 (MQ=255)

                                                            cATTTTCCGTGAAGTAGCGGGTAATGGTCGGGAACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTCGTAACCGTAGCCCg                                  <  2:2989768/93‑1 (MQ=255)

                                                             aTTTTCCGTGAAGTAGCGGGTAATGGTCGGGGACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTCGTAACCGTAGCCCGc                                 >  1:7134564/1‑93 (MQ=255)

                                                              ttttCCGTGAAGTAGCGGGTAATGGTCGGGGACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTCGTAACCGTAGCCCGCt                                <  2:3197538/93‑1 (MQ=255)

                                                              ttttCCGTGAAGTAGCGGGTAATGGTCGGGGACAGCGGCTGGAACAGATGATATTCACACAACAACTTATAGGATTCGTAACCGTACCCCGCt                                >  1:496632/1‑93 (MQ=255)

                                                              ttttCCGTGAAGTAGCGGGTAATGGTCGGGAACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTCGTAACCGTAGCCCGCt                                <  1:3204827/93‑1 (MQ=255)

                                                              ttttCCGTGAAGTAGCGGGTAATGGTCGGGAACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTCGTAACCGTAGCCCGCt                                <  2:4261845/93‑1 (MQ=255)

                                                                 tCCGTGAAGTAGCGGGTAATGGTCGGGAACAGCGGCTGGAACAGATGATATTCACAAAACAGCTTATAGGTTTCGTAACCGTAGCCCGCTTGt                             >  2:8470779/1‑93 (MQ=255)

                                                                   cGTGAAGTAGCGGGTAATGGTCGGGGACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTCGTAACCGTAGCCCGCTTGTAg                           <  1:4433183/93‑1 (MQ=255)

                                                                     tGAAGTAGCGGGTAATGGTCGGGGACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTCGTAACCGTAGCCCGCTTGTAGCa                         <  2:2353356/93‑1 (MQ=255)

                                                                     tGAAGTAGCGGGTAATGGTCGGGAACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTCGTAACCGTAGCCCGCTTGTAGCa                         <  1:8172146/93‑1 (MQ=255)

                                                                      gAAGTAGCGGGTAATGGTCGGGGACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTCGTAACCGTAGCCCGCTTGTAGCAg                        >  1:6494837/1‑93 (MQ=255)

                                                                      gAAGTAGCGGGTAATGGTCGGGAACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTCGTAACCGTAGCCCGCTTGTAGCAg                        >  2:2011298/1‑93 (MQ=255)

                                                                         gTAGCGGGTAATGGTCGGGGACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTCGTAACCGTAGCCCGCTTGTAGCAGttt                     >  2:8312780/1‑93 (MQ=255)

                                                                         gTAGCGGGTAATGGTCGGGAACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTCGTAACCGTAGCCCGCTTGTAGCAGttt                     <  2:7322460/93‑1 (MQ=255)

                                                                          tAGCGGGTAATGGTCGGGGACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTCGTAACCGTAGCCCGCTTGTAGCAGTTTa                    >  1:144896/1‑93 (MQ=255)

                                                                           aGCGGGTAATGGTCGGGGACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTCGTAACCGTAGCCCGCTTGTAGCAGTTTaa                   <  1:348173/93‑1 (MQ=255)

                                                                            gCGGGTAATGGTCGGGGACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTCGTAACCGTAGCCCGCTTGTAGCAGTTTAAg                  <  1:5687721/93‑1 (MQ=255)

                                                                            gCGGGTAATGGTCGGGGACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTCGTAACCGTAGCCCGCTTGTAGCAGTTTAAg                  >  1:521810/1‑93 (MQ=255)

                                                                              gggTAATGGTCGGGAACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTCGTAACCGTAGCCCGCTTGTAGCAGTTTAAGCg                >  1:866538/1‑93 (MQ=255)

                                                                                gTAATGGTCGGGGACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTCGTAACCGTAGCCCGCTTGTAGCAGTTTAAGCGAt              >  2:394136/1‑93 (MQ=255)

                                                                                gTAATGGTCGGGGACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTCGTAACCGTAGCCCGCTTGTAGCAGTTTAAGCGAt              <  2:2899685/93‑1 (MQ=255)

                                                                                 tAATGGTCGGGGACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTCGTAACCGTAGCCCGCTTGTAGCAGTTTAAGCGAtt             >  2:6955679/1‑93 (MQ=255)

                                                                                  aaTGGTCGGGGACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTCGTAACCGTAGCCCGCTTGTAGCAGTTTAAGCGAttc            <  2:6603598/93‑1 (MQ=255)

                                                                                  aaTGGTCGGGAACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTCGTAACCGTAGCCCGCTTGTAGCAGTTTAAGCGAttc            <  1:2031559/93‑1 (MQ=255)

                                                                                   aTGGTCGGGGACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTCGTAACCGTAGCCCGCTTGTAGCAGTTTAAGCGAttct           <  2:8224270/93‑1 (MQ=255)

                                                                                   aTGGTCGGGAACAGCGGCTGGAACAGATTATATTCACACAACAGCTTATAGGTTTCGTAACCGTAGCCCGCTTGTAGCAGTTTAAGCGAttct           <  1:1349230/93‑1 (MQ=255)

                                                                                     ggTCGGGGACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTCGTAACCGTAGCCCGCTTGTAGCAGTTTAAGCGAttcttc         <  2:4922363/93‑1 (MQ=255)

                                                                                     ggTCGGGAACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTCGTAACCGTCGCCCGCTTGTAGCAGTTTAAGCGAttgttc         >  2:3563252/1‑93 (MQ=255)

                                                                                     ggTCGGGAACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTCGTAACCGTAGCCCGCTTGTAGCAGTTTAAGCGAttcttc         >  1:4662385/1‑93 (MQ=255)

                                                                                      gTCGGGGACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTCGTAACCGTAGCCCGCTTGTAGCAGTTTAAGCGATTCTTCa        >  2:4069467/1‑93 (MQ=255)

                                                                                       tCGGGAACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTCGTAACCGTAGCCCGCTTGTAGCAGTTTAAGCGATTCTTCaa       >  2:2597652/1‑93 (MQ=255)

                                                                                       tCGGGAACAGCGGCTGGAACAGATGATATTCACACAACAGATTATAGGTTTCGTAACCGTAGCCCGCTTGTAGCAGTTTAAGCGATTCTTCaa       >  1:8718755/1‑93 (MQ=255)

                                                                                         ggggACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTCGTAACCGTAGCCCGCTTGTAGCAGTTTAAGCGATTCTTCAAAc     <  2:4710123/93‑1 (MQ=255)

                                                                                         gggAACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTCGTAACCGTAGCCCGCTTGTAGCAGTTTAAGCGATTCTTCAAAc     <  1:8282426/93‑1 (MQ=255)

GGATACGCGTATCGGTATTCTTCAGCACCTGTTCAATGATCCGCTCCATCGGGCTGTCGCCATTTTCCGTGAAGTAGCGGGTAATGGTCGGGAACAGCGGCTGGAACAGATGATATTCACACAACAGCTTATAGGTTTCGTAACCGTAGCCCGCTTGTAGCAGTTTAAGCGATTCTTCAAACAGG  >  NC_000913/158208‑158392

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