Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,189,684	A→C	6.6%	N170H (AAC→CAC)	PP_1937 →	helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,189,684	0	A	C	6.6%	169.3 / 4.5	61	N170H (AAC→CAC)	PP_1937	helicase
Reads supporting (aligned to +/- strand):  ref base A (29/28);  new base C (2/2);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.28e-01

ATCGAGCTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAG  >  NC_002947/2189537‑2189827                                                                                                                                                    |                                                                                                                                                aTCGAGCTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAAc                                                                                                                                               >  2:328123/1‑150 (MQ=255)     aGCTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTgc                                                                                                                                           >  2:884442/1‑150 (MQ=255)         gATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCtt                                                                                                                                       <  2:329092/150‑1 (MQ=255)                gcCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGaa                                                                                                                                                >  1:730859/1‑134 (MQ=255)                gcCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGaa                                                                                                                                                <  2:730859/134‑1 (MQ=255)                 ccGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTg                                                                                                                               <  2:1272113/150‑1 (MQ=255)                     tCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGa                                                                                                                           >  1:439032/1‑150 (MQ=255)                         gTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATATTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAaga                                                                                                                       <  1:1356128/150‑1 (MQ=255)                         gTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAaga                                                                                                                       >  2:1356128/1‑150 (MQ=255)                         gTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACCGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAaga                                                                                                                       <  1:489343/150‑1 (MQ=255)                         gTGCGAATTCGTTACCCGACAGGGTAAAGACGGCTACATGCAAAAATGTACTGCTCCGATGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACATACGCTTCACCTCTCGGGAAaga                                                                                                                       >  2:297556/1‑150 (MQ=255)                          tGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGCACGTTCGCTTCatgctga                                                                                                                               >  2:635310/1‑135 (MQ=255)                          tGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGCACGTTCGCTTCatgctga                                                                                                                               <  1:635310/141‑7 (MQ=255)                                   gTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTttg                                                                                                             <  1:728997/150‑1 (MQ=255)                                    ttACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACATGCGCTTCACCTCTGGGGAAAGAGCGCGCTttgt                                                                                                            <  1:346474/150‑1 (MQ=255)                                     tACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCACATAAACTAGAACAGGTTTACCAACTCAGCATAAAGCGAACGTAAGCTTCACCTGTGGGGAAAGCGCGCGCTttgcc                                                                                                           <  2:287789/150‑3 (MQ=255)                                     ctcccGACAGGGTACAGACGGCTACATGCAAAAATGTACCGCTCCGCTGGAGCAGTGACGCTTGGCTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGAGCATCCCCTCTGGGGAAAGAGCGCGCTttgtt                                                                                                           <  1:397765/148‑1 (MQ=255)                                        ccGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATCCTCAGGCAGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTgag                                                                                                        >  2:397765/1‑150 (MQ=255)                                        ccGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGTGCGCTCTGTTgag                                                                                                        >  2:807805/1‑150 (MQ=255)                                         cGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGCACGTTCGCTTCa                                                                                                                                     <  1:146447/120‑1 (MQ=255)                                         cGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGCACGTTCGCTTCa                                                                                                                                     >  2:146447/1‑120 (MQ=255)                                          gACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGGAAGAGCGCGCTTTGTTGAGAc                                                                                                      >  1:966467/1‑150 (MQ=255)                                               ggTACAGACGGCTACATGCACAACTGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGGAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAAtc                                                                                                 >  2:381921/1‑150 (MQ=255)                                                         gCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTGCTTACTCAGCATGAAGCTAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGa                                                                                       <  1:884442/150‑1 (MQ=255)                                                                 cAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCtttt                                                                               >  2:743631/1‑150 (MQ=255)                                                                 cAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCtttt                                                                               <  1:743631/150‑1 (MQ=255)                                                                  aaaaaTGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTg                                                                              >  1:419714/1‑150 (MQ=255)                                                                  aaaaaTGTACTGCTCCGCTGGAGCAGTGACGCTAGGGCCCCAAGGCTCAACGCATCGGGGGTTAGATACTCATCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGGGGTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTg                                                                              >  1:387698/1‑150 (MQ=255)                                                                    aaaTGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGa                                                                            <  1:644017/150‑1 (MQ=255)                                                                       tGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTTGACCAGGTTTACATACTCAGCATGAATCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCt                                                                         <  1:297556/150‑1 (MQ=255)                                                                        tcaCTGCTCCAGCGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCtt                                                                        >  2:1204210/3‑150 (MQ=255)                                                                          aCTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTgg                                                                      >  2:1098564/1‑150 (MQ=255)                                                                              cTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCttt                                                                                                              <  1:702380/106‑1 (MQ=255)                                                                              cTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCttt                                                                                                              >  2:702380/1‑106 (MQ=255)                                                                               tCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAAtctc                                                                                               >  1:1330392/1‑120 (MQ=255)                                                                               tCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAAtctc                                                                                               <  2:1330392/120‑1 (MQ=255)                                                                               tCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTcc                                                                                              <  1:95143/121‑1 (MQ=255)                                                                               tCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTcc                                                                                              >  2:1244443/1‑121 (MQ=255)                                                                               tCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTcc                                                                                              <  1:1244443/121‑1 (MQ=255)                                                                               tCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTcc                                                                                              >  2:95143/1‑121 (MQ=255)                                                                               tCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGgaga                                                                 >  1:1281310/1‑150 (MQ=255)                                                                                    cggAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGc                                                            <  1:1204210/149‑1 (MQ=255)                                                                                        gCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGaa                                                        >  2:1343104/1‑150 (MQ=255)                                                                                        gCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGaa                                                        <  2:439032/150‑1 (MQ=255)                                                                                            tGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACg                                                    <  1:1098564/150‑1 (MQ=255)                                                                                                gCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCttt                                                >  2:520167/1‑150 (MQ=255)                                                                                                       tCCAGCAACTACACTGGACCAGGTTTACTTTCTCAGCATGAAGCGAACGTGCGCTTCATCTCTGGGGAAAGAGCGTGCTTTGTTGAGACCAATCTCCTAGCTTAGTGCTTTTTTATCTGTGGGATAAACGCCGAAGACGCTTTCATTTAg                                         <  2:648853/150‑1 (MQ=255)                                                                                                           gCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAAc                                                              <  2:704032/125‑1 (MQ=255)                                                                                                           gCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAAc                                                              >  1:704032/1‑125 (MQ=255)                                                                                                           gCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAgctgc                                     <  2:1281310/150‑1 (MQ=255)                                                                                                                   acTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTAGagcag                             >  2:1044849/1‑150 (MQ=255)                                                                                                                    cTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTAGagcagc                            <  1:520167/150‑1 (MQ=255)                                                                                                                             gTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTAGAGCAGCTtctc                       >  2:620794/1‑146 (MQ=255)                                                                                                                             gTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTAGAGCAGCTtctc                       <  1:620794/146‑1 (MQ=255)                                                                                                                             gTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGt                   >  1:277015/1‑150 (MQ=255)                                                                                                                                   aTACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCATTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCa             >  2:908883/1‑150 (MQ=255)                                                                                                                                     aCTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCg                                                          <  1:952492/103‑1 (MQ=255)                                                                                                                                     aCTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCg                                                          >  2:952492/1‑103 (MQ=255)                                                                                                                                      cTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCt          <  1:1044849/150‑1 (MQ=255)                                                                                                                                       tCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTg         >  2:1153983/1‑150 (MQ=255)                                                                                                                                            aTGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCa             >  2:1239061/1‑141 (MQ=255)                                                                                                                                            aTGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCa             <  1:1239061/141‑1 (MQ=255)                                                                                                                                              gAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAagag  <  1:1343104/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                ATCGAGCTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAG  >  NC_002947/2189537‑2189827

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