BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,688,230	A→C	5.8%	K102Q (AAG→CAG)	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,230	0	A	C	5.8%	385.1 / 8.7	132	K102Q (AAG→CAG)	PP_3253	carboxylate‑amine ligase
Reads supporting (aligned to +/- strand): ref base A (66/57); new base C (4/4); total (70/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 136 total aligned reads displayed.

CGAGCTGGCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTA  >  NC_002947/3688082‑3688379

     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)

            gCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGaa                                                                                                                                                        >  2:1692998/1‑138 (MQ=255)

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

             ccGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACaa                                                                                                                                          >  2:2253396/1‑150 (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)

                    tCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGTCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACt                                                                                                                                   >  2:2046222/1‑150 (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)

                          aCCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAagc                                                                                                                             <  1:2049703/150‑1 (MQ=255)

                             tGTATGAAGCACGTGAGTTTCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCTCTGCCTCGCTGGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAagcagc                                                                                                                          >  1:2984943/1‑150 (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)

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

                                          tGAGTTTCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGTTGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTa                                                                                                            >  1:1329036/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)

                                              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‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATCCGCTGCGGCAAGCCACCCCTGCGCCGCG‑GGCTTCTGCAGCAGccacgc                                                                                                                                                 <  1:1606196/106‑5 (MQ=38)

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

                                                  tcCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCAtg                                                                                                    >  1:1887972/1‑150 (MQ=255)

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

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

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

                                                    cAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGcac                                                                                                                                            >  2:2582915/1‑109 (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)

                                                                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)

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

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

                                                                         tGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCa                                                                               >  1:577502/1‑150 (MQ=255)

                                                                           cgcgTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAAt                                                                             <  2:249585/150‑1 (MQ=255)

                                                                           cgcgTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAAt                                                                             <  2:2716002/150‑1 (MQ=255)

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

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

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

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

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

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

                                                                                   aTTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCa                                                                     <  1:2601777/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                                                                >  2:1589886/1‑150 (MQ=255)

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

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

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

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

                                                                                                    tGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTgc                                                    <  2:469566/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)

                                                                                                             aCGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCtt                                           <  1:2401215/150‑1 (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)

                                                                                                                 tGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGa                                       <  1:501384/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                                 <  2:3027157/147‑1 (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                         <  2:1329036/150‑1 (MQ=255)

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

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

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

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

                                                                                                                                             ctgctgCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTg           <  1:1500444/150‑1 (MQ=255)

                                                                                                                                              tgctgcAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCTGAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGc          >  1:2518306/1‑150 (MQ=255)

                                                                                                                                                ctgcAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGc                                                                >  1:2112329/1‑94 (MQ=255)

                                                                                                                                                ctgcAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGc                                                                <  2:2112329/94‑1 (MQ=255)

                                                                                                                                                 ttcAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCt       >  1:2413269/3‑150 (MQ=255)

                                                                                                                                                   cAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGc                                                                                  <  1:2188004/73‑1 (MQ=255)

                                                                                                                                                   cAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGc                                                                                  >  2:2188004/1‑73 (MQ=255)

                                                                                                                                                    aGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGt    <  2:1828689/150‑1 (MQ=255)

                                                                                                                                                     gAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGtt   >  2:1799388/1‑150 (MQ=255)

                                                                                                                                                      aaGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTa  >  2:2178389/1‑150 (MQ=255)

CGAGCTGGCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACC‑‑GTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACG‑ACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTA  >  NC_002947/3688082‑3688379

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: `‑`