BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,444,690	T→C	V1001A (GTA→GCA)	ytfN →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,444,690	0	T	C	98.0%	470.4 / ‑0.8	153	V1001A (GTA→GCA)	ytfN	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (2/1); new base C (83/67); total (85/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.13e-01

Only 100 of 160 total aligned reads displayed.

TCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATT  >  NC_000913/4444600‑4444781

tCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGcagcag                                                                                           >  1:1109140/1‑93 (MQ=255)

  gAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGc                                                                                         >  1:455127/1‑93 (MQ=255)

  gAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGc                                                                                         >  2:2004379/1‑93 (MQ=255)

   aGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCg                                                                                        <  2:2666239/93‑1 (MQ=255)

    ggCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGt                                                                                       >  1:897065/1‑93 (MQ=255)

     gCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGtc                                                                                      >  2:2148868/1‑93 (MQ=255)

       tACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTTCCGGAAAGCGCAGCAGGCGtctc                                                                                    <  1:1544994/93‑1 (MQ=255)

       tACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGtctc                                                                                    <  2:3404746/93‑1 (MQ=255)

       tACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCACTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGtctc                                                                                    >  1:31321/1‑93 (MQ=255)

         cacCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCa                                                                                  >  2:3255925/1‑93 (MQ=255)

          acCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAg                                                                                 <  2:107072/93‑1 (MQ=255)

           ccAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGc                                                                                <  2:2398580/93‑1 (MQ=255)

            cAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCg                                                                               >  2:1041131/1‑93 (MQ=255)

            cAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCg                                                                               >  2:1966996/1‑93 (MQ=255)

              aaCCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGAt                                                                             >  1:3664442/1‑93 (MQ=255)

                ccTGTTTACCCTCGATGGTCCGGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGAtgt                                                                           <  1:1282355/93‑1 (MQ=255)

                 cTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGAtgtg                                                                          <  1:517487/93‑1 (MQ=255)

                  tGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATgtgg                                                                         >  2:3066080/1‑93 (MQ=255)

                   gTTTTCCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCTCAGCAGGCGTCTCCAGCGATgtggt                                                                        <  2:607281/93‑1 (MQ=255)

                   gTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGAGTCTCCAGCGATgtggt                                                                        >  1:4710970/1‑93 (MQ=255)

                    tttACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATgtggtg                                                                       >  2:1801061/1‑93 (MQ=255)

                     ttACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGa                                                                      <  1:1489578/93‑1 (MQ=255)

                     ttACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGa                                                                      <  1:3967273/93‑1 (MQ=255)

                     ttACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGa                                                                      >  2:843839/1‑93 (MQ=255)

                      tACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGAt                                                                     <  2:4181505/93‑1 (MQ=255)

                          cTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGGGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCtt                                                                 >  2:972983/1‑93 (MQ=255)

                           tCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTa                                                                >  2:1630189/1‑93 (MQ=255)

                           tCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGAACTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTa                                                                >  1:2988989/1‑93 (MQ=255)

                            cGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTaa                                                               <  1:59788/93‑1 (MQ=255)

                             gATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAAc                                                              >  1:1056212/1‑93 (MQ=255)

                             gATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAAc                                                              >  2:3894744/1‑93 (MQ=255)

                             gATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAAc                                                              >  2:4557672/1‑93 (MQ=255)

                             gATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGGGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAAc                                                              >  2:2594027/1‑93 (MQ=255)

                             gATGGTCGCGTGGATGTCCCTNGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAAc                                                              <  1:213192/93‑1 (MQ=255)

                              aTGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACg                                                             >  1:3134257/1‑93 (MQ=255)

                              aTGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCGGCGATGTGGTGATGCTTAACg                                                             <  2:3011037/93‑1 (MQ=255)

                               tGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGa                                                            <  1:766758/93‑1 (MQ=255)

                                 gTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATa                                                          >  2:858923/1‑93 (MQ=255)

                                  tCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATaa                                                         <  2:3803048/93‑1 (MQ=255)

                                     cgTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCt                                                      >  2:3421835/1‑93 (MQ=255)

                                     cgTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCt                                                      >  1:164176/1‑93 (MQ=255)

                                      gTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTg                                                     <  1:18974/93‑1 (MQ=255)

                                       tGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGc                                                    <  2:2299328/93‑1 (MQ=255)

                                       gggATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGc                                                    <  2:2348457/92‑1 (MQ=255)

                                        ggATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCa                                                   >  2:3850173/1‑93 (MQ=255)

                                        ggATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCa                                                   >  2:420958/1‑93 (MQ=255)

                                        ggATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCa                                                   >  2:613783/1‑93 (MQ=255)

                                          aTGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAAc                                                 <  2:1237347/93‑1 (MQ=255)

                                             tCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCgg                                              <  1:2400146/93‑1 (MQ=255)

                                              cccTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGa                                             <  1:4120972/93‑1 (MQ=255)

                                              cccTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGa                                             <  2:714631/93‑1 (MQ=255)

                                               ccTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGaa                                            >  2:319252/1‑93 (MQ=255)

                                               ccTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGaa                                            >  2:28236/1‑93 (MQ=255)

                                                cTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAag                                           <  2:4594405/93‑1 (MQ=255)

                                                cTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAag                                           <  2:3643767/93‑1 (MQ=255)

                                                cTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGGCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAag                                           >  2:3585615/1‑93 (MQ=255)

                                                 ttGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAaga                                          >  2:1623559/1‑93 (MQ=255)

                                                 ttGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAaga                                          >  1:1364783/1‑93 (MQ=255)

                                                 ttGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAaga                                          >  1:2336618/1‑93 (MQ=255)

                                                 ttGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAaga                                          >  1:2602719/1‑93 (MQ=255)

                                                 ttGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAaga                                          >  1:2858908/1‑93 (MQ=255)

                                                    ggCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGcc                                       >  2:1527667/1‑93 (MQ=255)

                                                    ggCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGcc                                       <  2:3059380/93‑1 (MQ=255)

                                                     gcgcgcATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGCAGAGCCg                                      >  1:4609412/1‑93 (MQ=255)

                                                     gcgcgcATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCg                                      >  1:4609316/1‑93 (MQ=255)

                                                      cgcgcATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGa                                     <  1:2535058/93‑1 (MQ=255)

                                                        cgcATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGaaa                                   >  1:151356/1‑93 (MQ=255)

                                                         gcATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGaaaa                                  <  2:2676370/93‑1 (MQ=255)

                                                         gcATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGaaaa                                  >  2:3744493/1‑93 (MQ=255)

                                                         gcATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGaaaa                                  >  2:1538560/1‑93 (MQ=255)

                                                          cATCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAAc                                 >  2:816298/1‑93 (MQ=255)

                                                           aTCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACg                                >  1:1549822/1‑93 (MQ=255)

                                                           aTCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACg                                >  2:449897/1‑93 (MQ=255)

                                                            tCGTGGTGCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACgg                               >  1:4536122/1‑93 (MQ=255)

                                                              gtggtgCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCg                             >  1:1127236/1‑93 (MQ=255)

                                                               tggtgCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGt                            <  2:2540754/93‑1 (MQ=255)

                                                               tggtgCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGt                            <  2:2226862/93‑1 (MQ=255)

                                                                ggtgCACGATCTACCAGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTc                           <  2:751953/93‑1 (MQ=255)

                                                                  tgCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGa                         >  1:1003480/1‑93 (MQ=255)

                                                                  tgCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGa                         <  1:2459735/93‑1 (MQ=255)

                                                                   gCACGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGAt                        >  1:2655064/1‑93 (MQ=255)

                                                                      cGATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTcc                     <  2:3029001/93‑1 (MQ=255)

                                                                       gATCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCg                    >  1:4021889/1‑93 (MQ=255)

                                                                        aTCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGa                   >  2:1812863/1‑93 (MQ=255)

                                                                         tCTGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGAt                  >  1:234206/1‑93 (MQ=255)

                                                                           tGCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTa                >  1:22667/1‑93 (MQ=255)

                                                                            gCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTaa               >  2:1558857/1‑93 (MQ=255)

                                                                            gCCGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTaa               <  2:2477441/93‑1 (MQ=255)

                                                                             ccGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAAc              >  2:2780273/1‑93 (MQ=255)

                                                                              cGGAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAACa             <  1:1685900/93‑1 (MQ=255)

                                                                               ggAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAACAg            >  2:4305948/1‑93 (MQ=255)

                                                                               ggAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAACAg            >  1:3500610/1‑93 (MQ=255)

                                                                               ggAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAACAg            <  1:1214512/93‑1 (MQ=255)

                                                                                gtaaGCGCAGGAGGCGTGTCGAGAAATGTCGTTATGGTAAACGAAACCCTGCATCCGGAAGAGCCGAAACCGGCGTCGATTCCGATTAACact           <  2:4258563/91‑3 (MQ=255)

                                                                                gAAAGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAACAGt           <  1:4248421/93‑1 (MQ=255)

                                                                                 aaaGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAACAGTa          >  2:332575/1‑93 (MQ=255)

                                                                                  aaGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAACAGTaa         <  1:4017747/93‑1 (MQ=255)

                                                                                  aaGCGCAGCAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAACAGTaa         <  1:4427339/93‑1 (MQ=255)

                                                                                       cagcagGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGa    <  2:2768632/93‑1 (MQ=255)

                                                                                         gTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGAtt  <  2:1995521/93‑1 (MQ=255)

TCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATT  >  NC_000913/4444600‑4444781

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