BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,084,090	A→G	Q165Q (CAA→CAG)	fdhD →	affects formate dehydrogenase‑N

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,084,090	0	A	G	100.0%	634.3 / NA	198	Q165Q (CAA→CAG)	fdhD	affects formate dehydrogenase‑N
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (110/88); total (110/88)

Only 100 of 228 total aligned reads displayed.

GGTGCAGCCGCTACCGTTCACCCAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCATGTGGCGCT  >  NC_000913/4084000‑4084179

 gTGCAGCCGCTACCGTTCACCCAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAgctg                                                                                        >  1:3318567/1‑93 (MQ=255)

 gTGCAGCCGCTACCGTTCACCCAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAgctg                                                                                        >  1:1620013/1‑93 (MQ=255)

 gTGCAGCCGCTACCGTTCACCCAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAgctg                                                                                        <  1:2490809/93‑1 (MQ=255)

  tGCAGCCGCTACCGTTCACCCAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAgctga                                                                                       <  2:2627378/93‑1 (MQ=255)

  tGCAGCCGCTACCGTTCACCCAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAgctga                                                                                       <  1:3119370/93‑1 (MQ=255)

    gcgCCGCTACCGTTCACCCAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAgctgact                                                                                     <  1:3669636/91‑1 (MQ=255)

    cAGCCGCTACCGTTCACCCAGACGTTTGATCTCGACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAgctgact                                                                                     >  1:3221890/1‑93 (MQ=255)

    cAGCCGCTACCGTTCACCCAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAgctgact                                                                                     <  2:1107075/93‑1 (MQ=255)

    cAGCCGCTACCGTTCACCCAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAgctgact                                                                                     <  2:4607023/93‑1 (MQ=255)

    cAGCCGCTACCGTTCACCCAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAgctgact                                                                                     <  2:4590226/93‑1 (MQ=255)

    cAGCCGCTACCGTTCACCCAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAgctgact                                                                                     <  1:3044196/93‑1 (MQ=255)

     aGCCGCTACCGTTCACCCAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAgctgactg                                                                                    <  1:4053192/93‑1 (MQ=255)

      gCCGCTAACGTTCACACAGAAGATTGATATCAACAAACTGGATGATGCATTACGTCATCTCGACGATTTTAAGCCAGTGGGGCAGCTGACTgg                                                                                   >  2:2024508/1‑93 (MQ=255)

         gCTACCGTTCACCCAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTg                                                                                <  2:4061172/93‑1 (MQ=255)

           tACCGTTCACCCAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTa                                                                              <  1:1034828/93‑1 (MQ=255)

            aCCGTTCACCCAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTAc                                                                             <  1:1240179/93‑1 (MQ=255)

            aCCGTTCACCCAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTAc                                                                             <  1:2344417/93‑1 (MQ=255)

            aCCGTTCACCCAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTAc                                                                             >  1:4510506/1‑93 (MQ=255)

            aCCGTTCACCCAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTAc                                                                             <  1:2928059/93‑1 (MQ=255)

            aCCGTTCACCCAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTAc                                                                             <  2:1162355/93‑1 (MQ=255)

             ccGTTCACCCAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACt                                                                            <  2:2496762/93‑1 (MQ=255)

             ccGTTCACCCAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACt                                                                            <  2:3006425/93‑1 (MQ=255)

               gTTCACCCAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCa                                                                          >  1:1688315/1‑93 (MQ=255)

                ttCACCCAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCAc                                                                         <  2:2351717/93‑1 (MQ=255)

                ttCACCCAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCAc                                                                         >  2:2164679/1‑93 (MQ=255)

                   aCCCAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCAcgcc                                                                      >  1:2584914/1‑93 (MQ=255)

                    cccAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCAcgccg                                                                     >  2:695641/1‑93 (MQ=255)

                     ccAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGCGCAGCTGACTGGTTGTACTCAcgccgc                                                                    >  1:2141519/1‑93 (MQ=255)

                      cAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCt                                                                   <  2:3947311/93‑1 (MQ=255)

                      cAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCt                                                                   >  2:3830279/1‑93 (MQ=255)

                          cGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCt                                                               <  1:396415/93‑1 (MQ=255)

                           gTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTg                                                              >  2:3237925/1‑93 (MQ=255)

                            tttGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTgg                                                             >  1:4358421/1‑93 (MQ=255)

                             ttGTTCTCAACATACTGGATGATCCAGTGCGTCATCTCAACGATTTTCAGGCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGa                                                            >  2:3022050/1‑93 (MQ=255)

                               taTCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATg                                                          <  1:4614800/92‑1 (MQ=255)

                               gATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGCATg                                                          >  2:1580173/1‑93 (MQ=255)

                                aTCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGt                                                         >  2:161542/1‑93 (MQ=255)

                                   tcAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGc                                                      >  1:4678914/1‑93 (MQ=255)

                                    caacaaACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGcc                                                     >  2:2576248/1‑93 (MQ=255)

                                     aacaaACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCa                                                    <  2:2233739/93‑1 (MQ=255)

                                      acaaACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCAt                                                   <  1:1439489/93‑1 (MQ=255)

                                      acaaACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCAt                                                   >  2:1277820/1‑93 (MQ=255)

                                      acaaACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCAt                                                   >  2:1745938/1‑93 (MQ=255)

                                         aaCTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTg                                                <  1:412789/93‑1 (MQ=255)

                                          aCTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTgg                                               >  2:2337015/1‑93 (MQ=255)

                                           cTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGc                                              >  2:4237364/1‑93 (MQ=255)

                                            tGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTCCTCACGCCGCTGCCTGGATGTTGCCCTCTGGcc                                             >  2:2199329/1‑92 (MQ=255)

                                            tGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCg                                             >  1:1318512/1‑93 (MQ=255)

                                             ggATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGa                                            >  1:3411883/1‑93 (MQ=255)

                                              gatgatGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGaa                                           >  1:4378758/1‑93 (MQ=255)

                                               atgatgCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGaat                                          >  1:2692091/1‑92 (MQ=255)

                                               atgatgCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAAc                                          <  1:1426704/93‑1 (MQ=255)

                                                   tgCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGt                                      >  1:1009421/1‑93 (MQ=255)

                                                   tgCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGt                                      <  2:1944874/93‑1 (MQ=255)

                                                    gCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTc                                     <  2:3504198/93‑1 (MQ=255)

                                                    gCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTc                                     >  2:93000/1‑93 (MQ=255)

                                                    gCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTc                                     <  2:2820706/93‑1 (MQ=255)

                                                    gCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTc                                     <  1:1813687/93‑1 (MQ=255)

                                                    gCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTc                                     <  1:3589765/93‑1 (MQ=255)

                                                     cATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTcg                                    <  1:988280/93‑1 (MQ=255)

                                                      aTTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTcgg                                   >  1:1693309/1‑93 (MQ=255)

                                                      aTTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTcgg                                   >  1:3581392/1‑93 (MQ=255)

                                                       ttACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTcggg                                  >  2:883269/1‑92 (MQ=255)

                                                       ttACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTcggc                                  >  1:4434385/1‑93 (MQ=255)

                                                       ttACGTCATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTcggc                                  >  2:528838/1‑93 (MQ=255)

                                                             cATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCAt                            <  1:1940236/93‑1 (MQ=255)

                                                             cATCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCAt                            <  1:1723084/93‑1 (MQ=255)

                                                              aTCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATg                           <  1:243462/93‑1 (MQ=255)

                                                              aTCTCAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATg                           >  1:1349784/1‑93 (MQ=255)

                                                               tctcAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGa                          <  2:3787007/93‑1 (MQ=255)

                                                               tctcAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGa                          <  2:3146852/93‑1 (MQ=255)

                                                               tctcAACGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGa                          <  1:788266/93‑1 (MQ=255)

                                                                   aaCGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGAc                      >  1:1782246/1‑93 (MQ=255)

                                                                    aCGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACg                     >  1:3436744/1‑93 (MQ=255)

                                                                     cGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGAccg                    >  1:2446466/1‑91 (MQ=255)

                                                                     cGATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGt                    <  2:2773299/93‑1 (MQ=255)

                                                                      gATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTg                   >  1:1157247/1‑93 (MQ=255)

                                                                      gATTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTg                   >  1:3913702/1‑93 (MQ=255)

                                                                       tttttCAGCGAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTgg                  <  1:3883577/92‑1 (MQ=255)

                                                                       aTTTTCAGCCAGTGGGGCGGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTgg                  <  2:2515751/93‑1 (MQ=255)

                                                                       aTTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACgggg                  >  2:3032473/1‑93 (MQ=255)

                                                                       aTTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGtcg                  <  2:2956281/93‑3 (MQ=255)

                                                                       aTTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTgg                  >  2:1053135/1‑93 (MQ=255)

                                                                       aTTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTgg                  >  2:1858403/1‑93 (MQ=255)

                                                                       aTTTTCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTgg                  <  2:3855344/93‑1 (MQ=255)

                                                                        ttttCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTggg                 >  1:1845474/1‑93 (MQ=255)

                                                                        ttttCAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTggg                 >  2:565979/1‑93 (MQ=255)

                                                                          ttCAGCCGGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGGCGGCGGGCATGAAGACGTggggc               >  1:3101651/1‑91 (MQ=255)

                                                                            cAGCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGc             >  2:2108255/1‑93 (MQ=255)

                                                                              gCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCa           >  2:4525434/1‑93 (MQ=255)

                                                                              gCCAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCa           >  2:1462839/1‑93 (MQ=255)

                                                                               ccAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGGCGccct          >  1:3939782/1‑91 (MQ=255)

                                                                                cAGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGAGCATGAAGACGTGGGTCGCCAtg         >  2:858659/1‑93 (MQ=255)

                                                                                 aGTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCAtgt        <  2:3370229/93‑1 (MQ=255)

                                                                                  gTGGGGCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCAtgtg       >  1:2054385/1‑93 (MQ=255)

                                                                                     gggCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAAATGGTCGGCGGGCATGAAGACGTGGGTCGCCATGTGgcg    >  2:698427/1‑93 (MQ=255)

                                                                                      ggCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCATGTGgcgc   <  2:2712843/93‑1 (MQ=255)

                                                                                      ggCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCATGTGgcgc   >  1:4205745/1‑93 (MQ=255)

                                                                                      ggCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCATGTGgcgc   <  1:2642858/93‑1 (MQ=255)

                                                                                       gCAGCTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCATGTGGCGCt  <  2:4291332/93‑1 (MQ=255)

GGTGCAGCCGCTACCGTTCACCCAGACGTTTGATCTCAACAAACTGGATGATGCATTACGTCATCTCAACGATTTTCAGCCAGTGGGGCAACTGACTGGTTGTACTCACGCCGCTGCCTGGATGTTGCCATCTGGCGAACTGGTCGGCGGGCATGAAGACGTGGGTCGCCATGTGGCGCT  >  NC_000913/4084000‑4084179

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