BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,442,106	A→G	S28S (AGT→AGC)	b2325 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,442,106	0	A	G	99.3%	454.5 / ‑4.8	142	S28S (AGT→AGC)	b2325	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (84/57); total (84/58)
Fisher's exact test for biased strand distribution p-value = 4.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

Only 100 of 146 total aligned reads displayed.

GCTAACGGTCGTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCACTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCATCTTCTT  >  NC_000913/2442016‑2442195

gCTAACGGTCGTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCt                                                                                         <  2:2152461/93‑1 (MQ=255)

gCTAACGGTCGTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCt                                                                                         <  2:3535369/93‑1 (MQ=255)

 cTAACGGTCGTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTg                                                                                        <  1:2383317/93‑1 (MQ=255)

 cTAACGGTCGTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTg                                                                                        >  2:3308608/1‑93 (MQ=255)

 cTAACGGTCGTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTg                                                                                        >  2:3308607/1‑93 (MQ=255)

   aaCGGTCGTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGc                                                                                      >  1:830538/1‑93 (MQ=255)

      ggTCGTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCaac                                                                                   >  1:2620774/1‑92 (MQ=255)

      ggTCGTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCaaa                                                                                   >  2:3583915/1‑93 (MQ=255)

      ggTCGTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCaaa                                                                                   >  1:1267585/1‑93 (MQ=255)

       gTCGTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAaac                                                                                  >  1:2759124/1‑93 (MQ=255)

         cGTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAaacaa                                                                                <  1:3775696/93‑1 (MQ=255)

         cGTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAaacaa                                                                                <  1:378957/93‑1 (MQ=255)

          gTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAaacaac                                                                               >  2:1097477/1‑93 (MQ=255)

             tGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACtca                                                                            <  2:2354596/93‑1 (MQ=255)

             tGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACtca                                                                            >  2:2956768/1‑93 (MQ=255)

               cACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACtcatc                                                                          <  2:3980626/93‑1 (MQ=255)

                 ccGCCTGAGCGGAGTGGTTATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGACACGCAAATANNCGCTGGCAAACAACtcatnat                                                                        >  1:16275/1‑93 (MQ=255)

                 ccGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACtcatcat                                                                        >  1:1331148/1‑93 (MQ=255)

                 ccGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACtcatcat                                                                        >  1:2646433/1‑93 (MQ=255)

                 ccGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACtcatcat                                                                        >  2:1968500/1‑93 (MQ=255)

                 ccGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACtcatcat                                                                        >  2:416859/1‑93 (MQ=255)

                 ccGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACtcatcat                                                                        >  2:416860/1‑93 (MQ=255)

                    ccTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCAc                                                                     >  1:3578137/1‑93 (MQ=255)

                    ccTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCAc                                                                     <  1:3326294/93‑1 (MQ=255)

                       gAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATTCCCGCTGGCAAACAACTCATCATCACTgg                                                                  >  2:1793808/1‑93 (MQ=255)

                       gAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTgg                                                                  <  2:1546454/93‑1 (MQ=255)

                           ggAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATg                                                              >  2:149810/1‑93 (MQ=255)

                           ggAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATg                                                              >  2:727387/1‑93 (MQ=255)

                           attgTGGTCCTCGCCACTTTCCAGTTCTGCGATGGCTAATGACAGGTAACCACGCACCTGCCCGCTGACAACCACTTCATGGTCATAGGCATg                                                              <  2:1274694/90‑1 (MQ=255)

                            gAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACCACTCATCATCACTGGCAtcg                                                             >  2:566724/1‑91 (MQ=255)

                              gTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCGTCACTGGCATGGTc                                                           >  2:650825/1‑93 (MQ=255)

                              gTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTc                                                           >  1:1119862/1‑93 (MQ=255)

                              gTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTc                                                           >  2:4001307/1‑93 (MQ=255)

                                ggTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCaa                                                         >  2:2331986/1‑93 (MQ=255)

                                  tcatcaCCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAAcc                                                       >  1:4229602/1‑93 (MQ=255)

                                  tcatcaCCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAAcc                                                       >  1:4066469/1‑93 (MQ=255)

                                   catcaCCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCa                                                      >  2:2239109/1‑93 (MQ=255)

                                   catcaCCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCa                                                      >  2:1614097/1‑93 (MQ=255)

                                   catcaCCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCa                                                      >  2:2989893/1‑93 (MQ=255)

                                    atcaCCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCAt                                                     >  2:2809035/1‑93 (MQ=255)

                                    atcaCCACTTTCCAGTTCTGCGATGACTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCAt                                                     <  2:2264719/93‑1 (MQ=255)

                                     tcaCCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATa                                                    >  2:3984681/1‑93 (MQ=255)

                                     tcaCCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATa                                                    >  2:1300798/1‑93 (MQ=255)

                                     tcaCCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATa                                                    <  1:1866487/93‑1 (MQ=255)

                                     tcaCCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATa                                                    <  2:726027/93‑1 (MQ=255)

                                     tcaCCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATa                                                    <  1:1932420/93‑1 (MQ=255)

                                      caccacTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATAc                                                   <  1:4323351/93‑1 (MQ=255)

                                      caccacTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATAc                                                   <  1:2489437/93‑1 (MQ=255)

                                       accacTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATAcc                                                  >  1:4622308/1‑93 (MQ=255)

                                        ccacTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCg                                                 >  1:4552783/1‑93 (MQ=255)

                                             ttCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGAt                                            <  1:1202911/93‑1 (MQ=255)

                                              tCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGAtt                                           >  1:2192231/1‑93 (MQ=255)

                                               ccAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTa                                          >  1:3745612/1‑93 (MQ=255)

                                                cAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTaa                                         <  1:1095193/93‑1 (MQ=255)

                                                       gCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGa                                  <  1:388618/93‑1 (MQ=255)

                                                        cGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGaa                                 <  1:1641408/93‑1 (MQ=255)

                                                        cGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGaa                                 <  1:3033664/93‑1 (MQ=255)

                                                         gATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAAt                                >  1:3811722/1‑93 (MQ=255)

                                                          aTGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATg                               <  1:114669/93‑1 (MQ=255)

                                                           tGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGc                              >  2:2571916/1‑93 (MQ=255)

                                                            ggCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCg                             >  1:1793241/1‑93 (MQ=255)

                                                             gCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGt                            >  1:2766895/1‑93 (MQ=255)

                                                              cTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTg                           >  1:1306067/1‑93 (MQ=255)

                                                               tAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGTCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGa                          >  2:965929/1‑93 (MQ=255)

                                                               tAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGa                          >  2:1976376/1‑93 (MQ=255)

                                                                 aTGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGAtt                        >  2:1036942/1‑93 (MQ=255)

                                                                  tGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTc                       >  2:2957248/1‑93 (MQ=255)

                                                                    tCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCaa                     >  2:3861993/1‑93 (MQ=255)

                                                                     cAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCACCCATACCGTCGATTAATGCCAGAATGCGTGATTCaaa                    <  1:2221914/93‑1 (MQ=255)

                                                                     cAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCaaa                    <  2:3201255/93‑1 (MQ=255)

                                                                      aGGTGGCCACGCACATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAAc                   <  1:1936648/93‑1 (MQ=255)

                                                                      aGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAAc                   >  1:3144958/1‑93 (MQ=255)

                                                                      aGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAAc                   <  1:533381/93‑1 (MQ=255)

                                                                      aGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAAc                   >  1:708484/1‑93 (MQ=255)

                                                                      aGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAAc                   >  1:3144957/1‑93 (MQ=255)

                                                                       ggTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACt                  <  2:2203748/93‑1 (MQ=255)

                                                                       ggTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACt                  <  2:3710648/93‑1 (MQ=255)

                                                                       ggTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACt                  <  2:2203749/93‑1 (MQ=255)

                                                                           gCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACAATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCcgc              <  2:3447210/93‑1 (MQ=255)

                                                                           gCCACGCAAATACCCACTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCcgc              <  1:4018094/93‑1 (MQ=255)

                                                                            ccACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCcgcg             <  2:241924/93‑1 (MQ=255)

                                                                             cACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGa            >  1:3160771/1‑93 (MQ=255)

                                                                              aCGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGat           <  1:796986/93‑1 (MQ=255)

                                                                              aCGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGat           >  1:56510/1‑93 (MQ=255)

                                                                              aCGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGat           >  1:3075886/1‑93 (MQ=255)

                                                                               cGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGatc          <  2:904436/93‑1 (MQ=255)

                                                                               cGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGatc          <  1:3634061/93‑1 (MQ=255)

                                                                               cGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGatc          <  1:3463317/93‑1 (MQ=255)

                                                                                gCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGatca         >  2:61709/1‑93 (MQ=255)

                                                                                gCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGatca         <  2:528667/93‑1 (MQ=255)

                                                                                gCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGatca         <  1:4061252/93‑1 (MQ=255)

                                                                                gCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGatca         <  2:1077708/93‑1 (MQ=255)

                                                                                ccAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGatca         <  1:4061251/92‑1 (MQ=255)

                                                                                   aaTACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCAtct      >  2:3174206/1‑93 (MQ=255)

                                                                                   aaTACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCAtct      >  1:810582/1‑93 (MQ=255)

                                                                                   aaTACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCAtct      >  1:1729922/1‑93 (MQ=255)

                                                                                    aTACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCAtctt     <  2:2623860/93‑1 (MQ=255)

                                                                                      aCCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCAtcttct   >  2:2817546/1‑93 (MQ=255)

                                                                                      aCCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCAtcttct   >  2:2445988/1‑93 (MQ=255)

                                                                                       cccGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCATcttctt  <  2:3823478/93‑1 (MQ=255)

GCTAACGGTCGTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCACTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCATCTTCTT  >  NC_000913/2442016‑2442195

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