Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	401,309	C→T	8.7%	Q298* (CAG→TAG)	PP_0334 →	transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	401,309	0	C	T	8.7%	168.5 / 7.7	69	Q298* (CAG→TAG)	PP_0334	transposase
Reads supporting (aligned to +/- strand):  ref base C (31/32);  new base T (3/3);  total (34/35)
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.63e-01

ACAGCAGCGGAAACAAAATCAGATTGGGCAAGATGAGCAAGCGCGGTGACGGCTACATACGCAGCCTGATGATCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGA  >  NC_002947/401190‑401457                                                                                                                        |                                                                                                                                                      gccccggtggcattacccggttggtcagaggctcatagccgcctccggtttttcgttcggtttgcgttggtaaacctgatcgctctgtaacagggcccatataatccgcagatttCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGc                                                                                                                        <  1:2441401‑M1/35‑1 (MQ=255)                         tcagaggctcatagccgcctccggtttttcgttcggtttgcgttggtaaacctgatcgctctgtaacagggcccatataatccgcagattcCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAaggagg                                                                                                <  1:2465046‑M1/59‑1 (MQ=255)                         tcagaggctcatagccgcctccggtttttcgttcggtttgcgttggtaaacctgatcgctctgtaacagggcccagataatccgcagattcCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAaggagg                                                                                                >  1:2613510‑M1/92‑150 (MQ=255)                             aggctcatagccgcctccggtttttcgttcggtttgcgttggtaaacctgatcgctctgtaacagggcccatataatccgcagattcCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAggcggc                                                                                            >  2:2465046‑M1/88‑150 (MQ=255)                                            tccggtttttcgttcggtttgcgttggtaaacctgatcgctctgtaacagggcccagataatccgcagattcCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAaggag                                                                                                 >  2:282920‑M1/73‑130 (MQ=255)                                            tccggtttttcgttcggtttgcgttggtaaacctgatcgctctgtaacagggcccagataatccgcagattcCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAaggag                                                                                                 <  1:282920‑M1/58‑1 (MQ=255)                                             ccggtttttcgttcggtttgcgttggtaaacctgatcgctctgtaacagggcccagataatccgcagatttCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAAt                                                                            >  1:894227‑M1/72‑150 (MQ=255)                                                   tttcgttcggtttgcgttggtaaacctgatcgctctgtaacagggcccatataatccgcagatttCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAAt                                                                      >  2:2441401‑M1/66‑150 (MQ=255)                                                                gCCTGATGATCCAAGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCCCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCt                                                         <  1:1121499/150‑1 (MQ=14)                                                                 ccTGATGATCCCGGGCGCGCCTTAAGTTCTGAGACATTTTTGCCCGGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTg                                                        <  2:2335884/150‑1 (MQ=11)                                                                  cTGATGATCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTgg                                                       <  1:151622/150‑1 (MQ=32)                                                                  cTGATGATCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTgg                                                       >  1:1822071/1‑150 (MQ=32)                                                                   ttggtaaacctgatcgctctgtaacagggcccagataatccgcagatttCGATTAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAggcg                                                                                              <  1:677978‑M1/61‑1 (MQ=255)                                                                   ttggtaaacctgatcgctctgtaacagggcccagataatccgcagatttCGATTAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAggcg                                                                                              >  2:677978‑M1/50‑110 (MQ=255)                                                                   tgatgaTCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTggg                                                      >  2:1964482/1‑150 (MQ=32)                                                                    gatgatCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGc                                                     >  2:2904624/1‑150 (MQ=32)                                                                     atgatCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGATCGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGcc                                                    <  1:1463490/150‑1 (MQ=17)                                                                       taaacctgatcgctctgtaacagggcccatataatccgcagatttCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAggcg                                                                                              >  2:770337‑M1/46‑106 (MQ=255)                                                                       taaacctgatcgctctgtaacagggcccatataatccgcagatttCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAggcg                                                                                              <  1:770337‑M1/61‑1 (MQ=255)                                                                         tCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGt                                                >  1:250916/1‑150 (MQ=32)                                                                         tCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGt                                                >  2:1096649/1‑150 (MQ=32)                                                                         tCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGt                                                >  2:1064077/1‑150 (MQ=32)                                                                         tCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGt                                                >  2:1479001/1‑150 (MQ=32)                                                                         tCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGt                                                <  1:351248/150‑1 (MQ=32)                                                                         tCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGt                                                >  2:2320780/1‑150 (MQ=32)                                                                         aacctgatcgctctgtgacagggcccagataatccgcagatttCGATTAGCCAGCCTGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAaggagg                                                                                                <  2:1731921‑M1/60‑1 (MQ=255)                                                                         aacctgatcgctctgtgacagggcccagataatccgcagatttCGATTAGCCAGCCTGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAaggagg                                                                                                >  1:1731921‑M1/44‑103 (MQ=255)                                                                           cAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTa                                              <  2:1027335/150‑1 (MQ=32)                                                                             gggCGCGCATTCCGTGTTGAGCCAGCTAGGGCCTGATTCCGATCAGGGAGG‑GGAGCGCCGCCTCTTGCGCTGGCTAGCCCGCCTTGGGCGCAAGGAGGCGGCGGTCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACa                                            <  2:1352522/150‑1 (MQ=11)                                                                             gggCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTCGCCTGCCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACa                                            <  1:1725100/150‑1 (MQ=255)                                                                                   gcATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGa                                      <  1:1064077/150‑1 (MQ=32)                                                                                     aTTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATc                                    >  1:1564217/1‑150 (MQ=32)                                                                                     aTTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATc                                    >  2:240981/1‑150 (MQ=32)                                                                                      ttCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCa                                   <  1:738898/150‑1 (MQ=32)                                                                                      ttCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCa                                   <  1:1964482/150‑1 (MQ=32)                                                                                        ccGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAgg                                 >  1:409072/1‑150 (MQ=32)                                                                                        ccGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCAAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAgg                                 <  1:941770/150‑1 (MQ=17)                                                                                         cGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGc                                        >  2:2804375/1‑142 (MQ=32)                                                                                         cGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGc                                        <  1:2804375/142‑1 (MQ=32)                                                                                         cGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGt                                >  1:2035775/1‑150 (MQ=32)                                                                                           tCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGttt                              <  2:1603054/150‑1 (MQ=32)                                                                                           tCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGttt                              <  2:1822071/150‑1 (MQ=32)                                                                                           cagggcccatataatccgcagatttCGATTAGCCAGCCTGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAggc                                                                                               <  2:500986‑M1/61‑1 (MQ=255)                                                                                           cagggcccatataatccgcagatttCGATTAGCCAGCCTGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAggc                                                                                               >  1:500986‑M1/26‑86 (MQ=255)                                                                                            cTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTa                             <  2:250916/150‑1 (MQ=32)                                                                                             ttGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTAc                            <  1:1479001/150‑1 (MQ=32)                                                                                             ttGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TCATCGCCGCCCCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTAc                            <  1:2465882/150‑1 (MQ=17)                                                                                               gAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGACCAGGTTCACCa                          >  2:794241/1‑150 (MQ=17)                                                                                                aGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCaa                         <  2:191150/150‑1 (MQ=32)                                                                                                 gCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAAc                        <  1:1061729/150‑1 (MQ=32)                                                                                                  ccAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACg                       >  1:1042523/1‑150 (MQ=32)                                                                                                        aaGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAAcc                 >  2:2792061/1‑150 (MQ=32)                                                                                                         aGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTAc                            <  1:1085631/138‑1 (MQ=33)                                                                                                         aGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTAc                            >  2:1085631/1‑138 (MQ=33)                                                                                                           gCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACcgaa              >  2:2369777/1‑150 (MQ=32)                                                                                                            ccTGATTCCGATCAGCCAGA‑TGATCTCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGATTTACCAACGCAAACcgaac             >  2:1498951/1‑150 (MQ=17)                                                                                                                aTTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACcgaac             <  1:363328/146‑1 (MQ=32)                                                                                                                aTTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACcgaac             >  2:363328/1‑146 (MQ=32)                                                                                                                aTTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGaaa         >  1:989547/1‑150 (MQ=32)                                                                                                                    CGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAAcc     <  1:2911287‑M1/150‑1 (MQ=255)                                                                                                                     GATCAGCCAGA‑TGATCGCCGGCTCTTGACCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCg    <  2:372052‑M1/150‑1 (MQ=255)                                                                                                                     GATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGACCGAAAAACCg    <  1:1667402‑M1/150‑1 (MQ=255)                                                                                                                     GATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCg    >  1:2828835‑M1/1‑150 (MQ=255)                                                                                                                      aTCAGCCAGA‑TGACCGCAGGAAATTGCACTGGCTAGCCCGCCTTGAAAGCAACGCGGCGAGTAACGGAAAGTCTGATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCgg   <  1:55965/150‑1 (MQ=2)                                                                                                                      ATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAAc                  >  2:69297‑M1/1‑135 (MQ=255)                                                                                                                      ATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAAc                  <  1:69297‑M1/135‑1 (MQ=255)                                                                                                                      ATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCgg   >  2:2730967‑M1/1‑150 (MQ=255)                                                                                                                       TCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGa  >  1:1075983‑M1/1‑150 (MQ=255)                                                                                                                       TCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGa  <  2:2066693‑M1/150‑1 (MQ=255)                                                                                                                        |                                                                                                                                                      ACAGCAGCGGAAACAAAATCAGATTGGGCAAGATGAGCAAGCGCGGTGACGGCTACATACGCAGCCTGATGATCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGA  >  NC_002947/401190‑401457

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