BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,407,232	G→A	A176T (GCG→ACG)	prmA →	PrmA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,407,232	0	G	A	98.1%	536.3 / ‑2.7	155	A176T (GCG→ACG)	prmA	PrmA
Reads supporting (aligned to +/- strand): ref base G (2/1); new base A (57/95); total (59/96)
Fisher's exact test for biased strand distribution p-value = 5.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 171 total aligned reads displayed.

TCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGCGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGC  >  NC_000913/3407142‑3407321

tctctGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCgcgacg                                                                                         >  2:1537482/1‑93 (MQ=255)

 ctctGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGacgc                                                                                        <  1:1664044/93‑1 (MQ=255)

  tctGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCt                                                                                       <  2:2380153/93‑1 (MQ=255)

    tgtgCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGa                                                                                     >  1:590446/1‑93 (MQ=255)

      tgCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGaaa                                                                                   <  1:1191030/93‑1 (MQ=255)

        ccTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACt                                                                                 <  1:400146/93‑1 (MQ=255)

         cTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTg                                                                                <  2:1174955/93‑1 (MQ=255)

          tGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTgg                                                                               >  1:1121731/1‑93 (MQ=255)

          tGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTgg                                                                               <  2:3584620/93‑1 (MQ=255)

           gCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTggg                                                                              <  1:2638220/93‑1 (MQ=255)

           gCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTggg                                                                              <  1:4488872/93‑1 (MQ=255)

           gCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTggg                                                                              <  1:552581/93‑1 (MQ=255)

           gCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTggg                                                                              >  2:2756926/1‑93 (MQ=255)

            cAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGt                                                                             <  1:4643505/93‑1 (MQ=255)

             aaTGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTg                                                                            >  2:1011490/1‑93 (MQ=255)

              aTGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTgc                                                                           <  1:1729446/93‑1 (MQ=255)

              aTGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTgc                                                                           <  2:557771/93‑1 (MQ=255)

               tGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTgca                                                                          <  1:3685472/93‑1 (MQ=255)

                 gCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTgcagc                                                                        <  1:3217778/93‑1 (MQ=255)

                   tCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCaa                                                                      >  1:4724417/1‑93 (MQ=255)

                    tcaCAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCCGCGCTGAAACTGGGTGCAGCaaa                                                                     >  2:1981147/3‑93 (MQ=255)

                    cGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCaaa                                                                     >  1:1165647/1‑93 (MQ=255)

                      aCAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAg                                                                   >  1:1346171/1‑93 (MQ=255)

                      aCAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAg                                                                   >  2:1554170/1‑93 (MQ=255)

                        aGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGcc                                                                 >  1:2996683/1‑93 (MQ=255)

                        aGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGcc                                                                 <  1:153428/93‑1 (MQ=255)

                        aGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGcc                                                                 <  1:3472117/93‑1 (MQ=255)

                          ccTCGATTTATCCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCAt                                                               <  1:3033671/93‑1 (MQ=255)

                          ccTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCAt                                                               >  2:2309814/1‑93 (MQ=255)

                           cTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCAtt                                                              >  2:1392723/1‑93 (MQ=255)

                            tCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTg                                                             >  1:4244275/1‑93 (MQ=255)

                            tCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTg                                                             <  1:3889150/93‑1 (MQ=255)

                            tCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTg                                                             >  2:3615878/1‑93 (MQ=255)

                                tttAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTAt                                                         >  2:1982020/1‑93 (MQ=255)

                                tttAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTAt                                                         <  1:3346281/93‑1 (MQ=255)

                                 ttAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTAtt                                                        <  1:15909/93‑1 (MQ=255)

                                 ttAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTAtt                                                        >  2:1782254/1‑93 (MQ=255)

                                  tAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTg                                                       <  1:3640070/93‑1 (MQ=255)

                                   aaCCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGa                                                      >  1:1457/1‑93 (MQ=255)

                                   aaCCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGa                                                      <  2:2137838/93‑1 (MQ=255)

                                   aaCCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGa                                                      <  1:660954/93‑1 (MQ=255)

                                    aCCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGat                                                     <  2:3812890/93‑1 (MQ=255)

                                     ccGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGata                                                    <  1:4518840/93‑1 (MQ=255)

                                        gTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCg                                                 >  1:2078144/1‑93 (MQ=255)

                                         tAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGa                                                <  1:2947245/93‑1 (MQ=255)

                                          aaaaCGGTCCTCGCCTTTGGCTGTGGTTCCGGCATGCTGGCGATCGCGACGCTGCAACTGGGTGCAGCAAAAGCCATTGGCATTGATATCGAt                                               <  1:41124/93‑1 (MQ=255)

                                          aaaaCAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGAt                                               <  1:1266294/93‑1 (MQ=255)

                                           aaaCAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGAta                                              <  1:427699/93‑2 (MQ=255)

                                           aaaCAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATc                                              >  1:4149653/1‑93 (MQ=255)

                                           aaaCAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATc                                              <  1:76467/93‑1 (MQ=255)

                                            aaCAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATcc                                             <  1:1737750/93‑1 (MQ=255)

                                              aaGTCATCGACTTGGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAATCCATTGGTATTGATATCGATCCGc                                           <  1:3460605/92‑1 (MQ=255)

                                               aGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAACGCAATTGGAATTGATATCGATCCGCa                                          <  2:4227276/93‑1 (MQ=255)

                                               aGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCa                                          <  1:1909266/93‑1 (MQ=255)

                                                 tCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAgg                                        <  1:4188406/93‑1 (MQ=255)

                                                  cATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGc                                       <  1:215353/93‑1 (MQ=255)

                                                  cATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGc                                       >  1:2282934/1‑93 (MQ=255)

                                                   aTCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCg                                      >  1:3085356/1‑93 (MQ=255)

                                                     cgtcTTTGGCTGTGGTTCCGGTATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGAt                                    <  2:30001/90‑1 (MQ=255)

                                                      gACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGAtt                                   >  1:3746314/1‑93 (MQ=255)

                                                      caCTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGAtt                                   <  2:474824/92‑1 (MQ=255)

                                                       aCTTTTGCTGTGGTTCCGGCATTCTGGCGATGGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTc                                  <  2:3543137/93‑1 (MQ=255)

                                                       aCTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTc                                  >  2:2534870/1‑93 (MQ=255)

                                                       aCTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTc                                  >  2:2534871/1‑93 (MQ=255)

                                                        cTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAATCCATTGGTATTGATATCGATCCGCAGGCGATTCa                                 <  1:1039739/93‑1 (MQ=255)

                                                        cTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCa                                 <  1:1039737/93‑1 (MQ=255)

                                                        cTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCa                                 <  2:2229481/93‑1 (MQ=255)

                                                         tttGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAg                                <  2:4569344/93‑1 (MQ=255)

                                                         tttGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAg                                <  1:49225/93‑1 (MQ=255)

                                                           tGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGc                              <  2:2041243/93‑1 (MQ=255)

                                                            ggCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGcc                             <  2:1665627/93‑1 (MQ=255)

                                                               tgtgGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGc                          <  2:3062196/93‑1 (MQ=255)

                                                                gtgGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGcc                         >  1:2548319/1‑93 (MQ=255)

                                                                 tgGTTCCGGCATTCTGGCGATCGCGGCGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCcg                        <  2:3995456/93‑1 (MQ=255)

                                                                 tgGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCcg                        >  1:3264561/1‑93 (MQ=255)

                                                                    tttcGGCATTCTGGCGATCGCGATGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGa                     <  2:2588894/90‑1 (MQ=255)

                                                                    ttCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGa                     <  1:2467129/93‑1 (MQ=255)

                                                                    ttCAGGCATTCTGGCGATCGCGACGCTGAAGTTGGGTGCAGCAAAAGCCATTGTTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGa                     <  2:2765271/93‑1 (MQ=255)

                                                                     tCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGAt                    >  1:3823469/1‑93 (MQ=255)

                                                                       cgccATTCTGCCTACCGCGTCGCTGAAACTGGGGGCAGCAAAAGCCATTGGTATTGATATCTATCCGCAGGCGATTCAGGCCAGCCGCGATaa                  <  2:3189123/90‑1 (MQ=255)

                                                                       cGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATaa                  <  1:1518899/93‑1 (MQ=255)

                                                                       cGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATaa                  <  1:2726033/93‑1 (MQ=255)

                                                                        ggCATTCTGGCTATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAAc                 >  2:1456595/1‑93 (MQ=255)

                                                                        ggCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAAc                 <  2:1812505/93‑1 (MQ=255)

                                                                         gCATTCTGGCGATCGCGGCGCTGAAGCTGGGTGCGGCAAAAGCCATTGGTATTGATATCGACCCGCAGGCTGTTCAGGCCAGCCGCGATAACg                <  2:4282956/93‑1 (MQ=255)

                                                                          cATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGc               >  1:2666533/1‑93 (MQ=255)

                                                                           aTTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGcc              >  1:4388054/1‑93 (MQ=255)

                                                                           aTTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGcc              >  2:3664980/1‑93 (MQ=255)

                                                                            ttCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCg             >  2:3494458/1‑93 (MQ=255)

                                                                            ttCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCg             >  1:3414403/1‑93 (MQ=255)

                                                                            ttCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGGATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCg             >  1:408204/1‑93 (MQ=255)

                                                                             tCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGa            <  2:2691243/93‑1 (MQ=255)

                                                                               tGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAAc          <  2:3263000/93‑1 (MQ=255)

                                                                                   gATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTaa      <  1:2695329/93‑1 (MQ=255)

                                                                                   gATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTaa      >  1:672364/1‑93 (MQ=255)

                                                                                    atagcgacgCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAAt     <  2:4378405/90‑1 (MQ=255)

                                                                                    aTCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAAt     <  1:527528/93‑1 (MQ=255)

                                                                                     tCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATg    >  2:4186008/1‑93 (MQ=255)

                                                                                     tCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATg    <  2:1958387/93‑1 (MQ=255)

                                                                                      cgcgACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATgg   <  2:3089884/93‑1 (MQ=255)

TCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGCGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGC  >  NC_000913/3407142‑3407321

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