BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,853,930	G→A	D524N (GAT→AAT)	fhlA →	FhlA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,853,930	0	G	A	87.6%	757.5 / 60.0	259	D524N (GAT→AAT)	fhlA	FhlA
Reads supporting (aligned to +/- strand): ref base G (16/16); new base A (124/103); total (140/119)
Fisher's exact test for biased strand distribution p-value = 7.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.20e-01

Only 100 of 262 total aligned reads displayed.

GCTGCGTGTATTGCAGGAACAGGAGTTTGAACGTCTCGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGA  >  NC_000913/2853839‑2854021

gctgcGTGTATTGCAGGAACAGGAGTTTGAACGTCTCGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCGa                                                                                            <  1:6494977/93‑1 (MQ=255)

 ctgcGTGTATTGCAGGAACAGGAGTTTGAACGTCTCGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAAt                                                                                           <  2:4993607/93‑1 (MQ=255)

  tgcGTGTATTGCAGGAACAGGAGTTTGAACGTCTCGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATc                                                                                          <  1:5891047/93‑1 (MQ=255)

  tgcGTGTATTGCAGGAACAGGAGTTTGAACGTCTCGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATc                                                                                          >  2:5592123/1‑93 (MQ=255)

    cGTGTATTCCAGGAACAGGAGTTTGAACGTCTCGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTg                                                                                        >  1:320124/1‑93 (MQ=255)

     gtgtATTGCAGGAACAGGAGTTTGAACGTCTCGGCAGCAACAAAATCATTCAGACGGACGTGCGTTTAATCGCCGCGACTAACCGCAATCTGa                                                                                       >  2:2182677/1‑93 (MQ=255)

     gtgtATTGCAGGAACAGGAGTTTGAACGTCTCGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGa                                                                                       >  1:875279/1‑93 (MQ=255)

     gtgtATTGCAGGAACAGGAGTTTGAACGTCTCGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGa                                                                                       <  1:2635577/93‑1 (MQ=255)

      tgtATTGCAGGAACAGGAGTTTGAACGTCTCGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGaa                                                                                      >  1:4384245/1‑93 (MQ=255)

      tgtATTGCAGGAACAGGAGTTTGAACGTCTCGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGaa                                                                                      <  1:56153/93‑1 (MQ=255)

      tgtATTGCAGGAACAGGAGTTTGAACGTCTCGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGaa                                                                                      >  2:6751671/1‑93 (MQ=255)

      tgtATTGCAGGAACAGGAGTTTGAACGTCTCGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGaa                                                                                      <  1:2278297/93‑1 (MQ=255)

        tATTGCAGGAACAGGAGTTTGAACGTCTCGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCGATCTGaaaa                                                                                    >  2:2034308/1‑93 (MQ=255)

        tATTGCAGGAACAGGAGTTTGAACGTCTCGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGaaaa                                                                                    >  2:3149251/1‑93 (MQ=255)

        tATTGCAGGAACAGGAGTTTGAACGTCTCGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGaaaa                                                                                    >  2:1842362/1‑93 (MQ=255)

           tGCAGGAACAGGAGTTTGAACGTCTCGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGaaaaaaa                                                                                 <  1:3770423/93‑1 (MQ=255)

           tGCAGGAACAGGAGTTTGAACGTCTCGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGaaaaaaa                                                                                 <  2:4683336/93‑1 (MQ=255)

            gCAGGAACAGGAGTTTGAACGTCTCGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAAt                                                                                >  2:6355950/1‑93 (MQ=255)

            gCAGGAACAGGAGTTTGAACGTCTCGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAAt                                                                                <  2:719575/93‑1 (MQ=255)

            gCAGGAACAGGAGTTTGAACGTCTCGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAAt                                                                                <  2:3993335/93‑1 (MQ=255)

               ggAACAGGAGTTTGAACGTCTCGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGt                                                                             <  1:2921492/93‑1 (MQ=255)

                gAACAGGAGTTTGAACGTCTCGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTc                                                                            <  2:6849240/93‑1 (MQ=255)

                 aaCAGGAGTTTGAACGTCTCGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCg                                                                           >  1:6912191/1‑93 (MQ=255)

                  aCAGGAGTTTGAACGTCTCGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGc                                                                          <  1:63244/93‑1 (MQ=255)

                  aCAGGAGTTTGAACGTCTCGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGc                                                                          >  1:6707262/1‑93 (MQ=255)

                  aCAGGAGTTTGAACGTCTCGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGc                                                                          <  2:7308547/93‑1 (MQ=255)

                   cAGGAGTTTGAACGTCTCGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGcc                                                                         >  1:452150/1‑93 (MQ=255)

                   cAGGAGTTTGAACGTCTCGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGcc                                                                         <  2:1169011/93‑1 (MQ=255)

                     ggAGTTTGAACGTCTCGGCAGCAACAAAATCATTCAGACGGACGTTCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGa                                                                       >  1:3393103/1‑93 (MQ=255)

                     ggAGTTTGAACGTCTCGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGa                                                                       >  1:3669313/1‑93 (MQ=255)

                     ggAGTTTGAACGTCTCGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGa                                                                       >  1:186786/1‑93 (MQ=255)

                     ggAGTTTGAACGTCTCGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGa                                                                       >  1:5335386/1‑93 (MQ=255)

                     ggAGTTTGAACGTCTCGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGa                                                                       >  1:5020511/1‑93 (MQ=255)

                      gAGTTTGAACGTCTCGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGAc                                                                      <  2:2843205/93‑1 (MQ=255)

                         tttGAACGTCTCGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGt                                                                   <  2:5159840/93‑1 (MQ=255)

                           tGAACGTCTCGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGa                                                                 <  2:4362714/93‑1 (MQ=255)

                           tGAACGTCTCGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGa                                                                 >  2:3387603/1‑93 (MQ=255)

                           tGAACGTCTCGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGa                                                                 <  1:2010456/93‑1 (MQ=255)

                           tGAACGTCTCGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGa                                                                 >  2:7015390/1‑93 (MQ=255)

                              aCGTCTCGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAggt                                                              >  2:579344/1‑91 (MQ=255)

                              aCGTCTCGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGtt                                                              >  2:2756354/1‑93 (MQ=255)

                                gTCTCGGCAGCAACAAAATCATTCAGACGGCCGTGCGTCTAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTcc                                                            >  1:6456385/1‑93 (MQ=255)

                                gTCTCGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTcc                                                            >  1:3023749/1‑93 (MQ=255)

                                gTCTCGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTcc                                                            >  2:1550755/1‑93 (MQ=255)

                                   tcGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTa                                                         <  1:5853830/93‑1 (MQ=255)

                                   tcGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTa                                                         <  1:1565811/93‑1 (MQ=255)

                                      gcagcaACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCg                                                      <  2:1880616/93‑1 (MQ=255)

                                       cagcaACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGa                                                     <  1:1187480/93‑1 (MQ=255)

                                        agcaACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGAt                                                    >  2:4477974/1‑93 (MQ=255)

                                         gcaACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGAtc                                                   <  1:7259389/93‑1 (MQ=255)

                                          caacaaAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGAtct                                                  >  1:733549/1‑93 (MQ=255)

                                           aacaaAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGAtctc                                                 >  2:5492273/1‑93 (MQ=255)

                                           aacaaAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGAtctc                                                 >  1:6812263/1‑93 (MQ=255)

                                            acaaAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATctct                                                >  1:5562112/1‑93 (MQ=255)

                                             caaAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTa                                               <  2:6463597/93‑1 (MQ=255)

                                                aaTCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTa                                            <  1:3363578/93‑1 (MQ=255)

                                                 aTCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTAc                                           <  2:4764439/93‑1 (MQ=255)

                                                   cATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCg                                         <  1:4047214/93‑1 (MQ=255)

                                                    aTTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGc                                        >  1:4270928/1‑93 (MQ=255)

                                                     ttCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGcc                                       >  1:4978582/1‑93 (MQ=255)

                                                     ttCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGcc                                       >  1:4572337/1‑93 (MQ=255)

                                                           cggacgTGCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACg                                 >  1:4745081/1‑93 (MQ=255)

                                                             gacgTGCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTa                               <  1:2047264/93‑1 (MQ=255)

                                                               cgTGCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTAtt                             <  2:5266586/93‑1 (MQ=255)

                                                                gTGCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTc                            <  2:243469/93‑1 (MQ=255)

                                                                 tGCGTCTAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTcc                           <  1:211678/93‑1 (MQ=255)

                                                                 tGCGTCTAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTcc                           <  1:1384626/93‑1 (MQ=255)

                                                                 tGCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTcc                           <  1:2709796/93‑1 (MQ=255)

                                                                 tGCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTcc                           <  1:6099312/93‑1 (MQ=255)

                                                                 tGCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTcc                           >  2:429635/1‑93 (MQ=255)

                                                                  gCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTccc                          >  2:788184/1‑93 (MQ=255)

                                                                  gCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTccc                          <  2:3323769/93‑1 (MQ=255)

                                                                  gCGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTccc                          <  2:1705142/93‑1 (MQ=255)

                                                                   cGTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTCCCg                         >  2:7359484/1‑93 (MQ=255)

                                                                    gTCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTCCCGa                        <  2:4168858/93‑1 (MQ=255)

                                                                     tCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTCCCGAt                       >  1:4163415/1‑93 (MQ=255)

                                                                     tCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTCCCGAt                       >  1:2587435/1‑93 (MQ=255)

                                                                     tCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTCCCGAt                       >  2:7170109/1‑93 (MQ=255)

                                                                     tCTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTCCCGAt                       >  1:796358/1‑93 (MQ=255)

                                                                      cTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTCCCGAtt                      >  1:1674174/1‑93 (MQ=255)

                                                                      cTAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTCCCGAtt                      <  1:3593177/93‑1 (MQ=255)

                                                                       tAATCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTCCCGATTc                     <  2:2688856/93‑1 (MQ=255)

                                                                        aaTCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTCCCGATTCa                    <  1:5119793/93‑1 (MQ=255)

                                                                        aaTCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTCCCGATTCa                    >  1:3879436/1‑93 (MQ=255)

                                                                        aaTCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTCCCGATTCa                    >  2:1290041/1‑93 (MQ=255)

                                                                         aTCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTCCCGATTCAc                   <  2:5084833/93‑1 (MQ=255)

                                                                         aTCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTCCCGATTCAc                   <  1:1016584/93‑1 (MQ=255)

                                                                          tCGCCGCGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTCCCGATTCAcc                  >  1:2387163/1‑93 (MQ=255)

                                                                           cgccgcGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTCCCGATTCACCt                 <  1:7359552/93‑1 (MQ=255)

                                                                             ccgcGACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTCCCGATTCACCTgc               >  1:1940629/1‑93 (MQ=255)

                                                                               gcgACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTCCCGATTCACCTgccg             >  2:589647/1‑93 (MQ=255)

                                                                               gcgACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTCCCGATTCACCTgccg             >  2:7261651/1‑93 (MQ=255)

                                                                                 gACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTCCCGATTCACCTgccgcc           >  2:6014842/1‑93 (MQ=255)

                                                                                 gACTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTCCCGATTCACCTgccgcc           <  2:1752245/93‑1 (MQ=255)

                                                                                  aCTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTCCCGATTCACCTGCCGCCa          >  1:2185669/1‑93 (MQ=255)

                                                                                   cTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTCCCGATTCCCCTGCCGCCAc         <  1:2792001/93‑1 (MQ=255)

                                                                                   cTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTCCCGATTCACCTGCCGCCAc         <  2:1043471/93‑1 (MQ=255)

                                                                                   cTAACCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTCCCGATTCACCTGCCGCCAc         <  1:624860/93‑1 (MQ=255)

                                                                                      aCCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTAc      >  2:2997120/1‑93 (MQ=255)

                                                                                      aCCGCAATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTAc      >  2:3481636/1‑93 (MQ=255)

GCTGCGTGTATTGCAGGAACAGGAGTTTGAACGTCTCGGCAGCAACAAAATCATTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGA  >  NC_000913/2853839‑2854021

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