BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,185,017	Δ1 bp	coding (38/339 nt)	b2107 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,185,014	0	T	.	96.5%	614.8 / 14.0	143	coding (35/339 nt)	b2107	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (2/3); new base . (70/67); total (72/71)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 159 total aligned reads displayed.

TATTGTATGGCAGGAAAGATATGCGACTGATATTACAGATCCCCAAAGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGCTTTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGTAGT  >  NC_000913/2184922‑2185102

tatTGTATGGCAGGAAAGATATGCGACTGATATTACAGATCCCCAAAGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGCt                                                                                          <  2:4117359/93‑1 (MQ=255)

tatTGTATGGCAGGAAAGATATGCGACTGATATTACAGATCCCCAAAGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGCt                                                                                          >  2:3938206/1‑93 (MQ=255)

  ttGTATGGCAGGAAAGATATGCGACTGATATTACAGATCCCCAAAGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGCttt                                                                                        <  2:3100490/93‑1 (MQ=255)

    gTATGGCAGGAAAGATATGCGACTGATATTACAGATCCCCAAAGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGCtttgc                                                                                      <  1:3806110/93‑3 (MQ=255)

     tATGGCAGGAAAGATATGCGACTGATATTACAGATCCCCAAAGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGCtttgct                                                                                     <  2:4234378/93‑4 (MQ=255)

     tATGGCAGGAAAGATATGCGACTGATATTACAGATCCCCAAAGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGCtttgct                                                                                     >  1:2891126/1‑90 (MQ=255)

      aTGGCAGGAAAGATATGCGACTGATATTACAGATCCCCAAAGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGCtttgctg                                                                                    <  2:4008681/93‑5 (MQ=255)

       tGGCAGGAAAGATATGCGACTGATATTACAGATCCCCAAAGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTgg                                                                                  >  1:1241802/1‑93 (MQ=255)

         gCAGGAAAGATATGCGACTGATATTACAGATCCCCAAAGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGtt                                                                                <  1:2089706/93‑1 (MQ=255)

          cAGGAAAGATATGCGACTGATATTACAGATCCCCAAAGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTa                                                                               <  2:4709976/93‑1 (MQ=255)

           aGGAAAGATATGCGACTGATATTACAGATCCCCAAAGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTAc                                                                              >  2:614401/1‑93 (MQ=255)

           aGGAAAGATATGCGACTGATATTACAGATCCCCAAAGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTAc                                                                              <  1:2360496/93‑1 (MQ=255)

           aGGAAAGATATGCGACTGATATTACAGATCCCCAAAGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTAc                                                                              <  1:282044/93‑1 (MQ=255)

            ggAAAGATATGCGACTGATATTACAGATCCCCAAAGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGCTTTTGCTGGTTAc                                                                              >  2:1391829/1‑93 (MQ=255)

            ggAAAGATATGCGACTGATATTACAGATCCCCAAAGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTAcc                                                                             >  1:3698607/1‑93 (MQ=255)

            ggAAAGATATGCGACTGATATCACAGATCCCCAAAGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTAcc                                                                             >  2:1264631/1‑93 (MQ=255)

              aaaGATATGCGACTGATATTACAGATCCCCAAAGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAg                                                                           <  1:3222214/93‑1 (MQ=255)

              aaaGATATGCGACTGATATTACAGATCCCCAAAGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAg                                                                           <  2:2365951/93‑1 (MQ=255)

                 gATATGCGACTGATATTACAGATCCCCAAAGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTgc                                                                        <  2:946597/93‑1 (MQ=255)

                   tatGCGACTGATATTACAGATCCCCAAAGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTgccg                                                                      >  1:480176/1‑93 (MQ=255)

                   tatGCGACTGATATTACAGATCCCCAAAGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTgccg                                                                      >  1:4030583/1‑93 (MQ=255)

                     tGCGACTGATATTACAGATCCCCAAAGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTgccgcc                                                                    >  1:2429588/1‑93 (MQ=255)

                      gCGACTGATATTACAGATCCCCAAAGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCt                                                                   <  1:4578712/93‑1 (MQ=255)

                        gACTGATATTACAGATCCCCAAAGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTgg                                                                 >  1:1869589/1‑93 (MQ=255)

                        gACTGATATTACAGATCCCCAAAGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTgg                                                                 >  1:4639364/1‑93 (MQ=255)

                         aCTGATATTACAGATCCCCAAAGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGCGCCGCCTggg                                                                >  2:3289637/1‑93 (MQ=255)

                           tGATATTACAGATCCCCAAAGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGcc                                                              >  2:1477058/1‑93 (MQ=255)

                                  aCAGATCCCCAAAGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAg                                                       <  1:3859629/93‑1 (MQ=255)

                                  aCAGATCCCCAAAGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAg                                                       <  2:3511432/93‑1 (MQ=255)

                                  aCAGATCCCCAAAGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAg                                                       >  2:4356291/1‑93 (MQ=255)

                                    aGATCCCCAAAGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGcc                                                     >  2:1336970/1‑93 (MQ=255)

                                     gATCCCCAAAGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGCTTTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGcc                                                     >  2:3320733/1‑93 (MQ=255)

                                     gATCCCCAAAGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCa                                                    <  2:3073659/93‑1 (MQ=255)

                                      aTCCCCAAAGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGATGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCAc                                                   <  1:885039/93‑1 (MQ=255)

                                       tCCCCAAAGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCAcc                                                  <  2:1970439/93‑1 (MQ=255)

                                       tCCCCAAAGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCAcc                                                  >  1:350094/1‑93 (MQ=255)

                                        ccccAAAGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACcg                                                 >  2:2034008/1‑93 (MQ=255)

                                         cccAAAGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACcgc                                                <  2:3454371/93‑1 (MQ=255)

                                         cccAAAGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACcgc                                                <  2:2570331/93‑1 (MQ=255)

                                          ccAAAGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACcgcg                                               >  1:2582650/1‑93 (MQ=255)

                                           cAAAGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCgg                                              <  2:2972884/93‑1 (MQ=255)

                                           cAAAGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCgg                                              >  2:303619/1‑93 (MQ=255)

                                            aaaGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCggg                                             >  2:2871038/1‑93 (MQ=255)

                                             aaGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGt                                            >  1:531738/1‑93 (MQ=255)

                                             aaGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGt                                            <  1:605098/93‑1 (MQ=255)

                                              aGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGtt                                           <  1:2921049/93‑1 (MQ=255)

                                              aGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGtt                                           <  2:1997808/93‑1 (MQ=255)

                                                tGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCGGCCTGGGCCGCACCAGCCACCGCGGGTGCg                                         >  2:2587872/1‑93 (MQ=255)

                                                tGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCg                                         <  2:4033961/93‑1 (MQ=255)

                                                 ggAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGa                                        >  2:686926/1‑93 (MQ=255)

                                                 ggAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGa                                        >  1:3377286/1‑93 (MQ=255)

                                                   agagTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGAcc                                      >  1:1821439/1‑93 (MQ=255)

                                                      gTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAAt                                   <  1:3283318/93‑1 (MQ=255)

                                                      gTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAAt                                   >  1:4506083/1‑93 (MQ=255)

                                                       tttATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATa                                  <  2:4746365/93‑1 (MQ=255)

                                                        ttATGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATAc                                 <  2:1565049/93‑1 (MQ=255)

                                                          aTGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCt                               <  1:977928/93‑1 (MQ=255)

                                                          aTGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCt                               >  2:633173/1‑93 (MQ=255)

                                                          aTGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCt                               <  1:4372423/93‑1 (MQ=255)

                                                           tGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTc                              >  2:4610193/1‑93 (MQ=255)

                                                           tGACCATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTc                              >  2:4610194/1‑93 (MQ=255)

                                                               cATTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGa                          >  2:188545/1‑93 (MQ=255)

                                                                aTTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGaa                         >  2:4080772/1‑93 (MQ=255)

                                                                aTTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGaa                         >  2:2904314/1‑93 (MQ=255)

                                                                aTTAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGaa                         >  2:2164044/1‑93 (MQ=255)

                                                                 ttAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAAt                        <  2:2651190/93‑1 (MQ=255)

                                                                  tAAAAATAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAAtt                       <  2:3121753/93‑1 (MQ=255)

                                                                   aaaaaTAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAAttt                      >  2:1364753/1‑93 (MQ=255)

                                                                     aaaTAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCt                    >  2:1095394/1‑93 (MQ=255)

                                                                     aaaTAAGAGGTTGCTGGGTGCG‑ATTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCt                    <  1:3641027/93‑1 (MQ=255)

                                                                      aaTAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTa                   >  1:3908059/1‑93 (MQ=255)

                                                                      aaTAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTa                   >  1:3992151/1‑93 (MQ=255)

                                                                      aaTAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGACGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTa                   >  1:4279011/1‑93 (MQ=255)

                                                                       aTAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTaa                  >  2:2181710/1‑93 (MQ=255)

                                                                        tAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAg                 >  1:2730398/1‑93 (MQ=255)

                                                                        tAAGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAg                 >  1:789226/1‑93 (MQ=255)

                                                                          aGATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTa               <  1:3724251/93‑1 (MQ=255)

                                                                          aGATGTTGCTGGCTGCGC‑TTTGCTGGTTGCCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTa               <  1:209344/93‑1 (MQ=255)

                                                                           gATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTAt              <  2:1239431/93‑1 (MQ=255)

                                                                           gATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTAt              <  1:1765260/93‑1 (MQ=255)

                                                                           gATGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTAt              <  1:1982870/93‑1 (MQ=255)

                                                                            aTGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATg             <  1:3834510/93‑1 (MQ=255)

                                                                             tGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGa            >  2:1714308/1‑93 (MQ=255)

                                                                             tGTTGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACACATACCTCGGGAATTTCTAAGTATGa            >  2:1307848/1‑93 (MQ=255)

                                                                              gTTGCTGGGTGCGCTTTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGa            <  1:1765296/93‑1 (MQ=255)

                                                                               ttGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGt          <  2:2494624/93‑1 (MQ=255)

                                                                                tGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGtt         <  1:874853/93‑1 (MQ=255)

                                                                                tGCTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGtt         <  2:728647/93‑1 (MQ=255)

                                                                                  cTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTaa       <  2:277452/93‑1 (MQ=255)

                                                                                  cTGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTaa       <  1:4096624/93‑1 (MQ=255)

                                                                                  cTGGGTGCGC‑TGTGCTGGTTACCAGTGCCGCCTGGGCCGCGCCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTaa       <  2:233024/93‑1 (MQ=255)

                                                                                   tGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAag      >  2:2136043/1‑93 (MQ=255)

                                                                                   tGGGTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAag      <  2:1073101/93‑1 (MQ=255)

                                                                                    gggTGCGC‑TTTGCTGTTTACCAGTGCCGCCTGGTCCGAGCCAGCATCCGCGGGTTCGTCCACTACCTCGGGAATTTCTAAGTATGAGTTAagc     <  2:4715516/93‑2 (MQ=255)

                                                                                    gggTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAagc     <  2:4230574/93‑2 (MQ=255)

                                                                                    gggTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAagc     >  1:465260/1‑92 (MQ=255)

                                                                                     ggTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAagca    <  2:3190003/93‑3 (MQ=255)

                                                                                     ggTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAagca    <  1:1367631/93‑3 (MQ=255)

                                                                                      gTGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAagcag   <  1:80280/93‑1 (MQ=255)

                                                                                       tGCGC‑TTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAagcagt  >  2:129933/1‑93 (MQ=255)

TATTGTATGGCAGGAAAGATATGCGACTGATATTACAGATCCCCAAAGTGGAGAGTTTATGACCATTAAAAATAAGATGTTGCTGGGTGCGCTTTTGCTGGTTACCAGTGCCGCCTGGGCCGCACCAGCCACCGCGGGTTCGACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGTAGT  >  NC_000913/2184922‑2185102

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