BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,688,218	Δ1 bp	6.1%	coding (292/1107 nt)	PP_3253 →	carboxylate‑amine ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,688,218	0	T	.	6.1%	398.8 / 19.8	132	coding (292/1107 nt)	PP_3253	carboxylate‑amine ligase
Reads supporting (aligned to +/- strand): ref base T (63/60); new base . (4/4); total (68/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.72e-01

Only 100 of 133 total aligned reads displayed.

GCCAGATCGAGCTGGCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTT  >  NC_002947/3688075‑3688367

    gATCGAGCTGGCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCa                                                                                                                                               <  2:1965870/150‑1 (MQ=255)

    gATCGAGCTGGCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCa                                                                                                                                               <  1:486956/150‑1 (MQ=255)

       cGAGCTGGCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGaa                                                                                                                                            <  1:506640/150‑1 (MQ=255)

            tGGCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑cg                                                                                                                                      >  1:2716002/1‑150 (MQ=255)

                 tCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGca                                                                                                                                 <  1:649136/150‑1 (MQ=255)

                  cGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCCGTGAAGCCTGGACTACAACTAC‑CG‑GGAAGTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGTGCCGCG‑GGCTTCTGCAGCAGCCACGCGGcgc                                                                                                                                <  1:1135338/150‑3 (MQ=255)

                   gCCGGTGTCCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGaa                                                                                                                                            <  1:1692998/138‑1 (MQ=255)

                      ggTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGTGCCGCG‑GGCTTCTGCAGCAGCCACGCGGcgcaggg                                                                                                                            <  1:2634992/150‑7 (MQ=38)

                      ggTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAgg                                                                                                                            <  2:2844785/150‑1 (MQ=255)

                        tgtTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGc                                                                                                                          >  1:2540860/1‑150 (MQ=255)

                          ttCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCAc                                                                                                                        >  2:3044303/1‑150 (MQ=255)

                             aCCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGgctgct                                                                                                                                                  >  2:774259/1‑122 (MQ=255)

                             aCCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGgctgct                                                                                                                                                  <  1:774259/122‑1 (MQ=255)

                             aCCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACtat                                                                                                                     <  2:1307775/150‑1 (MQ=255)

                                aaCCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAag                                                                                                                  >  2:2601777/1‑150 (MQ=255)

                                 aCCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAagc                                                                                                                 >  2:1341940/1‑150 (MQ=255)

                                   cTGTATGAAGCACGTGAGTTTCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAagcag                                                                                                               >  1:1956493/1‑150 (MQ=255)

                                    tGTATGAAGCACGTGAGTTTCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCTCTGCCTCGCTGGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAagcagc                                                                                                              >  1:2984943/1‑150 (MQ=255)

                                        tGAAGCACGTGAGTTTCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCAc                                                                                                                        >  1:2481336/1‑136 (MQ=255)

                                        tGAAGCACGTGAGTTTCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCAc                                                                                                                        <  2:2481336/136‑1 (MQ=255)

                                        tGAAGCACGTGAGTTTCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGtt                                                                                                          <  1:1155454/150‑1 (MQ=255)

                                        tGAAGCACGTGAGTTTCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGtt                                                                                                          <  1:1341940/150‑1 (MQ=255)

                                         gAAGCACGTGAGTTTCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATCCGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTc                                                                                                         >  1:1144228/1‑150 (MQ=255)

                                            gCACGTGAGTTTCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTcgac                                                                                                      <  2:68197/150‑1 (MQ=255)

                                            gCACGTGAGTTTCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTcgac                                                                                                      <  1:863284/150‑1 (MQ=255)

                                            gCACGTGAGTTTCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTcgac                                                                                                      <  2:2419202/150‑1 (MQ=255)

                                              aCGTGAGTTTCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATaa                                                                                                                   >  2:2167454/1‑135 (MQ=255)

                                              aCGTGAGTTTCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGAGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCAAAAGGGCACTATAAGCAGCTGTTccccga                                                                                                    >  2:1706539/1‑145 (MQ=255)

                                                gTGAGTTTCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACt                                                                                                 >  1:349869/1‑150 (MQ=255)

                                                 tGAGTTTCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTa                                                                                                >  1:1109898/1‑150 (MQ=255)

                                                  gAGTTTCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTAc                                                                                               >  1:1875213/1‑150 (MQ=255)

                                                   aGTTTCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTAcc                                                                                              >  1:2145054/1‑150 (MQ=255)

                                                    gTTTCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCg                                                                                             >  1:1424666/1‑150 (MQ=255)

                                                    gTTTCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCGGCAGAGGCC‑CGCGGCACAAGAGCACGATAAGCAGCTGTTCGACG‑ACTACCg                                                                                             >  2:2039201/1‑150 (MQ=255)

                                                     tttCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGc                                                                                            >  2:2401215/1‑150 (MQ=255)

                                                     tttCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGc                                                                                            >  2:2064196/1‑150 (MQ=255)

                                                     tttCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAACC‑CGCGGCACAAAGGCACTATAAGCAGCTGTTCTACG‑ACTACCGc                                                                                            >  2:501384/1‑150 (MQ=255)

                                                         tcCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGTGCCGCG‑GGCTTCTGCAGCAGCCACGCGGCGCAggggtgg                                                                                                                        >  1:1481949/1‑116 (MQ=25)

                                                         tcCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGTGCCGCG‑GGCTTCTGCAGCAGCCACGCGGCGCAggggtgg                                                                                                                        <  2:1481949/119‑4 (MQ=25)

                                                         tcCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAg                                                                                                                                              <  2:877407/98‑1 (MQ=255)

                                                         tcCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGAGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCAtg                                                                                        >  2:1646561/1‑150 (MQ=255)

                                                         tcCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCG‑GGCTTCTGCAGCAGccacgc                                                                                                                                     >  2:1606196/1‑102 (MQ=255)

                                                          ccAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCAtgt                                                                                       <  2:2540860/150‑1 (MQ=255)

                                                           cAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGcac                                                                                                                                <  1:2582915/109‑1 (MQ=255)

                                                           cAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCt                                                                                                             >  2:2686593/1‑128 (MQ=255)

                                                           cAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCt                                                                                                             <  1:2686593/128‑1 (MQ=255)

                                                             gcgcAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGc                                                                                    >  1:386805/1‑150 (MQ=255)

                                                                  aCC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGc                                                                               <  2:1887972/150‑1 (MQ=255)

                                                                   cc‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCTTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCg                                                                              >  2:2199635/1‑150 (MQ=255)

                                                                       gTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTa                                                                            >  2:941467/1‑150 (MQ=255)

                                                                        tCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAg                                                                           >  1:249585/1‑150 (MQ=255)

                                                                         tgacGGTTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGc                                                                          <  2:2791883/147‑1 (MQ=255)

                                                                         cAGCGGCTGCGCGTGGAGTTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGc                                                                          <  1:2111826/150‑1 (MQ=255)

                                                                          aGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGcc                                                                         >  2:431741/1‑150 (MQ=255)

                                                                            cGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCAc                                                                                                                        >  2:955159/1‑102 (MQ=255)

                                                                            cGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCAc                                                                                                                        <  1:955159/102‑1 (MQ=255)

                                                                              gCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCt                                                                     <  1:1665892/150‑1 (MQ=255)

                                                                                     gTGGAATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGc                                                              <  2:851456/150‑1 (MQ=255)

                                                                                      tGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCAtgt                                                                                       <  2:2868130/124‑1 (MQ=255)

                                                                                      tGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATg                                                                <  1:592208/147‑1 (MQ=255)

                                                                                      tGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATg                                                                >  2:592208/1‑147 (MQ=255)

                                                                                       ggCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCt                                                            >  2:2765138/1‑150 (MQ=255)

                                                                                         cATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGc                                                          >  2:2563027/1‑150 (MQ=255)

                                                                                          aTTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGgc                                                                                                                                                      <  2:2781729/59‑1 (MQ=255)

                                                                                          aTTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGgc                                                                                                                                                      >  1:2781729/1‑59 (MQ=255)

                                                                                          aTTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTa                                                                            >  1:638383/1‑131 (MQ=255)

                                                                                          aTTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCa                                                         <  1:2601777/150‑1 (MQ=255)

                                                                                           ttGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCAt                                                        <  1:2064196/150‑1 (MQ=255)

                                                                                             ggCCGAAGAGGGGATGGGGCCATAAGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGTCTGCAtgt                                                      <  1:8332/150‑1 (MQ=255)

                                                                                               ccGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCa                                                                                                                         <  1:513483/82‑1 (MQ=255)

                                                                                               ccGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCa                                                                                                                         >  2:513483/1‑82 (MQ=255)

                                                                                               ccGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGc                                                    >  1:469566/1‑150 (MQ=255)

                                                                                                  aaGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACg                                                 >  2:1769280/1‑150 (MQ=255)

                                                                                                      ggggATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCAtgt                                                                                       <  1:2909769/108‑1 (MQ=255)

                                                                                                      ggggATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCAtgt                                                                                       >  2:2909769/1‑108 (MQ=255)

                                                                                                       gggATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGc                                            >  2:1855171/1‑150 (MQ=255)

                                                                                                           tGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGGCG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTgc                                        <  1:236099/150‑1 (MQ=255)

                                                                                                           tGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGGCG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTgc                                        <  1:236059/150‑1 (MQ=255)

                                                                                                             gggCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTgccg                                      >  1:267977/1‑150 (MQ=255)

                                                                                                               gCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTgccgcc                                    >  2:428969/1‑150 (MQ=255)

                                                                                                               gCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTgccgcc                                    >  2:1500444/1‑150 (MQ=255)

                                                                                                                ccATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCa                                   >  2:2604466/1‑150 (MQ=255)

                                                                                                                    aCGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGTAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCtt                               <  1:2039201/150‑1 (MQ=255)

                                                                                                                    aCGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTg                                                                       <  2:1996589/110‑1 (MQ=255)

                                                                                                                    aCGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTg                                                                       >  1:1996589/1‑110 (MQ=255)

                                                                                                                      gCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTtgt                             <  1:2765138/150‑1 (MQ=255)

                                                                                                                      gCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTtgt                             <  1:431741/150‑1 (MQ=255)

                                                                                                                        tGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGa                           <  1:2253396/150‑1 (MQ=255)

                                                                                                                           ggCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTtgt                             <  1:1760765/145‑1 (MQ=255)

                                                                                                                           ggCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCg                        >  1:2779175/1‑150 (MQ=255)

                                                                                                                                 ccACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCAt                     >  1:3027157/1‑147 (MQ=255)

                                                                                                                                  cACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCaa                 <  1:948866/150‑1 (MQ=255)

                                                                                                                                     ccttgtgccgcG‑GGCTTCTGCAGCAGCCACGCGGCGC‑AGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGa             >  1:404791/7‑150 (MQ=37)

                                                                                                                                      cttgtgccgcG‑GGCTTCTGCAGCAGCCACGCGGCGCAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGa             >  2:1736895/6‑150 (MQ=38)

                                                                                                                                      cttGCGCCGCGTGGCTGCTGCAAAAGCC‑CGGGGCACAAGAGAACTCTACGCAGCAGTTGTCTGTTCTGTCG‑AAGTTGGCAAGGCTTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGa             <  1:55091/148‑1 (MQ=255)

                                                                                                                                      ccTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGa             <  1:1855171/150‑1 (MQ=255)

                                                                                                                                            ccgcGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATg                                                                >  2:1759363/1‑93 (MQ=255)

                                                                                                                                                gTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAAc        >  2:1688681/1‑145 (MQ=255)

                                                                                                                                                gTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAAc        <  1:1688681/145‑1 (MQ=255)

                                                                                                                                                 tGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCtt  <  2:2507382/150‑1 (MQ=255)

GCCAGATCGAGCTGGCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTT  >  NC_002947/3688075‑3688367

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