Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,000,920	G→A	6.1%	Q26* (CAG→TAG)	PP_2624 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,000,920	0	G	A	6.1%	195.1 / 2.4	65	Q26* (CAG→TAG)	PP_2624	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (31/30);  new base A (2/2);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

GTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCAACGGATCCCGTATTTTTAGACATCGGTACGCAACTCCCTGTTGGGCTGTAGGATGTTCAGGTGCTAGTTGTGTTTATTATTGTCTCGGTCTT  >  NC_002947/3000774‑3001064                                                                                                                                                   |                                                                                                                                                  gTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGc                                                                                                                                                  >  8:110722/1‑148 (MQ=255)  tCCAGGCTCACATTGGCGTCGTTCAAGACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGACGCTCCCCCTCAGCCTgat                                                                                                                                                 >  6:201957/1‑146 (MQ=255)  tCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCCGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCt                                                                                                                                                 >  6:206658/1‑148 (MQ=255)   ccAGGCTCACATTGGCGTCGATCAACACGCTGTTGAATGTATGTTTATAGATGCGCATGACCTGACGGTCCTCAAGCGCACGTCAGTCATCCAGGCCGAAGTCGCCCAGCACCAGCAGCTTGTGCGGCAGCTACACCTCAGCCTGCt                                                                                                                                                 >  6:289617/1‑147 (MQ=21)         tCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTg                                                                                                                                         >  1:313223/1‑149 (MQ=255)                   gtcgtTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGa                                                                                                                               >  3:109635/1‑149 (MQ=255)                       ttCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGa                                                                                                                               <  4:307808/145‑1 (MQ=255)                       ttCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGa                                                                                                                               >  3:307808/1‑145 (MQ=255)                       ttCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATa                                                                                                                           >  7:18616/1‑149 (MQ=255)                       ttCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGAAGCGGGGACATa                                                                                                                           <  6:262001/149‑1 (MQ=255)                         cAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATAc                                                                                                                          <  5:272688/148‑1 (MQ=255)                           acacGCTGTTGAAGGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTg                                                                                                                       <  3:266637/149‑1 (MQ=255)                            cacGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGa                                                                                                                      <  3:206920/149‑1 (MQ=255)                                      aaTGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑Acgc                                                                                                            <  6:144429/148‑1 (MQ=255)                                                 aTCCATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGgc                                                                                                 <  7:82736/148‑1 (MQ=255)                                                  tCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGgcg                                                                                                >  3:299067/1‑148 (MQ=255)                                                    gATGCGCATGACCTGGCGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCaa                                                                                             >  6:221579/1‑149 (MQ=255)                                                    gATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCaa                                                                                             >  7:189926/1‑149 (MQ=255)                                                     aTGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCaag                                                                                            >  5:263884/1‑148 (MQ=255)                                                        cgcATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCAACgg                                                                                          <  4:109635/148‑1 (MQ=255)                                                          cATGACCTGGCGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCAACGGATc                                                                                       <  5:221579/149‑1 (MQ=255)                                                           aTGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCAACGGATcc                                                                                      >  6:197455/1‑149 (MQ=255)                                                            tGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCAACGGATccc                                                                                     >  5:142353/1‑149 (MQ=255)                                                                    cGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCAACGGATCCCGTAttttt                                                                             <  4:170948/149‑1 (MQ=255)                                                                    cGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCAACGGATCCCGTAttttt                                                                             <  6:59118/149‑1 (MQ=255)                                                                       tCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTACTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCAACGGATCCCGTATTTTTAGa                                                                          >  6:133440/1‑149 (MQ=255)                                                                       tCCTCAACCGCACGGCTGTCATCCACGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGACTACTCATCCCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCAAAGGATCCCCTATTTTTAGa                                                                          >  6:133444/1‑149 (MQ=255)                                                                         cTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCAACGGATCCCGTATTTTTAGACa                                                                        >  6:208358/1‑149 (MQ=255)                                                                          tCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCAACGGATCCCGTATTTTTAGACAt                                                                       >  5:323585/1‑149 (MQ=255)                                                                          tCAAGCGCACGGCTGTCATCCAGGCCGAAGTCACCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCAACGGATCCCGTAtt                                                                                <  1:234961/140‑1 (MQ=255)                                                                          tCAAGCGCACGGCTGTCATCCAGGCCGAAGTCACCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCAACGGATCCCGTAtt                                                                                >  2:234961/1‑140 (MQ=255)                                                                            aaGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑Ac                                                                                                              >  6:166974/1‑108 (MQ=255)                                                                            aaGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑Ac                                                                                                              <  5:166974/108‑1 (MQ=255)                                                                                gcACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCAACGGATCCCGTAttttt                                                                             >  3:194676/1‑137 (MQ=255)                                                                                gcACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCAACGGATCCCGTAttttt                                                                             <  2:194263/137‑1 (MQ=255)                                                                                gcACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCAACGGATCCCGTAttttt                                                                             >  1:194263/1‑137 (MQ=255)                                                                                gcACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCAACGGATCCCGTAttttt                                                                             <  4:194676/137‑1 (MQ=255)                                                                                  aCGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCAACGGATCCCGTATTTTTAGACATCGGTACGc                                                               <  8:18616/149‑1 (MQ=255)                                                                                    ggCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCAACGGATCCCGTATTTTTAGACATCGGTACGCaa                                                             <  4:117601/149‑1 (MQ=255)                                                                                           aTCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCAACGGATCCCGTATTTTTAGACATCGGTACGCAACTCCCTg                                                      >  7:30157/1‑149 (MQ=255)                                                                                           aTCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCAACGGATCCCGTATTTTTAGACATCGGTACGCAACTCCCTg                                                      >  7:196530/1‑149 (MQ=255)                                                                                           aTCCAGGACGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCAACGGATCCCGTATTTTTAGACATCGGTACGCAACTCCCTg                                                      >  7:213098/1‑149 (MQ=255)                                                                                              cAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCAACGGATCCCGTATTTTTAGACATCGGTACGCAACTCCCTGtt                                                    <  6:142353/148‑1 (MQ=255)                                                                                                ggCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCAACGGATCCCGTATTTTTAGACATCGGTACGCAACTCCCTGTTggg                                                 <  4:299067/149‑1 (MQ=255)                                                                                                ggCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCAACGGATCCCGTATTTTTAGACATCGGTACGCAACTCCCTGTTggg                                                 <  7:110722/149‑1 (MQ=255)                                                                                                    gAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCAACGGATCCCGTATTTTTAGACATCGGTACGCAACTCCCTGTTGGGCTg                                              >  2:25607/1‑148 (MQ=255)                                                                                                     aaGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCAACGGATCCCGTATTTTTAGACATCGGTACGCAACTCCCTGTTGGGCTGTa                                            <  1:25607/149‑1 (MQ=255)                                                                                                          gCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCAACGGATCCCGTATTTTTAGACATCGGTACGCAACTCCCTGTTGGGCTGTAGGATg                                       <  8:30157/149‑1 (MQ=255)                                                                                                               gcacctgcaTCTTGTGCGGCAGCTCCACCTCAGCCTACTCATCGCTGGTAGCGGGGTCATACTTGTTTTTG‑ACGCGTTCTTTCGGCGCAAGGGATCCCGTATTTTTATACATCGGTACGCAACTCCCTGTTGGGCTGTAGGATGTTCat                                  <  5:133444/149‑2 (MQ=255)                                                                                                               gcaccagcaTCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCAACGGATCCCGTATTTTTAGACATCGGTACGCAACTcc                                                         <  5:6797/126‑1 (MQ=255)                                                                                                               gcaccagcaTCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCAACGGATCCCGTATTTTTAGACATCGGTACGCAACTcc                                                         >  6:6797/1‑126 (MQ=255)                                                                                                               gcaccagcaTCTTGTGCGGCAGCTCCACCTCAGCCTACTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCAACGGATCCCGTATTTTTAGACATCGGTACGCAACTCCCTGTTGGGCTGTAGGATGTTCAg                                  <  5:133440/149‑1 (MQ=255)                                                                                                                   cagcaTCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCAACGGATCCCGTATTTTTAGACATCGGTACGCAACTCCCTGTTGGGCTGTAGGATg                                       <  8:12106/140‑1 (MQ=255)                                                                                                                   cagcaTCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCAACGGATCCCGTATTTTTAGACATCGGTACGCAACTCCCTGTTGGGCTGTAGGATg                                       >  7:12106/1‑140 (MQ=255)                                                                                                                   cagcaTCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCAACGGATCCCGTATTTTTAGACATCGGTACGCAACTCCCTGTTGGGCTGTAGGATGTTCAGGTGc                              >  4:122253/1‑149 (MQ=255)                                                                                                                    agcaTCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCAACGGATCCCGTATTTTTAGACATCGGTACGCAACTCCCTGtt                                                    <  3:98417/126‑1 (MQ=255)                                                                                                                    agcaTCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCAACGGATCCCGTATTTTTAGACATCGGTACGCAACTCCCTGtt                                                    >  4:98417/1‑126 (MQ=255)                                                                                                                        tCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTGAACGCGTTCTTTCGGCGCAACGGATCCCGTATTTTTAGACATCGGTACGCAACTCCCTGTTGGGCTGTAGGATTTTCAGGTGCTAGt                          >  1:94918/1‑149 (MQ=255)                                                                                                                          ttGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCAACGGATCCCGTATTTTTAGACATCGGTACGCAACTCCCTGTTGGGCTGTAGGATGTTCAGGTGCTAGTtgt                       >  5:105913/1‑149 (MQ=255)                                                                                                                            gtgCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCAACGGATCCCGTATTTTTAGACATCGGTACGCAACTCCCTGTTGGGCTGTAGGATGTTCAGGTGCTAGtt                         <  5:208358/145‑1 (MQ=255)                                                                                                                             tgCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCAACGGATCCCGTATTTTTAGACATCGGTACGCAACTCCCTGTTGGGCTGTAGGATGTTCAGGTGCTAGTTgtgt                     <  6:105913/148‑1 (MQ=255)                                                                                                                                  cAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACATGATATTG‑ACGCGTTCTTTCGGCGCAACGGATCCCGTATTTTTAGACATCGGTACGCAACTCCCTGTTGGGCTGTAGGATGTTCAGGTGCTAGTTGTGTttatta               >  2:72526/1‑149 (MQ=255)                                                                                                                                    gCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCAACGGATCCCGTATTTTTAGACATCGGTACGCAACTCCCTGTTGGGCTGTAGGATGTTCAGGTGCTAGTTGTGTTTattatt             <  5:197455/149‑1 (MQ=255)                                                                                                                                         aCCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCAACGGATCCCGTATTTTTAGACATCGGTACGCAACTCCCTGTTGGGCTGTAGGATGTTCAGGTGCTAGTTGTGTTTATTATTGtctc        >  6:39322/1‑149 (MQ=255)                                                                                                                                              aGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCAACGGATCCCGTATTTTTAGACATCGGTACGCAACTCCCTGTTGGGCTGTAGGATGTTCAGGTGCTAGTTGTGTTTATTATTGTCTCGGTc    <  2:94918/148‑1 (MQ=255)                                                                                                                                               gCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCAACGGATCCCGTATTTTTAGACATCGGTACGCAACTCCCTGTTGGGCTGTAGGATGTTCAGGTGCTAGTTGTGTTTATTATTGTCTCGGTCtt  <  2:313223/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                  GTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCCTGCTCATCGCTGGTAGCGGGGACATACTTGATATTG‑ACGCGTTCTTTCGGCGCAACGGATCCCGTATTTTTAGACATCGGTACGCAACTCCCTGTTGGGCTGTAGGATGTTCAGGTGCTAGTTGTGTTTATTATTGTCTCGGTCTT  >  NC_002947/3000774‑3001064

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