BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,637,944	T→C	H68H (CAT→CAC)	uspA →	universal stress protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,637,944	0	T	C	99.8%	1336.7 / ‑5.2	421	H68H (CAT→CAC)	uspA	universal stress protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (243/177); total (243/178)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.98e-01

Only 100 of 454 total aligned reads displayed.

CGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCGATGC  >  NC_000913/3637854‑3638033

 gTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCa                                                                                        >  1:1292112/1‑93 (MQ=255)

 gTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCa                                                                                        >  2:4225663/1‑93 (MQ=255)

  tAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCAc                                                                                       <  1:7498376/93‑1 (MQ=255)

  tAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCAc                                                                                       >  1:6411832/1‑93 (MQ=255)

   aGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACt                                                                                      >  2:6929554/1‑93 (MQ=255)

   aGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACt                                                                                      <  2:2038922/93‑1 (MQ=255)

    gATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTg                                                                                     >  2:2673285/1‑93 (MQ=255)

     aTGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGa                                                                                    <  2:6457993/93‑1 (MQ=255)

      tGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGAc                                                                                   >  2:6059044/1‑93 (MQ=255)

      tGTAAACTACTCTGACCCATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGAc                                                                                   >  1:2302765/1‑93 (MQ=255)

       gTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGAcc                                                                                  >  1:5736380/1‑93 (MQ=255)

         aaaCTACCCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGGGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGa                                                                                <  1:7095744/93‑1 (MQ=255)

           actactCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGc                                                                              >  1:6488662/1‑93 (MQ=255)

           actactCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGc                                                                              <  2:1356174/93‑1 (MQ=255)

                 ctGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCa                                                                        >  1:3451597/1‑93 (MQ=255)

                  tGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCAc                                                                       >  1:4661718/1‑93 (MQ=255)

                  tGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCAc                                                                       <  1:4355889/93‑1 (MQ=255)

                  tGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCAc                                                                       <  1:5837092/93‑1 (MQ=255)

                  tGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCAc                                                                       >  1:7214748/1‑93 (MQ=255)

                   gACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACt                                                                      >  1:1703236/1‑93 (MQ=255)

                     ccTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTaa                                                                    <  1:546463/93‑1 (MQ=255)

                     ccTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTaa                                                                    <  1:740371/93‑1 (MQ=255)

                       tataCACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATg                                                                  >  1:7862315/1‑93 (MQ=255)

                       tataCACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATg                                                                  >  1:7556529/1‑93 (MQ=255)

                       tataCACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATg                                                                  >  1:2644063/1‑93 (MQ=255)

                        ataCACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGc                                                                 >  2:3672238/1‑93 (MQ=255)

                         taCACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCa                                                                >  2:1731867/1‑93 (MQ=255)

                          acacCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAg                                                               >  1:1651020/1‑93 (MQ=255)

                           cacCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAgg                                                              >  2:195050/1‑93 (MQ=255)

                           cacCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAgg                                                              >  2:2594066/1‑93 (MQ=255)

                            acCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGc                                                             >  2:7274403/1‑93 (MQ=255)

                            acCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGc                                                             <  1:5254420/93‑1 (MQ=255)

                            acCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGc                                                             <  2:4905536/93‑1 (MQ=255)

                               gggCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTAc                                                          >  1:6259968/1‑93 (MQ=255)

                                ggCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTAATGCAGGCTAcc                                                         <  2:5884704/93‑1 (MQ=255)

                                ggCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTAcc                                                         >  1:6156277/1‑93 (MQ=255)

                                ggCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTAcc                                                         >  2:1568895/1‑93 (MQ=255)

                                ggCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTAcc                                                         >  1:2887924/1‑93 (MQ=255)

                                 gCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTAccc                                                        <  2:5547560/93‑1 (MQ=255)

                                   ttATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCaa                                                      >  2:7622379/1‑93 (MQ=255)

                                   ttATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCaa                                                      >  2:3108448/1‑93 (MQ=255)

                                    tATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAAt                                                     <  1:396567/93‑1 (MQ=255)

                                      ttGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCa                                                   <  1:3571502/93‑1 (MQ=255)

                                      ttGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCa                                                   >  1:814395/1‑93 (MQ=255)

                                      ttGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCa                                                   >  2:5075420/1‑93 (MQ=255)

                                      ttGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCa                                                   <  1:5305684/93‑1 (MQ=255)

                                        gATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACt                                                 >  1:4074288/1‑93 (MQ=255)

                                        gATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGAACTGACTGAGCTTTCCACTAATGCAGGCTACCCAATCACt                                                 >  1:7609342/1‑93 (MQ=255)

                                          tgtgAATCTGGGTGATATGCAGAAACGGATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGa                                               >  2:495417/1‑93 (MQ=255)

                                          tgtgAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGa                                               <  2:3007744/93‑1 (MQ=255)

                                           gtgAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGaa                                              >  2:4396045/1‑93 (MQ=255)

                                            tgAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGaaa                                             <  2:4694428/93‑1 (MQ=255)

                                            tgAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACACACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGaaa                                             <  1:3374023/93‑1 (MQ=255)

                                             gAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAAc                                            >  2:3310135/1‑93 (MQ=255)

                                             gAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAAc                                            >  2:1502769/1‑93 (MQ=255)

                                             gAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAAc                                            <  1:3842217/93‑1 (MQ=255)

                                               aTCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAAccc                                          <  1:5108129/93‑1 (MQ=255)

                                                tCTGGGTGATATGCAGAAACGTATCTTTGAAGAGACACATCACGCACTGACCGTTCTTTCCACTAATGCAGGCTACCCAATTATTTAAACCCt                                         >  2:7041630/1‑93 (MQ=255)

                                                tCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCt                                         >  1:7356930/1‑93 (MQ=255)

                                                tCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCt                                         >  2:6748298/1‑93 (MQ=255)

                                                tCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCt                                         >  1:2527773/1‑93 (MQ=255)

                                                  tGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGa                                       >  1:4452374/1‑93 (MQ=255)

                                                    ggTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGagc                                     <  2:7823917/93‑1 (MQ=255)

                                                     gTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGagcg                                    >  2:5218438/1‑93 (MQ=255)

                                                      tGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGagcgg                                   >  1:1933455/1‑93 (MQ=255)

                                                        atatGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGagcggca                                 <  1:7730459/93‑1 (MQ=255)

                                                         tatGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGagcggcag                                <  2:6212249/93‑1 (MQ=255)

                                                         tatGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGagcggcag                                <  1:938772/93‑1 (MQ=255)

                                                               gAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGa                          <  2:81278/93‑1 (MQ=255)

                                                               gAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGa                          >  1:6384376/1‑93 (MQ=255)

                                                               gAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGa                          >  1:6887198/1‑93 (MQ=255)

                                                                aaaCGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGAc                         <  2:2892869/93‑1 (MQ=255)

                                                                aaaCGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGAc                         <  1:2118016/93‑1 (MQ=255)

                                                                 aaCGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGAcc                        >  2:6343853/1‑93 (MQ=255)

                                                                 aaCGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGAcc                        <  2:5710500/93‑1 (MQ=255)

                                                                 aaCGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGAcc                        >  2:123121/1‑93 (MQ=255)

                                                                 aaCGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGAcc                        >  2:2226404/1‑93 (MQ=255)

                                                                  aCGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCt                       >  1:3203029/1‑93 (MQ=255)

                                                                   cGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTg                      >  1:6602603/1‑93 (MQ=255)

                                                                   cGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTg                      <  2:5365769/93‑1 (MQ=255)

                                                                    gCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTgg                     <  2:3399584/93‑1 (MQ=255)

                                                                    gCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTgg                     >  2:2640024/1‑93 (MQ=255)

                                                                       tctcTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGcg                  >  1:5674027/1‑92 (MQ=255)

                                                                        ctctGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCa                 <  2:3820734/93‑1 (MQ=255)

                                                                         tctGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAg                >  1:4971395/1‑93 (MQ=255)

                                                                          ctGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAgg               >  1:4844849/1‑93 (MQ=255)

                                                                             aaGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTc            >  1:5191833/1‑93 (MQ=255)

                                                                             aaGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTc            >  1:3999343/1‑93 (MQ=255)

                                                                                 gaCACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGt        >  2:7414822/1‑93 (MQ=255)

                                                                                 gaCACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGt        <  1:121349/93‑1 (MQ=255)

                                                                                 gaCACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGt        >  1:1035899/1‑93 (MQ=255)

                                                                                    acaTCACTCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCGa     <  1:4408270/93‑1 (MQ=255)

                                                                                    acaTCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCGa     >  2:4094658/1‑93 (MQ=255)

                                                                                     catcatGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAGCCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCGAt    <  2:3574794/93‑1 (MQ=255)

                                                                                     catcacGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCGAt    <  1:4325988/93‑1 (MQ=255)

                                                                                     catcacGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCGAt    <  1:3000576/93‑1 (MQ=255)

                                                                                     catcacGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCGAt    <  2:2388938/93‑1 (MQ=255)

                                                                                     catcacGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCGAt    >  1:1513845/1‑93 (MQ=255)

                                                                                       tcacGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCGATGc  <  1:6696443/93‑1 (MQ=255)

                                                                                       tcacGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCGATGc  >  2:283146/1‑93 (MQ=255)

CGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCGATGC  >  NC_000913/3637854‑3638033

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