Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,315,017	T→G	10.0%	V392G (GTA→GGA)	aspC →	aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,315,017	0	T	G	10.0%	154.2 / 10.0	60	V392G (GTA→GGA)	aspC	aminotransferase
Reads supporting (aligned to +/- strand):  ref base T (26/28);  new base G (3/3);  total (29/31)
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.90e-01

CAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGATTCCCTATCAGGTGGCCGGGCGTCATCAGGTACGCTTTTCCGTGGCGTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGAGCAGCAGTTTCCCTATTTCGAGGTGGAGGACTTCGAGCGCAAGCGACGCGAATATGGCCTTCCTGATGCACCGCAGCT  >  NC_002947/4314868‑4315166                                                                                                                                                       |                                                                                                                                                      cAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGATTCCCTATCAGGTGGCCGGGCGTCATCAGGTACGCTTTTCCGTGGCGTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGt                                                                                                                                                       >  1:468345/1‑150 (MQ=255)           tCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGATTCCCTATCAGGTGGCCGGGCGTCATCAGGTACGCTTTTCCGTGGCGTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGc                                                                                                                                             <  2:897430/150‑1 (MQ=255)           tCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGATTCCCTATCAGGTGGCCGGGCGTCATCAGGTACGCTTTTCCGTGGCGTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGc                                                                                                                                             <  2:432473/150‑1 (MQ=255)              cAATACCTGGCCACTGAGCTGGGTATCCTGACGATTCCCTATCAGGTGGCCGGGCGTCATCAGGTACGCTTTTCCGTGGCGTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCtt                                                                                                                                          >  1:65045/1‑150 (MQ=255)                   ccTGGCCACTGAGCTGGGTATCCTGACGATTCCCTATCAGGTGGCCGGGCGTCATCAGGTACGCTTTTCCGTGGCGTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCga                                                                                                                                     <  2:468345/150‑1 (MQ=255)                     tGGCCACTGAGCTGGGTATCCTGACGATTCCCTATCAGGTGGCCGGGCGTCATCAGGTACGCTTTTCCGTGGCGTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCgaga                                                                                                                                   <  2:1004634/150‑1 (MQ=255)                      ggCCACTGAGCTGGGTATCCTGACGATTCCCTATCAGGTGGCCGGGCGTCATCAGGTACGCTTTTCCGTGGCGTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCGagag                                                                                                                                  <  2:831785/150‑1 (MQ=255)                      ggCCACTGAGCTGGGTATCCTGACGATTCCCTATCAGGTGGCCGGGCGTCATCAGGTACGCTTTTCCGTGGCGTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCGagag                                                                                                                                  <  2:65045/150‑1 (MQ=255)                       gCCACTGAGCTGGGTATCCTGACGATTCCCTATCAGGTGGCCGGGCGTCATCAGGTACGCTTTTCCGTGGCGTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCGAGAgg                                                                                                                                 <  2:962158/150‑1 (MQ=255)                           cTGAGCTGGGTATCCTGACGATTCCCTATCAGGTGGCCGGGCGTCATCAGGTACGCTTTTCCGTGGCGTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCGAGAGGGTAt                                                                                                                             <  2:1120166/150‑1 (MQ=255)                                cTGGGTATCCTGACGATTCCCTATCAGGTGGCCGGGCGTCATCAGGTACGCTTTTCCGTGGCGTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCgag                                                                                                                                    >  2:40115/1‑138 (MQ=255)                                cTGGGTATCCTGACGATTCCCTATCAGGTGGCCGGGCGTCATCAGGTACGCTTTTCCGTGGCGTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCgag                                                                                                                                    <  1:40115/138‑1 (MQ=255)                                        ccTGACGATTCCCTATCAGGTGGCCGGGCGTCATCAGGTACGCTTTTCCGTGGCGTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCACCTGTCCAATGTACCGTTCAGCCTTGCCGAGAGGGTATGGGCATGACTGtc                                                                                                                <  2:764528/150‑1 (MQ=255)                                          tGACGATTCCCTATCAGGTGGCCGGGCGTCATCAGGTACGCTTTTCCGTGGCGTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGcc                                                                                                                                       <  1:802373/125‑1 (MQ=255)                                          tGACGATTCCCTATCAGGTGGCCGGGCGTCATCAGGTACGCTTTTCCGTGGCGTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGcc                                                                                                                                       >  2:802373/1‑125 (MQ=255)                                          tGACGATTCCCTATCAGGTGGCCGGGCGTCATCAGGTACGCTTTTCCGTGGCGTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGtctc                                                                                                              >  2:601743/1‑150 (MQ=255)                                           gACGATTCCCTATCAGGTGGCCGGGCGTCATCAGGTACGCTTTTCCGTGGCGTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGtctct                                                                                                             >  1:1126221/1‑150 (MQ=255)                                           gACGATTCCCTATCAGGTGGCCGGGCGTCATCAGGTACGCTTTTCCGTGGCGTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGtctct                                                                                                             >  2:1300126/1‑150 (MQ=255)                                                    cTATCAGGTGGCCGGGCGTCATCAGGTACGCTTTTCCGTGGCGTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCGagag                                                                                                                                  <  1:1339445/120‑1 (MQ=255)                                                    cTATCAGGTGGCCGGGCGTCATCAGGTACGCTTTTCCGTGGCGTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCGagag                                                                                                                                  >  2:1339445/1‑120 (MQ=255)                                                          ggTGGCCGGGCGTCATCAGGTACGCTTTTCCGTGGCGTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCGACTGTCCAATGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGAc                                                                                              >  1:1297887/1‑150 (MQ=255)                                                           gTGGCCGGGCGTCATCAGGTACGCTTTTCCGTGGCGTTGTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACa                                                                                             >  1:30800/1‑150 (MQ=255)                                                             ggCCGGGCGTCCTTAGGTACGCTTTTCCGTGGCGTTCTGCGGCGAAGCCGAAAAAGT‑CGCAGCCTTTTCAGTGCAGCAACCGTCCAATGTACAGTTCAGCCTTGCCGAGGGTGTATGGGCATGACTGTCTCTCTTCGCAATGTCGACAgt                                                                                           <  1:324676/150‑1 (MQ=255)                                                                cGGGCGTCATCAGGTACGCTTTTCCGTGGCGTTCTGCGGCGAAGCCGAAACCGTGTTCGACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTc                                                                                        >  1:464129/1‑150 (MQ=255)                                                                     gTCATCAGGTACGCTTTTCCGTGGCGTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTc                                                                                     >  1:601416/1‑148 (MQ=255)                                                                     gTCATCAGGTACGCTTTTCCGTGGCGTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTc                                                                                     <  2:601416/148‑1 (MQ=255)                                                                        atcaGGTACGCTTTTCCGTGGCGTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGa                                                                                >  1:1331675/1‑150 (MQ=255)                                                                           aGGTACGCTTTTCCGTGGCGTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGagca                                                                             >  2:1091477/1‑150 (MQ=255)                                                                           aGGTACGCTTTTCCGTGGCGTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGagca                                                                             >  1:496978/1‑150 (MQ=255)                                                                              tACGCTTTTCCGTGGCGTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGAgcagca                                                                          <  2:651314/150‑1 (MQ=255)                                                                                  cTTTTCCGTGGCGTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGAGCAGCAGttt                                                                      <  1:1300126/150‑1 (MQ=255)                                                                                   ttttCCGTGGCGTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGAGCAGCAGTTTc                                                                     >  1:454679/1‑150 (MQ=255)                                                                                        cGTGGCGTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGAGCAGCAGTTTCCCTAt                                                                >  2:627360/1‑150 (MQ=255)                                                                                              gTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGAGCAGCAGTTTCCCTATTTCGAg                                                          <  1:1091477/150‑1 (MQ=255)                                                                                              gTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGAGCAGCAGTTTCCCTATTTCGAg                                                          <  1:481209/150‑1 (MQ=255)                                                                                               ttCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGAGCAGCAGTTTCCCTATTTCGAgg                                                         >  1:511207/1‑150 (MQ=255)                                                                                                 cTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGAGCAGCAGTTTCCCTATTTCGAGGTg                                                       <  2:1126221/150‑1 (MQ=255)                                                                                                 cTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGAGCAGCAGTTTCCCTATTTCGAGGTg                                                       <  2:454679/150‑1 (MQ=255)                                                                                                  tGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCGAGGGGGTATCGGCATGACTGTCTCTCTTTGCAAGGTTTACAGTGTCCTCGGCGAGCAGCAGTTTCCCTATTTCGAGGTgg                                                      >  1:1386760/1‑150 (MQ=255)                                                                                                   gcggcgAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGAGCAGCAGTTTCCCTATTTCGAGGTgga                                                     <  2:1331675/150‑1 (MQ=255)                                                                                                          aGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGAGCAGCAGTTTCCCTATTTCGAGGTGGAGGACTTc                                              <  2:30800/150‑1 (MQ=255)                                                                                                           gCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCtt                                                                                                                                          >  1:890626/1‑57 (MQ=255)                                                                                                           gCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCtt                                                                                                                                          <  2:890626/57‑1 (MQ=255)                                                                                                            ccGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGAGCAGCAGTTTCCCTATTTCGAGGTGGAGGACTTCGa                                            >  2:793120/1‑150 (MQ=255)                                                                                                              gAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGGCAGTGTCCTCGGCGAGCAGCAGTTTCCCTATTTCGAGGTGGAGGACTTCGAgc                                          >  2:1004791/1‑150 (MQ=255)                                                                                                              gAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGAGCAGCAGTTTCTCTATTTCGAGGTGGAGGACTTCGAgc                                          >  1:204641/1‑150 (MQ=255)                                                                                                                          aaCCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCGAAAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGAGCAGCAGTTTCCCTATTTCGAGGGGGAGGAATTCGGGCGCAAGCGAcgcg                              >  1:477008/1‑150 (MQ=255)                                                                                                                                  gCGT‑CAGCGACTGTCCAATGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGAGCAGCAGTTTCCCTATTTCGAGGTGGAGGACTTCGAGCGCAAGCGACGCGAATATGGc                      <  2:1297887/150‑1 (MQ=255)                                                                                                                                    cagcAGCAACTGTCCAATGTACAGTTCAGCGTTGCCGAGAGGATATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGAGTCCTCGGCGAGCAGCAGTTTCCCTATTTCGAGGTGGAGGACTTCGAGCGCAAGCGACGCGAATATGGcc                     <  2:477008/147‑1 (MQ=255)                                                                                                                                         tgaaCTGTACATTGGACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGAGCAGCAGTTTCCCTATTTCGAGGTgga                                                     >  2:238887/3‑113 (MQ=255)                                                                                                                                         tgaaCTGTACATTGGACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGAGCAGCAGTTTCCCTATTTCGAGGTgga                                                     <  1:238887/111‑1 (MQ=255)                                                                                                                                         tgaaCTGTACATTGGACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGAGCAGCAGTTTCCCTATTTCGAGGTGGAGGACTTCGAGCGCAAGCGACGCGAATATg                        <  1:150208/140‑1 (MQ=255)                                                                                                                                         tgaaCTGTACATTGGACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGAGCAGCAGTTTCCCTATTTCGAGGTGGAGGACTTCGAGCGCAAGCGACGCGAATATg                        >  2:150208/3‑142 (MQ=255)                                                                                                                                         tgaaCTGTACATTGGACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGAGCAGCAGTTTCCCTATTTCGAGGTGGAGGACTTCGAGCGCAAGCGACGCGAATATGGCCTTcc                 <  1:1189221/147‑1 (MQ=255)                                                                                                                                         tgaaCTGTACATTGGACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGAGCAGCAGTTTCCCTATTTCGAGGTGGAGGACTTCGAGCGCAAGCGACGCGAATATGGCCTTcc                 >  2:1189221/3‑149 (MQ=255)                                                                                                                                           aaCTGTCCAATGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGAGCAGCAGTTTCCCTATTTCGAGGTGGAGGACTTCGAGCGCAAGCGACGCGAATATGGCCTTCCTGa              <  2:464129/150‑1 (MQ=255)                                                                                                                                            aCTGTCCAATGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGAGCAGCAGTTTCCCTATTTCGAGGTGGAGGACTTCGAGCGCAAGCGACGCGAATATGGCCTTCCTGAt             <  1:290446/150‑1 (MQ=255)                                                                                                                                               gTCCAATGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGAGCAGCAGTTTCCCTATTTCGAGGTggagga                                                  >  1:608642/1‑110 (MQ=255)                                                                                                                                               gTCCAATGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGAGCAGCAGTTTCCCTATTTCGAGGTggagga                                                  <  2:608642/110‑1 (MQ=255)                                                                                                                                                tCCAATGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGAGCAGCAGTTTCCCTATTTCGAGGTGGAGGACTTCGAGCGCAAGCGACGCGAATATGGCCTTCCTGATGCAc         >  1:942391/1‑150 (MQ=255)                                                                                                                                                     tGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGAGCAGCAGTTTCCCTATTTCGAGGTGGAGGACTTCGAGCGCAAGCGACGCGAATATGGCCTTCCTGATGCACCGCAg    >  2:940331/1‑150 (MQ=255)                                                                                                                                                       tACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGAGCAGCAGTTTCCCTATTTCGAGGTGGAGGACTTCGAGCGCAAGCGACGCGAATATGGCCTTCCTGATGCACCGCAGCt  >  2:700303/1‑150 (MQ=255)                                                                                                                                                       tACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGAGCAGCAGTTTCCCTATTTCGAGGTGGAGGACTTCGAGCGCAAGCGACGCGAATATGGCCTTCCTGATGCACCGCAGCt  >  2:524120/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                      CAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGATTCCCTATCAGGTGGCCGGGCGTCATCAGGTACGCTTTTCCGTGGCGTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGAGCAGCAGTTTCCCTATTTCGAGGTGGAGGACTTCGAGCGCAAGCGACGCGAATATGGCCTTCCTGATGCACCGCAGCT  >  NC_002947/4314868‑4315166

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