BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,025,101	T→C	A318A (GCA→GCG)	b1956 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,025,101	0	T	C	96.5%	771.7 / 10.1	228	A318A (GCA→GCG)	b1956	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (3/5); new base C (106/114); total (109/119)
Fisher's exact test for biased strand distribution p-value = 7.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 231 total aligned reads displayed.

ACGCACCAGCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACT  >  NC_000913/2025011‑2025190

   cACCAGCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTg                                                                                      <  2:5997811/93‑1 (MQ=255)

     ccAGCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTgcg                                                                                    >  2:5354950/1‑93 (MQ=255)

      cAGCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTgcgc                                                                                   >  2:679176/1‑93 (MQ=255)

      cAGCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTgcgc                                                                                   >  1:3824692/1‑93 (MQ=255)

        gCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCa                                                                                 <  2:2206420/93‑1 (MQ=255)

         cacaCCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAg                                                                                <  2:3479751/93‑1 (MQ=255)

           cacCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCa                                                                              >  2:5270529/1‑93 (MQ=255)

            cccgccgAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCaa                                                                             <  1:1601609/92‑1 (MQ=255)

            cccGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCaa                                                                             <  1:469419/92‑1 (MQ=255)

             ccGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAg                                                                            >  2:1666044/1‑93 (MQ=255)

               gTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGc                                                                          >  2:4310457/1‑93 (MQ=255)

                tCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGcc                                                                         >  1:5027668/1‑93 (MQ=255)

                tCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGcc                                                                         >  1:224729/1‑93 (MQ=255)

                tCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGcc                                                                         <  1:2671017/93‑1 (MQ=255)

                  gAAATGGTCCAGACGTTCCCAGCTAACATAGCTACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAg                                                                       <  1:3224973/93‑1 (MQ=255)

                  gAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAg                                                                       <  1:5758555/93‑1 (MQ=255)

                   aaaTGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGt                                                                      >  1:4306063/1‑93 (MQ=255)

                    aaTGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGtt                                                                     <  2:15325/93‑1 (MQ=255)

                       ggTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCAc                                                                  >  2:3435903/1‑93 (MQ=255)

                        gTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACg                                                                 >  2:2217849/1‑93 (MQ=255)

                           cAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAgg                                                              >  2:6111599/1‑93 (MQ=255)

                           cAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAgg                                                              >  1:1786943/1‑93 (MQ=255)

                            aggcGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGa                                                             <  2:2124429/90‑1 (MQ=255)

                            aGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGa                                                             <  1:2753340/93‑1 (MQ=255)

                            aGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGa                                                             <  2:131006/93‑1 (MQ=255)

                            aGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGa                                                             <  1:3421403/93‑1 (MQ=255)

                            aGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGa                                                             <  1:6140816/93‑1 (MQ=255)

                            aGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGa                                                             <  1:2086445/93‑1 (MQ=255)

                            aGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGa                                                             <  1:5889867/93‑1 (MQ=255)

                             gACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGAt                                                            >  2:1745979/1‑93 (MQ=255)

                             gACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGAt                                                            <  1:1987492/93‑1 (MQ=255)

                              aCGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATc                                                           >  1:3290390/1‑93 (MQ=255)

                                 ttCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCaaa                                                        >  2:4871575/1‑93 (MQ=255)

                                 ttCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCaaa                                                        <  2:63179/93‑1 (MQ=255)

                                  tCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCaaaa                                                       <  1:3702187/93‑1 (MQ=255)

                                   cccAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCaaaaa                                                      >  2:5291649/1‑93 (MQ=255)

                                   cccAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCaaaaa                                                      <  1:2022619/93‑1 (MQ=255)

                                   cccAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCaaaaa                                                      <  1:5417906/93‑1 (MQ=255)

                                     cAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAata                                                    <  2:304922/93‑1 (MQ=255)

                                      aGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAatat                                                   >  1:1046497/1‑93 (MQ=255)

                                      aGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAatat                                                   >  1:2291368/1‑93 (MQ=255)

                                       gCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGTTCAAAAATAtt                                                  <  1:3664514/93‑1 (MQ=255)

                                       gCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATAtt                                                  >  2:3315326/1‑93 (MQ=255)

                                        cTAACATAGCGCCTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTc                                                 <  1:2567357/93‑1 (MQ=255)

                                        cTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTc                                                 <  2:5946122/93‑1 (MQ=255)

                                        cTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTc                                                 <  2:1788548/93‑1 (MQ=255)

                                         tAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTcc                                                >  2:2237717/1‑93 (MQ=255)

                                         tAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTcc                                                >  1:3025227/1‑93 (MQ=255)

                                          aaCATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTccc                                               <  1:4269369/93‑1 (MQ=255)

                                          aaCATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTccc                                               >  1:1287632/1‑93 (MQ=255)

                                           aCATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCt                                              <  2:3886452/93‑1 (MQ=255)

                                            cATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCACGGCCAGTTCACGAGGATCAAAAATATTCCCTg                                             >  1:1446857/1‑93 (MQ=255)

                                            cATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTg                                             <  2:2342632/93‑1 (MQ=255)

                                             aTAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGt                                            <  1:518033/93‑1 (MQ=255)

                                              tAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGtt                                           >  2:5781429/1‑93 (MQ=255)

                                              aaGCGACTGGCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGtt                                           <  1:1730526/92‑1 (MQ=255)

                                               aGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTg                                          <  1:1052832/93‑1 (MQ=255)

                                                 cGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGa                                        <  1:5509314/93‑1 (MQ=255)

                                                  gACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGAt                                       >  2:3820848/1‑93 (MQ=255)

                                                   aCTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGATGATCAAAAATATTCCCTGTTGGATg                                      >  2:4445894/1‑93 (MQ=255)

                                                   aCTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGTTCAAAAATATTCCCTGTTGGATg                                      >  1:5656797/1‑93 (MQ=255)

                                                   aCTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATg                                      >  1:4883592/1‑93 (MQ=255)

                                                     tGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGt                                    <  2:1508600/93‑1 (MQ=255)

                                                     tGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGt                                    >  2:4181188/1‑93 (MQ=255)

                                                        ccATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAg                                 >  1:423305/1‑93 (MQ=255)

                                                         cATACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAgg                                >  2:2032653/1‑93 (MQ=255)

                                                          aTACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAggt                               >  2:4375662/1‑92 (MQ=255)

                                                           tACGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGAt                              >  2:4090222/1‑93 (MQ=255)

                                                           aaCGAATGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGAt                              <  1:4508139/92‑1 (MQ=255)

                                                               aaTGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGa                          >  2:4053945/1‑93 (MQ=255)

                                                                aTGCCGCCCCGCTTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGaa                         <  1:5545514/93‑1 (MQ=255)

                                                                aTGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGaa                         <  1:4653574/93‑1 (MQ=255)

                                                                 tGCCGCCCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAag                        <  1:5354050/93‑1 (MQ=255)

                                                                  gccgccCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAaga                       <  1:2925055/93‑1 (MQ=255)

                                                                   ccgccCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAagag                      <  2:2569395/93‑1 (MQ=255)

                                                                     gccCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAgga                    >  2:4894817/1‑92 (MQ=255)

                                                                     gccCCGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGt                    >  1:2958158/1‑93 (MQ=255)

                                                                      ccccGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTc                   >  2:5550781/1‑93 (MQ=255)

                                                                       cccGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCa                  >  2:3672791/1‑93 (MQ=255)

                                                                       cccGCGTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCa                  >  2:3129961/1‑93 (MQ=255)

                                                                          gcgTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCaaaa               <  2:5629897/93‑1 (MQ=255)

                                                                           cgTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCaaaaa              <  2:6128422/93‑1 (MQ=255)

                                                                           cgTGTCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCaaaaa              >  1:4750500/1‑93 (MQ=255)

                                                                            gtgtCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAAt             >  2:4815189/1‑93 (MQ=255)

                                                                            gtgtCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAAt             >  1:4994906/1‑93 (MQ=255)

                                                                            gtgtCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAAt             >  2:3542157/1‑93 (MQ=255)

                                                                             tgtCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAAtc            >  2:4679144/1‑93 (MQ=255)

                                                                             tgtCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAAtc            <  2:2814155/93‑1 (MQ=255)

                                                                             tgtCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAAtc            <  1:4426803/93‑1 (MQ=255)

                                                                             tgtCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAAtc            >  2:1562228/1‑93 (MQ=255)

                                                                               tCATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAAtctc          <  2:3770138/93‑1 (MQ=255)

                                                                                cATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCa         <  1:937772/93‑1 (MQ=255)

                                                                                cATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCa         <  1:4851385/93‑1 (MQ=255)

                                                                                cATGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCa         <  1:4097262/93‑1 (MQ=255)

                                                                                 aTGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAg        >  1:285932/1‑93 (MQ=255)

                                                                                  tGTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGc       >  1:4463048/1‑93 (MQ=255)

                                                                                   gTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCt      <  2:3493936/93‑1 (MQ=255)

                                                                                   gTTCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCt      <  1:1330738/93‑1 (MQ=255)

                                                                                    ttCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCtt     >  2:2954123/1‑93 (MQ=255)

                                                                                    ttCCATCGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCtt     <  2:110091/93‑1 (MQ=255)

ACGCACCAGCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACT  >  NC_000913/2025011‑2025190

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