Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	755,850	C→A	6.5%	P148Q (CCG→CAG)	PP_0646 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	755,850	0	C	A	6.5%	190.6 / 2.2	61	P148Q (CCG→CAG)	PP_0646	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (25/32);  new base A (2/2);  total (27/34)
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

CGGTTGCCGTACGCACTAGCCTCACACCCGCCGGCCCTGGGAGCTCTAGCGCTTTGGCAGGGGCTTCGCCCGCCATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGATAACGGTTATGCCTGCACCCAGGAGG  >  NC_002947/755709‑755999                                                                                                                                              |                                                                                                                                                      cGGTTGCCGTACGCACTAGCCTCACACCCGCCGGCCCTGGGAGCTCTAGCGCTTTGGCAGGGGCTTCGCCCGCCATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTc                                                                                                                                               >  2:1368045/1‑150 (MQ=255)    ttGCCGTACGCACTAGCCTCACACCCGCCGGCCCTGGGAGCTCTAGCGCTTTGGCAGGGGCTTCGCCCGCCATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGt                                                                                                                                            >  2:1199837/1‑150 (MQ=255)       ccGAACGCACTAGCCTCACACCCGCCGGCCCTGGGAGCTCTAGCGCTTTGGCAGGGGCTTCGCCCGCCATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCTCAGATTCTGTTgg                                                                                                                                         >  1:761182/1‑150 (MQ=255)        cGTACGCACTAGCCTCACACCCGCCGGCCCTGGGAGCTCTAGCGCTTTGGCAGGGGCTTCGCCCGCCATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGa                                                                                                                                        <  1:1199837/150‑1 (MQ=255)        cGTACGCACTAGCCTCACACCCGCCGGCCCTGGGAGCTCTAGCGCTTTGGCAGGGGCTTCGCCCGCCATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGa                                                                                                                                        >  2:841602/1‑150 (MQ=255)          tACGCACTAGCCTCACACCCGCCGGCACTGGGAGCTCTAGCGCTTTGGCAGGGGCTTCGCCCGCCATGGTTGGCGGGTTGCCTTGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCTTTCCCGCAGATTCTGTTGTAGc                                                                                                                                      >  2:77927/1‑150 (MQ=255)             gCACTAGCCTCACACCCGCCGGCCCTGGGAGCTCTAGCGCTTTGGCAGGGGCTTCGCCCGGCATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTgg                                                                                                                                   <  2:1705639/150‑1 (MQ=255)                cTAGCCTCACACCCGCCGGCCCTGGGAGCTCTAGCGCTTTGGCAGGGGCTTCGCCCGCCATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTa                                                                                                                                >  2:1117646/1‑150 (MQ=255)                  aGCCTCACACCCGCCGGCCCTGGGAGCTCTAGCGCTTTGGCAGGGGCTTCGCCCGCCATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGc                                                                                                                              <  2:1882563/150‑1 (MQ=255)                   gCCTCACACCCGCCGGCCCTGGGAGCTCTAGCGCTTTGGCAGGGGCTTCGCCCGCCATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGcc                                                                                                                             >  2:1126954/1‑150 (MQ=255)                      tCACACCCGCCGGCCCTGGGAGCTCTAGCGCTTTGGCAGGGGCTTCGCCCGCCATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATg                                                                                                                          >  1:816320/1‑150 (MQ=255)                          acCCGCCGGCCCTGGGAGCTCTAGCGCTTTGGCAGGGGCTTCGCCCGCCATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCGCCGGAGCATGTCGTTCCCGCAGCTTCTGTTGGAGCTGGCTAGCCATGATTc                                                                                                                      >  2:799845/1‑150 (MQ=255)                          acCCGCCGGCCCTGGGAGCTCTAGCGCTTTGGCAGGGGCTTCGCCCGCCATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTCGAGCTGGCTAGCCATGATTc                                                                                                                      >  2:370719/1‑150 (MQ=255)                           cccGCCGGCCCTGGGAGCTCTAGCGCTTTGGCAGGGGCTTCGCCCGCCATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCg                                                                                                                     >  2:345801/1‑150 (MQ=255)                            ccgccgGCCCTGGGTGCTCTAGCGCTTTGGCAGGGGCTTCGCCCGCCATGGTTAGCAAATCCTCTCACGTGATTGTCGCGGGATTCTATGTGGTGCCCCGGCGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGa                                                                                                                    <  1:924769/150‑1 (MQ=255)                             cgccgGCCCTGGGAGTTCTAGCGCTTTGGCAGGGGCTTCGACCGCCATGGTTGGCGGGTTGCCTGGCCAGATGCAGCCCGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGAc                                                                                                                   <  2:168059/150‑1 (MQ=255)                              gccgGCCCTGGGAGCTCTAGCGCTTTGGCAGGGGCTTCGCCCGCCATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGAcc                                                                                                                  <  1:1681734/150‑1 (MQ=255)                                       gggAGCTCTAGCGCTTTGGCAGGGGCTTCGCCCGCCATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCAtt                                                                                                         <  2:726431/150‑1 (MQ=255)                                        ggAGCTCTAGCGCTTTGGCAGGGGCTTCGCCCGCCATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCAttt                                                                                                        <  1:841602/150‑1 (MQ=255)                                         gAGCTCTAGCGCTTTGGCAGGGGCTTCGCCCGCCATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTc                                                                                                       <  1:345801/150‑1 (MQ=255)                                            ctctAGCGCTTTGGCAGGGGCTTCGCCCGCCATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTa                                                                                                    <  1:1236741/150‑1 (MQ=255)                                               tAGCGCTTTGGCAGGGGCTTCGCCCGCCATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCAGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCtt                                                                                                                <  2:1239496/135‑1 (MQ=255)                                               tAGCGCTTTGGCAGGGGCTTCGCCCGCCATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCAGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCtt                                                                                                                >  1:1239496/1‑135 (MQ=255)                                                   gcTTTGGCAGGGGCTTCGCCCGCCATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCaa                                                                                             <  2:816320/150‑1 (MQ=255)                                                       tGGCAGGGGCTTCGCCCGCCATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGcc                                                                                         <  2:1065533/150‑1 (MQ=255)                                                         gCAGGGGCTTCGCCCGCCATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCtt                                                                                       <  1:1233138/150‑1 (MQ=255)                                                           aGGGGCTTCGCCCGCCATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGACTCGACCTTGCCCATTTCGTACCTGCAACGCCTTaa                                                                                     <  1:1117646/150‑1 (MQ=255)                                                                    gcccgccATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGTTTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACttt                                                                            <  1:210670/150‑1 (MQ=255)                                                                     cccgccATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCAGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTc                                                                           >  1:937710/1‑150 (MQ=255)                                                                      ccgccATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTAc                                                                                                   >  1:734667/1‑125 (MQ=255)                                                                      ccgccATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTAc                                                                                                   <  2:734667/125‑1 (MQ=255)                                                                          cATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCa                                                                      >  2:84121/1‑150 (MQ=255)                                                                             ccatGGCGGGCTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAggg                                                                   >  1:1920647/4‑150 (MQ=255)                                                                              gTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCAGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAgggg                                                                  <  2:937710/150‑1 (MQ=255)                                                                               ttGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGt                                                                 >  2:524685/1‑150 (MQ=255)                                                                                     ggTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAgg                                                           <  1:1029572/150‑1 (MQ=255)                                                                                         gCCTGGGCAGATGCAGCCGGGATTCTATGTGGCGCACCGTAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAgg                                                       <  1:593276/150‑1 (MQ=255)                                                                                          ccTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAggg                                                      >  1:1853561/1‑150 (MQ=255)                                                                                              ggCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCACCTTCTGTTGGACCAAAATCTGCATGATTCGACCGTGCCCATTTCGTACCTTCAACGCCTTAACCCGACGTTCGAGCAGGGGTTCAAGGCAGGGGaaa                                                  <  1:1304162/150‑1 (MQ=255)                                                                                               gCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCAGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAg                                                                     >  1:183890/1‑130 (MQ=255)                                                                                               gCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCAGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAg                                                                     <  2:183890/130‑1 (MQ=255)                                                                                                   aTGCAGCCGTGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTc                                             <  1:84121/150‑1 (MQ=255)                                                                                                   aTGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTc                                             <  1:1615885/150‑1 (MQ=255)                                                                                                   aTGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGCTCTCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACATACACCGCTATCTCCCGACTTTCGAGCAGGCGTTCAAGGCCGGGGAGATCTTc                                             <  2:761182/150‑1 (MQ=255)                                                                                                         ccGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATc                                       <  1:703896/150‑1 (MQ=255)                                                                                                          cGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCg                                      >  1:955305/1‑150 (MQ=255)                                                                                                           gggATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCgg                                     >  2:155074/1‑150 (MQ=255)                                                                                                                    tgtgGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGa                            <  2:955305/150‑1 (MQ=255)                                                                                                                     gtggtgCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCt                                                                                                 >  1:1677686/1‑80 (MQ=255)                                                                                                                     gtggtgACCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCt                                                                                                 <  2:1677686/80‑1 (MQ=255)                                                                                                                         tgCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGATCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGCCCTCTGCCCACTGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGATAACg                       >  1:1574663/1‑150 (MQ=255)                                                                                                                            cccGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGATAACGGtt                    >  2:608537/1‑150 (MQ=255)                                                                                                                              cGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAggg                                                      >  1:1436637/1‑114 (MQ=255)                                                                                                                              cGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAggg                                                      <  2:1436637/114‑1 (MQ=255)                                                                                                                              cGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGACCAGGTGTTCAAGTCAGTGGAAATCTTCGTCATCGGCGATCCGGATAACGGTTAt                  <  1:530943/150‑1 (MQ=255)                                                                                                                                gAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGATAACGGTTATGc                <  1:524685/150‑1 (MQ=255)                                                                                                                                 aGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGATAACGGTTATGcc               <  1:1126954/150‑1 (MQ=255)                                                                                                                                         gTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGATAACGGTTATGCCTGCACCCa       <  2:1208528/150‑1 (MQ=255)                                                                                                                                          ttCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGATAACGGTTATGc                >  2:1838243/1‑140 (MQ=255)                                                                                                                                          ttCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGATAACGGTTATGc                <  1:1838243/140‑1 (MQ=255)                                                                                                                                          ttCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGATAACGGTTATGCCTGCACCCag      >  2:608232/1‑150 (MQ=255)                                                                                                                                            cccGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGATAACGGTTATGCCTGCACCCagga    <  2:1853561/150‑1 (MQ=255)                                                                                                                                             ccGCAGATTCTGTTGGATCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGATAACGGTTATGCCTGCACCCaggag   >  1:383380/1‑150 (MQ=255)                                                                                                                                             ccGCAGATTCTGTTGGATCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGATAACGGTTATGCCTGCACCCaggag   >  1:936755/1‑150 (MQ=255)                                                                                                                                              cGCAGCTTCTGTTGGAGCTGGCTAGCCATGCTGCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACGTTAGAGCAGGGCTTTAAGGCAGGGGAAATCTTCGTCATCGGCGAACCGCATAACGGTTATGCCGGCACCCcggcgg  >  1:41325/1‑150 (MQ=255)                                                                                                                                              cGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGATAACGGTTATGCCTGCACCCaggagg  <  1:512490/150‑1 (MQ=255)                                                                                                                                              cGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGACGTTGGCGATCCGGATAACGGTTATGCCTGCACCCaggagg  >  1:509937/1‑150 (MQ=255)                                                                                                                                              |                                                                                                                                                      CGGTTGCCGTACGCACTAGCCTCACACCCGCCGGCCCTGGGAGCTCTAGCGCTTTGGCAGGGGCTTCGCCCGCCATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGATAACGGTTATGCCTGCACCCAGGAGG  >  NC_002947/755709‑755999

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