Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,634,508	T→G	13.1%	intergenic (+33/+350)	putA → / ← PP_5718	proline dehydrogenase/1‑pyrroline‑5‑carboxylate dehydrogenase/membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,634,508	0	T	G	13.1%	168.8 / 21.0	69	intergenic (+33/+350)	putA/PP_5718	proline dehydrogenase/1‑pyrroline‑5‑carboxylate dehydrogenase/membrane protein
Reads supporting (aligned to +/- strand):  ref base T (32/28);  new base G (5/4);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.80e-01

GTCCAGCGGGGATCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGCCTGGGCCGGTCC  >  NC_002947/5634367‑5634655                                                                                                                                              |                                                                                                                                                    gTCCAGCGGGGATCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCtctc                                                                                                                                             <  1:119278/150‑1 (MQ=255) gTCCAGCGGGGACCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGACTCGGCtctc                                                                                                                                             >  2:725298/1‑150 (MQ=255)      gCGGGGATCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCtttt                                                                                                                                        <  2:423082/150‑1 (MQ=255)      gCGGGGATCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCtttt                                                                                                                                        >  2:569226/1‑150 (MQ=255)           gATCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATg                                                                                                                                   >  2:1483418/1‑150 (MQ=255)           gATCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATg                                                                                                                                   >  1:979184/1‑150 (MQ=255)            aTCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGc                                                                                                                                  >  1:1589011/1‑150 (MQ=255)              cACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGt                                                                                                                                <  2:25181/150‑1 (MQ=255)                 caaaTTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCa                                                                                                                             >  1:1645306/4‑150 (MQ=255)                         gcTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTAt                                                                                                                                    >  2:560044/1‑135 (MQ=255)                         gcTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTAt                                                                                                                                    <  1:560044/135‑1 (MQ=255)                          cTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGaa                                                                                                                    <  1:1483418/150‑1 (MQ=255)                          cTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGaa                                                                                                                    <  1:725298/150‑1 (MQ=255)                                cgcCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAggg                                                                                                              <  1:911101/150‑1 (MQ=255)                                 gcCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAgggg                                                                                                             <  1:1085899/150‑1 (MQ=255)                                     ggTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCGAGGCtctc                                                                                                                                             >  2:570810/1‑114 (MQ=255)                                     ggTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCGAGGCtctc                                                                                                                                             <  1:570810/114‑1 (MQ=255)                                      gTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGa                                                                                                        >  1:1326494/1‑150 (MQ=255)                                             aGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAg                                                                                                 <  2:1188133/150‑1 (MQ=255)                                              gcgcgcGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAgg                                                                                                <  1:1541902/150‑1 (MQ=255)                                                        aGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCtttt                                                                                                                                        <  1:1275953/100‑1 (MQ=255)                                                        aGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCtttt                                                                                                                                        >  2:1275953/1‑100 (MQ=255)                                                         gCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGt                                                                                                                                <  1:1200044/107‑1 (MQ=255)                                                         gCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGt                                                                                                                                >  2:1200044/1‑107 (MQ=255)                                                         gCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGa                                                                                     >  1:1325119/1‑150 (MQ=255)                                                             gAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGt                                                                                                                                >  1:852003/1‑103 (MQ=255)                                                             gAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGt                                                                                                                                <  2:852003/103‑1 (MQ=255)                                                                    gCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGc                                                                                                                                  <  2:404297/94‑1 (MQ=255)                                                                    gCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGc                                                                                                                                  >  1:404297/1‑94 (MQ=255)                                                                     cAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGc                                                                         >  2:427088/1‑150 (MQ=255)                                                                      aGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCtctc                                                                                                                                             >  1:680881/1‑81 (MQ=255)                                                                      aGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCtctc                                                                                                                                             <  2:680881/81‑1 (MQ=255)                                                                          gctggtggtAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGCTGCTTGAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAGGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACg                                                                    >  1:861587/4‑150 (MQ=255)                                                                                     cGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGgt                                                         <  1:427088/150‑1 (MQ=255)                                                                                        cAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACt                                                                                                      >  2:995711/1‑102 (MQ=255)                                                                                        cAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACt                                                                                                      <  1:995711/102‑1 (MQ=255)                                                                                        cAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGGCTTGGGTGtc                                                      >  1:1486290/1‑150 (MQ=255)                                                                                        cAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGtc                                                      >  1:535119/1‑150 (MQ=255)                                                                                            cTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCAt                                                                                                                            >  2:577933/1‑76 (MQ=255)                                                                                            cTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCAt                                                                                                                            <  1:577933/76‑1 (MQ=255)                                                                                            cTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTgg                                                                                                         >  1:82747/1‑95 (MQ=255)                                                                                            cTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTgg                                                                                                         <  2:82747/95‑1 (MQ=255)                                                                                                   ccATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCAt                                                                           >  1:1342245/1‑118 (MQ=255)                                                                                                   ccATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCAt                                                                           <  2:1342245/118‑1 (MQ=255)                                                                                                    cATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCtttt                                                                                                                                        <  2:1514387/56‑1 (MQ=255)                                                                                                    cATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCtttt                                                                                                                                        >  1:1514387/1‑56 (MQ=255)                                                                                                    cATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGAGTGCATGTCCCGAAGGAGGAGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTAc                                          >  1:609423/1‑150 (MQ=255)                                                                                                      ttGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCg                                        >  1:945469/1‑150 (MQ=255)                                                                                                          cTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTa                                    >  2:1308578/1‑150 (MQ=255)                                                                                                           tAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTAt                                   >  1:696245/1‑150 (MQ=255)                                                                                                              tGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATccc                                <  1:1308578/150‑1 (MQ=255)                                                                                                              ggggttaGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATccc                                <  2:1326494/144‑1 (MQ=255)                                                                                                                           ttgAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCg                   <  2:696245/150‑1 (MQ=255)                                                                                                                                 aagaagGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGtc                                                      >  1:928522/1‑109 (MQ=255)                                                                                                                                 aagaagGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGtc                                                      <  2:928522/109‑1 (MQ=255)                                                                                                                                 aagaagGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGCGCTCACCAGAACGCGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATGCCCTGTAAATTCTCGGCTGcc             >  1:279318/1‑150 (MQ=255)                                                                                                                                  aaaagGAGAGCCGAGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGa                                                                                                        >  2:190046/3‑58 (MQ=25)                                                                                                                                  aaaagGAGAGCCGAGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGa                                                                                                        <  1:190046/56‑1 (MQ=37)                                                                                                                                   gaagGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGCCTg           <  2:1645306/150‑1 (MQ=255)                                                                                                                                   aaagGAGAGCCGAGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGt                                                                  <  2:707104/94‑1 (MQ=255)                                                                                                                                   aaagGAGAGCCGAGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGt                                                                  >  1:707104/2‑95 (MQ=255)                                                                                                                                    aagGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGgt                                                         >  1:21933/1‑103 (MQ=255)                                                                                                                                    aagGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGgt                                                         <  2:21933/103‑1 (MQ=255)                                                                                                                                    aagGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGCCTgg          >  2:107781/1‑150 (MQ=255)                                                                                                                                      ggAGAGCCTAGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGCCTGGGc        >  1:1630377/1‑150 (MQ=255)                                                                                                                                      ggAGAGCCGAGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGCCTGGGc        >  2:605604/1‑150 (MQ=255)                                                                                                                                       gagaGCCGAGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGTCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATtctc                    <  1:17290/137‑1 (MQ=255)                                                                                                                                       gagaGCCGAGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATtctc                    >  2:17290/1‑137 (MQ=255)                                                                                                                                            ccTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGCCTGGGCCGGTcc  <  2:979184/150‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                    GTCCAGCGGGGATCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGCCTGGGCCGGTCC  >  NC_002947/5634367‑5634655

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