BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,710,261	T→C	E101E (GAA→GAG)	b3548 ←	putative lipase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,710,261	0	T	C	85.2%	754.8 / 120.9	324	E101E (GAA→GAG)	b3548	putative lipase
Reads supporting (aligned to +/- strand): ref base T (25/23); new base C (157/119); total (182/142)
Fisher's exact test for biased strand distribution p-value = 5.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.97e-02

Only 100 of 326 total aligned reads displayed.

GTCAGCGTGGTAGTACTGTCGCCAAACTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGTTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGTCCCAGCGCCCT  >  NC_000913/3710169‑3710353

  cAGCGTGGTAGTACTGTCGCCAAACTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTc                                                                                            <  1:6416369/93‑1 (MQ=255)

   aGCGTGGTAGTACTGTCGCCAAACTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGGCGGTTGTAACCTGCTCa                                                                                           <  2:6861460/93‑1 (MQ=255)

   aGCGTGGTAGTACTGTCGCCAAACTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGTTCa                                                                                           <  2:1844574/93‑1 (MQ=255)

   aGCGTGGTAGTACTGTCGCCAAACTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCa                                                                                           <  2:2474609/93‑1 (MQ=255)

     cGTGGTAGTACTGTCGCCAAACTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCata                                                                                         >  2:2785954/1‑93 (MQ=255)

      gTGGTAGTACTGTCGCCAAACTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCatat                                                                                        >  2:4249049/1‑93 (MQ=255)

      gTGGTAGTACTGTCGCCAAACTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCatat                                                                                        >  1:5719357/1‑93 (MQ=255)

      gTGGTAGTACTGTCGCCAAACTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCatat                                                                                        <  2:6489996/93‑1 (MQ=255)

       tGGTAGTACTGTCGCCAAACTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCAtata                                                                                       >  1:7126096/1‑93 (MQ=255)

        ggTAGTACTGTCGCCAAACTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATaa                                                                                      <  1:4583626/93‑1 (MQ=255)

         gtagtaCTGTCGCCAAACTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATaaa                                                                                     >  1:819735/1‑93 (MQ=255)

         gtagtaCTGTCGCCAAACTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATaaa                                                                                     >  1:3469964/1‑93 (MQ=255)

          tagtaCTGTCGCCAAACTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAg                                                                                    <  1:48727/93‑1 (MQ=255)

           agtaCTGTCGCCAAACTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAgg                                                                                   <  1:4968697/93‑1 (MQ=255)

            gtaCTGTCGCCAAACTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGa                                                                                  >  1:2517061/1‑93 (MQ=255)

             taCTGTCGCCAAACTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGaa                                                                                 >  2:1539677/1‑93 (MQ=255)

             taCTGTCGCCAAACTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGaa                                                                                 <  1:4174859/93‑1 (MQ=255)

              aCTGTCGCCAAACTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGaaa                                                                                <  2:3745405/93‑1 (MQ=255)

              aCTGTCGCCAAACTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGaaa                                                                                <  2:2303929/93‑1 (MQ=255)

               cTGTCGCCAAACTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGaaaa                                                                               <  1:1034583/93‑1 (MQ=255)

                tGTCGCCAAACTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGTTCATATAAAGGAAAAt                                                                              <  2:3842823/93‑1 (MQ=255)

                tGTCGCCAAACTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAAt                                                                              >  1:3340110/1‑93 (MQ=255)

                  tCGCCAAACTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTc                                                                            <  1:1858211/93‑1 (MQ=255)

                   cGCCAAACTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTcc                                                                           <  2:3066363/93‑1 (MQ=255)

                    gCCAAACTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCa                                                                          <  2:4994274/93‑1 (MQ=255)

                    gCCAAACTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCa                                                                          >  1:5037150/1‑93 (MQ=255)

                     ccAAACTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAg                                                                         >  2:4324923/1‑93 (MQ=255)

                     ccAAACTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAg                                                                         >  1:4504725/1‑93 (MQ=255)

                     ccAAACTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAg                                                                         >  2:4595953/1‑93 (MQ=255)

                     ccAAACTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAg                                                                         <  1:6029135/93‑1 (MQ=255)

                      cAAACTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAGc                                                                        >  1:203965/1‑93 (MQ=255)

                       aaaCTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGTTCATATAAAGGAAAACTCCAGcc                                                                       >  2:448542/1‑93 (MQ=255)

                        aaCTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAGccc                                                                      >  1:39484/1‑93 (MQ=255)

                        aaCTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAGccc                                                                      >  2:1197142/1‑93 (MQ=255)

                        aaCTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAGccc                                                                      >  1:1384226/1‑93 (MQ=255)

                         aCTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAGCCca                                                                     >  2:1128996/1‑93 (MQ=255)

                          cTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAGCCcat                                                                    >  1:296967/1‑93 (MQ=255)

                           tCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAGCCcata                                                                   <  1:682494/93‑1 (MQ=255)

                             gTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAGCCcatact                                                                 >  1:920383/1‑93 (MQ=255)

                              ttATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAGCCcatactc                                                                >  1:5297170/1‑93 (MQ=255)

                              ttATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAGCCcatactc                                                                >  2:1658533/1‑93 (MQ=255)

                                ataCATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAGCCcatactcat                                                              >  2:2062543/1‑93 (MQ=255)

                                ataCATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAGCCcatactcat                                                              >  2:2697798/1‑93 (MQ=255)

                                 taCATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAGCCcatactcata                                                             >  2:3705657/1‑93 (MQ=255)

                                   cATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGCCCGGGCCGGTTGTACCCTGCTCATATAAAGGAAAATTCCAGCCcatactcatact                                                           >  1:6923480/1‑93 (MQ=255)

                                   cATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGTTCATATAAAGGAAAATTCCAGCCcatactcatact                                                           >  1:2878305/1‑93 (MQ=255)

                                   cATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGTTCATATAAAGGAAAATTCCAGCCcatactcatact                                                           >  2:568769/1‑93 (MQ=255)

                                    aTCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAGCCCATACTCATActg                                                          >  1:3112965/1‑93 (MQ=255)

                                      cGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAGCCCATACTCATActgct                                                        <  2:4110624/93‑1 (MQ=255)

                                         ggtggtTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGt                                                     <  2:3200104/93‑1 (MQ=255)

                                          gtggtTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTc                                                    >  1:5877054/1‑93 (MQ=255)

                                           tggtTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTca                                                   >  1:2093610/1‑93 (MQ=255)

                                            ggtTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGTTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTcac                                                  <  1:3688904/93‑1 (MQ=255)

                                             gtTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTcaca                                                 >  2:2605479/1‑93 (MQ=255)

                                              ttGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACAt                                                >  1:5225372/1‑93 (MQ=255)

                                               tgtgCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATc                                               <  1:6319887/93‑1 (MQ=255)

                                                 tgCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCca                                             <  2:5944125/93‑1 (MQ=255)

                                                   ccATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGTTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCcaca                                           <  1:7052294/93‑1 (MQ=255)

                                                   ccATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGTTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCcaca                                           >  1:6137077/1‑93 (MQ=255)

                                                   ccATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCcaca                                           >  2:614294/1‑93 (MQ=255)

                                                   ccATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCcaca                                           <  1:1461768/93‑1 (MQ=255)

                                                    cATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCCACAg                                          >  1:273145/1‑93 (MQ=255)

                                                      tCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGTTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCCACAggg                                        <  1:2931968/93‑1 (MQ=255)

                                                      tCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCCACAggg                                        <  2:1061711/93‑1 (MQ=255)

                                                      tCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCCACAggg                                        >  2:76217/1‑93 (MQ=255)

                                                        gtaatgtaatGCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCCACAGGGGt                                      >  2:6878044/1‑93 (MQ=255)

                                                           atgtaatgCCGCGACCGGGCCGGTTGTAACCTGTTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCCACAGGGGTttc                                   >  2:774126/1‑93 (MQ=255)

                                                            tgtaatgCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCCACAGGGGTttct                                  >  2:5202400/1‑93 (MQ=255)

                                                            tgtaatgCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCCACAGGGGTttct                                  <  2:671600/93‑1 (MQ=255)

                                                              taatgCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCCACAGGGGTttcttc                                >  1:3971131/1‑93 (MQ=255)

                                                              taatgCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCCACAGGGGTttcttc                                <  2:7413581/93‑1 (MQ=255)

                                                              taatgCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCCACAGGGGTttcttc                                <  2:3640423/93‑1 (MQ=255)

                                                               aatgCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCCACAGGGGTTTCTTcc                               <  2:62950/93‑1 (MQ=255)

                                                                  gCCGCGACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCCACAGGGGTTTCTTCCGcc                            >  2:2158907/1‑93 (MQ=255)

                                                                     gcgACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCCACAGGGGTTTCTTCCGCCAGa                         >  1:7439090/1‑93 (MQ=255)

                                                                      cgACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCCACAGGGGTTTCTTCCGCCAGAt                        <  2:2936524/93‑1 (MQ=255)

                                                                       gACCGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCCACAGGGGTTTCTTCCGCCAGAtt                       >  2:3405970/1‑93 (MQ=255)

                                                                        aCCGGGCCGGTTGTAACCTGTTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCCACAGGGGTTTCTTCCGCCAGATTg                      >  1:3535756/1‑93 (MQ=255)

                                                                         ccGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACAGCCACAGGGGTTTCTTCCGCCAGATTGa                     <  2:6416985/93‑1 (MQ=255)

                                                                          cGGGCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGAc                    >  2:7049358/1‑93 (MQ=255)

                                                                           gggCCGGTTGTAACCTGCTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACt                   >  1:3215483/1‑93 (MQ=255)

                                                                              ccGGTTGTAACCTGTTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTcc                >  1:567301/1‑93 (MQ=255)

                                                                              ccGGTTGTAACCTGCTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTcc                <  2:4867109/93‑1 (MQ=255)

                                                                               cGGTTGTAACCTGCTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCa               <  1:4095310/93‑1 (MQ=255)

                                                                               cGGTTGTAACCTGCTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCa               <  1:3034620/93‑1 (MQ=255)

                                                                                 gTTGTAACCTGTTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCCACAGGGGTTTCTACCGCCAGATTGACTTCCAGt             >  1:5522887/1‑93 (MQ=255)

                                                                                 gTTGTAACCTGCTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGt             >  2:4929407/1‑93 (MQ=255)

                                                                                 gTTGTAACCTGCTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGt             >  2:3999209/1‑93 (MQ=255)

                                                                                  ttGTAACCTGCTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGTc            >  2:5389738/1‑93 (MQ=255)

                                                                                  ttGTAACCTGCCCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGTc            >  2:4655318/1‑93 (MQ=255)

                                                                                    gTAACCTGCTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGTccc          >  1:113921/1‑93 (MQ=255)

                                                                                    gTAACCTGCTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGTccc          >  2:509816/1‑93 (MQ=255)

                                                                                    gTAACCTGCTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGTccc          >  1:6180266/1‑93 (MQ=255)

                                                                                    gTAACCTGCTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGTccc          >  1:4790964/1‑93 (MQ=255)

                                                                                    gTAACCTGCTCATATAAAAGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGTccc          >  1:4294450/1‑93 (MQ=255)

                                                                                     tAACCTGCTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGTCCCa         <  2:5770238/93‑1 (MQ=255)

                                                                                      aaCCTGCTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGTCCCAg        <  2:5860184/93‑1 (MQ=255)

                                                                                          tGTTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGTCCCAGCGcc    >  1:6638293/1‑93 (MQ=255)

                                                                                           gTTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCCACAGGGGTTTCTTCCGTCAGATTGACTTCCAGTCCCAGCGccc   >  1:1714654/1‑93 (MQ=255)

                                                                                            ttCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGTCCCAGCGCCCt  <  1:97539/93‑1 (MQ=255)

GTCAGCGTGGTAGTACTGTCGCCAAACTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGTTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGTCCCAGCGCCCT  >  NC_000913/3710169‑3710353

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