Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,315,014	A→T	9.8%	N391I (AAT→ATT)	aspC →	aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,315,014	0	A	T	9.8%	177.4 / 10.0	61	N391I (AAT→ATT)	aspC	aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (27/28);  new base T (3/3);  total (30/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 = 2.75e-01

CCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGATTCCCTATCAGGTGGCCGGGCGTCATCAGGTACGCTTTTCCGTGGCGTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGAGCAGCAGTTTCCCTATTTCGAGGTGGAGGACTTCGAGCGCAAGCGACGCGAATATGGCCTTCCTGATGCAC  >  NC_002947/4314867‑4315159                                                                                                                                                     |                                                                                                                                                  ccAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGATTCCCTATCAGGTGGCCGGGCGTCATCAGGTACGCTTTTCCGTGGCGTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATg                                                                                                                                                 >  1:651314/1‑150 (MQ=255)  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:65045/150‑1 (MQ=255)                       ggCCACTGAGCTGGGTATCCTGACGATTCCCTATCAGGTGGCCGGGCGTCATCAGGTACGCTTTTCCGTGGCGTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCGagag                                                                                                                           <  2:831785/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                                                                                                                                >  2:802373/1‑125 (MQ=255)                                           tGACGATTCCCTATCAGGTGGCCGGGCGTCATCAGGTACGCTTTTCCGTGGCGTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGcc                                                                                                                                <  1:802373/125‑1 (MQ=255)                                           tGACGATTCCCTATCAGGTGGCCGGGCGTCATCAGGTACGCTTTTCCGTGGCGTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGtctc                                                                                                       >  2:601743/1‑150 (MQ=255)                                            gACGATTCCCTATCAGGTGGCCGGGCGTCATCAGGTACGCTTTTCCGTGGCGTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGtctct                                                                                                      >  2:1300126/1‑150 (MQ=255)                                            gACGATTCCCTATCAGGTGGCCGGGCGTCATCAGGTACGCTTTTCCGTGGCGTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGtctct                                                                                                      >  1:1126221/1‑150 (MQ=255)                                                     cTATCAGGTGGCCGGGCGTCATCAGGTACGCTTTTCCGTGGCGTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCGagag                                                                                                                           >  2:1339445/1‑120 (MQ=255)                                                     cTATCAGGTGGCCGGGCGTCATCAGGTACGCTTTTCCGTGGCGTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCGagag                                                                                                                           <  1:1339445/120‑1 (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                                                                              <  2:601416/148‑1 (MQ=255)                                                                      gTCATCAGGTACGCTTTTCCGTGGCGTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTc                                                                              >  1:601416/1‑148 (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                                                                                                                                   <  2:890626/57‑1 (MQ=255)                                                                                                            gCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCtt                                                                                                                                   >  1:890626/1‑57 (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                                              <  1:238887/111‑1 (MQ=255)                                                                                                                                          tgaaCTGTACATTGGACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGAGCAGCAGTTTCCCTATTTCGAGGTgga                                              >  2:238887/3‑113 (MQ=255)                                                                                                                                          tgaaCTGTACATTGGACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGAGCAGCAGTTTCCCTATTTCGAGGTGGAGGACTTCGAGCGCAAGCGACGCGAATATg                 >  2:150208/3‑142 (MQ=255)                                                                                                                                          tgaaCTGTACATTGGACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGAGCAGCAGTTTCCCTATTTCGAGGTGGAGGACTTCGAGCGCAAGCGACGCGAATATg                 <  1:150208/140‑1 (MQ=255)                                                                                                                                          tgaaCTGTACATTGGACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGAGCAGCAGTTTCCCTATTTCGAGGTGGAGGACTTCGAGCGCAAGCGACGCGAATATGGCCTTcc          >  2:1189221/3‑149 (MQ=255)                                                                                                                                          tgaaCTGTACATTGGACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGAGCAGCAGTTTCCCTATTTCGAGGTGGAGGACTTCGAGCGCAAGCGACGCGAATATGGCCTTcc          <  1:1189221/147‑1 (MQ=255)                                                                                                                                            aaCTGTCCAATGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGAGCAGCAGTTTCCCTATTTCGAGGTGGAGGACTTCGAGCGCAAGCGACGCGAATATGGCCTTCCTGa       <  2:464129/150‑1 (MQ=255)                                                                                                                                             aCTGTCCAATGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGAGCAGCAGTTTCCCTATTTCGAGGTGGAGGACTTCGAGCGCAAGCGACGCGAATATGGCCTTCCTGAt      <  1:290446/150‑1 (MQ=255)                                                                                                                                                gTCCAATGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGAGCAGCAGTTTCCCTATTTCGAGGTggagga                                           <  2:608642/110‑1 (MQ=255)                                                                                                                                                gTCCAATGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGAGCAGCAGTTTCCCTATTTCGAGGTggagga                                           >  1:608642/1‑110 (MQ=255)                                                                                                                                                 tCCAATGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGAGCAGCAGTTTCCCTATTTCGAGGTGGAGGACTTCGAGCGCAAGCGACGCGAATATGGCCTTCCTGATGCAc  >  1:942391/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                  CCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGATTCCCTATCAGGTGGCCGGGCGTCATCAGGTACGCTTTTCCGTGGCGTTCTGCGGCGAAGCCGAAACCGTGTTCAACCGTTTGCGT‑CAGCAACTGTCCAATGTACAGTTCAGCCTTGCCGAGAGGGTATCGGCATGACTGTCTCTCTTCGCAATGTCGACAGTGTCCTCGGCGAGCAGCAGTTTCCCTATTTCGAGGTGGAGGACTTCGAGCGCAAGCGACGCGAATATGGCCTTCCTGATGCAC  >  NC_002947/4314867‑4315159

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