Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,338,538	A→T	18.0%	intergenic (+114/+28)	PP_3812 → / ← PP_5632	lipase/allo‑isoleucine metabolism protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,338,538	0	A	T	18.0%	198.3 / 35.2	84	intergenic (+114/+28)	PP_3812/PP_5632	lipase/allo‑isoleucine metabolism protein
Reads supporting (aligned to +/- strand):  ref base A (38/30);  new base T (7/8);  total (46/38)
Fisher's exact test for biased strand distribution p-value = 5.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

GGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATCGA  >  NC_002947/4338394‑4338684                                                                                                                                                 |                                                                                                                                                   gggCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGGCCGGCCCTTTCGCGGGTTGACCcg                                                                                                                                               >  1:1646593/1‑150 (MQ=255)  ggCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCGAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCcgc                                                                                                                                              <  1:1853191/150‑1 (MQ=255)    cAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGa                                                                                                                                            >  1:1529293/1‑150 (MQ=255)      aCGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTTTACCCGCGaaa                                                                                                                                          <  2:684819/150‑1 (MQ=255)          aCTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTTTACCCGCGa                                                                                                                                            <  2:1261095/144‑1 (MQ=255)          aCTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTTTACCCGCGa                                                                                                                                            >  1:1261095/1‑144 (MQ=255)              aaaGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTTTACCCGCGAAAGGGcc                                                                                                                                     >  1:1212239/1‑147 (MQ=255)              aaaGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTTTACCCGCGAAAGGGcc                                                                                                                                     <  2:1212239/147‑1 (MQ=255)                 gTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGGCCGGCCCTTTCGCGGGTTTACCCGCGAAAGGGCCggacct                                                                                                                               <  2:1646593/150‑5 (MQ=255)                      tCTGCAGTGAGGCGATGAGCAATAGACCTGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTg                                                                                                                          <  1:183269/150‑1 (MQ=255)                      tCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTg                                                                                                                          >  1:1281147/1‑150 (MQ=255)                       cggCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGc                                                                                                                         <  1:241881/148‑1 (MQ=255)                                   gATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTTTACCCGCGAAAGGGCCg                                                                                                                                    >  2:291895/1‑127 (MQ=255)                                   gATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTTTACCCGCGAAAGGGCCg                                                                                                                                    <  1:291895/127‑1 (MQ=255)                                    agaaGCAAAAGACCCGCCCAGGCGGGTCTATTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGc                                                                                                                                      >  1:1629075/4‑124 (MQ=38)                                    agaaGCAAAAGACCCGCCCAGGCGGGTCTATTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGc                                                                                                                                      <  2:1629075/121‑1 (MQ=38)                                    agaaGCAAAAGACCCGCCCAGGCGGGTCTATTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGc                                                                                                                                      <  2:1628273/121‑1 (MQ=38)                                    agaaGCAAAAGACCCGCCCAGGCGGGTCTATTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGc                                                                                                                                      >  1:1628273/4‑124 (MQ=38)                                    aTGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAAc                                                                                                            >  2:1154753/1‑150 (MQ=255)                                       aGCAAAAGACCCGCCCAGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGaaa                                                                                                                                          <  2:1137818/117‑1 (MQ=39)                                       aGCAAAAGACCCGCCCAGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGaaa                                                                                                                                          >  1:1137818/1‑117 (MQ=39)                                        gCAAAAGACCCGCCCAGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAg                                                                                                                     >  1:1025577/1‑137 (MQ=255)                                        gCAAAAGACCCGCCCAGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAg                                                                                                                     <  2:1025577/137‑1 (MQ=255)                                         cAATAGACCCGCCTGGGCGGGACTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCGGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCACATCCCGCTCAGGGCTCAACCCCtcc                                                                                                       >  1:696656/1‑148 (MQ=255)                                           aaaGACCCGCCCAAGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCca                                                                                                     >  1:396392/3‑150 (MQ=255)                                              gACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACaa                                                                                                  <  2:1764519/150‑1 (MQ=255)                                               aCCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCgg                                                                                                                                   <  1:1012712/116‑1 (MQ=255)                                               aCCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCgg                                                                                                                                   >  2:1012712/1‑116 (MQ=255)                                               aCCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAAt                                                                                                 >  2:1714923/1‑150 (MQ=255)                                               aCCCGCCCAAGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAAt                                                                                                 <  2:396392/150‑1 (MQ=255)                                                 ccGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCgg                                                                                                                                   <  2:1178261/114‑1 (MQ=255)                                                 ccGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCgg                                                                                                                                   >  1:1178261/1‑114 (MQ=255)                                                 ccGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGa                                                                                               >  2:1490528/1‑150 (MQ=255)                                                   gCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATc                                                                                             <  2:1586714/150‑1 (MQ=255)                                                            ggTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCt                                                                                                         >  2:790266/1‑129 (MQ=255)                                                            ggTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCt                                                                                                         <  1:790266/129‑1 (MQ=255)                                                            ggTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCAc                                                                                           <  1:402986/143‑1 (MQ=255)                                                            ggTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCAc                                                                                           >  2:402986/1‑143 (MQ=255)                                                              tCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAggg                                                                                                                   >  1:133400/1‑117 (MQ=255)                                                              tCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAggg                                                                                                                   <  2:133400/117‑1 (MQ=255)                                                                ttttGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGaaa                                                                                                                                          <  2:1111108/92‑1 (MQ=255)                                                                ttttGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGaaa                                                                                                                                          >  1:1111108/1‑92 (MQ=255)                                                                 tttGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAAGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGc                                                                               <  1:344253/150‑1 (MQ=255)                                                                    gCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCAc                                                                            <  2:292032/150‑1 (MQ=255)                                                                     cTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCAcc                                                                           >  1:1042733/1‑150 (MQ=255)                                                                      ttCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCAccc                                                                          >  1:1764595/1‑150 (MQ=255)                                                                      ttCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAACCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCAccc                                                                          >  2:1649472/1‑150 (MQ=255)                                                                                ggAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGaaa                                                                                                                                          <  2:1351127/76‑1 (MQ=255)                                                                                ggAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGaaa                                                                                                                                          >  1:1351127/1‑76 (MQ=255)                                                                                      aGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCa                                                          >  2:29474/1‑150 (MQ=255)                                                                                      aGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTCCTCAGCGCTCAAACCCTTCCACCATGCTCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGcc                                                          >  1:729674/1‑149 (MQ=255)                                                                                       gCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCATAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCAt                                                         <  2:1529293/150‑1 (MQ=255)                                                                                       gCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCAt                                                         <  1:1490528/150‑1 (MQ=255)                                                                                        cccGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATg                                                        <  2:1684438/150‑1 (MQ=255)                                                                                         ccGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAAGAACGCGCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGc                                                       >  1:1304282/1‑150 (MQ=255)                                                                                                 ttCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGt                                               >  1:1661693/1‑150 (MQ=255)                                                                                                 ttCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCA                                                                                       <  1:412256‑M1/110‑1 (MQ=255)                                                                                                 ttCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCA                                                                                       >  2:412256‑M1/1‑110 (MQ=255)                                                                                                  tCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTa                                              <  1:41951/150‑1 (MQ=255)                                                                                                  tCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTa                                              >  2:1785792/1‑150 (MQ=255)                                                                                                     aTGAATCTGCCAGCGGTACGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGCGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAg                                           >  2:2836/1‑150 (MQ=255)                                                                                                         aTCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCATCGCCAAGGGGTGGCGAAggctgaggtgatcattgtggaagggt                                       >  1:159755‑M2/1‑124 (MQ=255)                                                                                                            tGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAg                                    <  2:1281147/150‑1 (MQ=255)                                                                                                            tGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCGCAGCCTTCGCCTCCCCTTGTCGACCGCTGCTTGCCTCGTTGTACAGCGCGGAg                                    >  2:1069730/1‑150 (MQ=255)                                                                                                             gCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGc                                   >  2:1491393/1‑150 (MQ=255)                                                                                                                    gTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACa                            >  1:394641/1‑150 (MQ=255)                                                                                                                        ggCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATc                                                                                             <  1:221062/81‑1 (MQ=255)                                                                                                                        ggCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATc                                                                                             >  2:221062/1‑81 (MQ=255)                                                                                                                        ggCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCAcc                                                                                          <  2:215391/84‑1 (MQ=255)                                                                                                                        ggCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCAcc                                                                                          >  1:215391/1‑84 (MQ=255)                                                                                                                                 ggCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTgc               >  2:26195/1‑150 (MQ=255)                                                                                                                                   cccTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTgcgc             <  1:26195/150‑1 (MQ=255)                                                                                                                                   cccTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAACCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTgcgc             <  1:1649472/150‑1 (MQ=255)                                                                                                                                    ccTTTCGCGGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGcccg            >  1:774515/1‑150 (MQ=255)                                                                                                                                    ccTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCATTCCCCAATGAGCACCTCAGCCTTCGCCACCCCTTGGCGATGGCGGCATGCCGCTTGGTACAGCGCGGAGCTGTGACAGGCGACGGCCTGcgcg            >  1:665357/1‑150 (MQ=255)                                                                                                                                     cTTTCGCGGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGt           >  1:763914/1‑150 (MQ=255)                                                                                                                                      tttCGCGGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAg                                    >  1:1233141/1‑124 (MQ=255)                                                                                                                                      tttCGCGGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAg                                    <  2:1233141/124‑1 (MQ=255)                                                                                                                                      tttCGCGGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTa          <  2:763914/150‑1 (MQ=255)                                                                                                                                      tttCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTa          >  2:28385/1‑150 (MQ=255)                                                                                                                                       ttCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTAc         <  2:1661693/150‑1 (MQ=255)                                                                                                                                        tCGCGGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTgc               >  2:1614980/1‑143 (MQ=255)                                                                                                                                        tCGCGGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTgc               <  1:1614980/143‑1 (MQ=255)                                                                                                                                              gTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATCGa  >  1:1223166/1‑150 (MQ=255)                                                                                                                                                 |                                                                                                                                                   GGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATCGA  >  NC_002947/4338394‑4338684

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