Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	616,401	G→A	6.0%	noncoding (121/142 nt)	PP_mr11 ←	FMN regulatory region

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	616,401	0	G	A	6.0%	200.0 / 2.0	66	noncoding (121/142 nt)	PP_mr11	FMN regulatory region
Reads supporting (aligned to +/- strand):  ref base G (32/30);  new base A (2/2);  total (34/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 = 7.89e-01

GCGGCGGAAACATTGGGGAATTGCGCGTGGTGCTGGGTGAACATGAAACGCTCCTCGTGAATGATCGATGAACGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGT  >  NC_002947/616258‑616540                                                                                                                                                 |                                                                                                                                            gcggcgGAAACATTGGGGAATTGCGCGTGGTGCTGGGTGAACATGAAACGCTCCTCGTGAATGATCGATGAACGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTAt                                                                                                                                        <  3:304303/149‑1 (MQ=255)   ggcggAAACATTGGGGAATTGCGCGTGGTGCTGGGTGAACATGAAACGCTCCTCGTGAATGATCGATGAACGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGa                                                                                                                                      <  6:167642/149‑1 (MQ=255)      ggAAACATTGGGGAATTGCGCGTGGTGCTGGGTGAACATGAAACGCTCCTCGTGAATGATCGATGAACGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGa                                                                                                                                      <  1:177165/146‑1 (MQ=255)      ggAAACATTGGGGAATTGCGCGTGGTGCTGGGTGAACATGAAACGCTCCTCGTGAATGATCGATGAACGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGa                                                                                                                                      >  2:177165/1‑146 (MQ=255)      ggAAACATTGGGGAATTGCGCGTGGTGCTGGGTGAACATGAAACGCTCCTCGTGAATGATCGATGAACGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCg                                                                                                                                   >  8:194608/1‑149 (MQ=255)        aaaCATTGGGGAATTGCGCGTGGTGCTGGGTGAACATGAAACGCTCCTCGTGAATGATCGATGAACGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTc                                                                                                                                 >  4:121817/1‑149 (MQ=255)        aaaCATTGGGGAATTGCGCGTGGTGCTGGGTGAACATGAAACGCTCCTCGTGAATGATCGATGAACGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTc                                                                                                                                 >  8:201225/1‑149 (MQ=255)         aaCATTGGGGAATTGCGCGTGGTGCTGGGTGAACATGAAACGCTCCTCGTGAATGATCGATGAACGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTCg                                                                                                                                >  2:263453/1‑149 (MQ=255)             ttGGGGAATTGCGCGTGGTGCTGGGTGAACATGAAACGCTCCTCGTGAATGATCGATGAACGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTCGGctc                                                                                                                            >  2:27467/1‑149 (MQ=255)                           gtggtgCTGGGTGAACATGAAACGCTCCTCGTGAATGATCGATGAACGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCAc                                                                                                                 >  6:197691/1‑146 (MQ=255)                                   gggTGAACATGAAACGCTCCTCGTGAATGATCGATGAACGTCTCGGGGCTAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGACCGCCTTCTTTCATCCGGACTATGACCGTCGGctct                                                                                                                           <  1:202469/129‑1 (MQ=255)                                   gggTGAACATGAAACGCTCCTCGTGAATGATCGATGAACGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTCGGctct                                                                                                                           >  2:202469/1‑128 (MQ=255)                                    ggTGAACATGAAACGCTCCTCGTGAATGATCGATGAACGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATctgct                                                                                                       <  7:201225/147‑1 (MQ=255)                                     gTGAACATGAAACGCTCCTCGTGAATGATCGATGAACGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGAc                                                                                                    <  1:27467/149‑1 (MQ=255)                                            tGAAACGCTCCTCGTGAATGATCGATGAACGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTCGGCTCAGGAGTTGCACCAGATCTGCTGACCCCCGAc                                                                                             <  1:25788/149‑1 (MQ=255)                                             gAAACGCTCCTCGTGAATGATCGATGAACGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGAACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACa                                                                                            >  7:112126/1‑149 (MQ=255)                                             gAAACGCTCCTCGTGAATGATCGATGAACGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGAATGA‑CGCCTTCTTTCATCCGAACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACa                                                                                            >  3:197982/1‑149 (MQ=255)                                              aaaCGCTCCTCGTGAATGATCGATGAACGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTCCACCAGATCTGCTGACCCCCGACa                                                                                            <  7:194608/148‑1 (MQ=255)                                                            aaTGATCGATGAACGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTcg                                                                             >  8:288425/1‑149 (MQ=255)                                                                aTCGATGAACGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCggg                                                                         >  2:28822/1‑149 (MQ=255)                                                                 tCGATGAACGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCggg                                                                         <  6:156880/148‑1 (MQ=255)                                                                 tCGATGAACGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGc                                                                        >  7:281597/1‑149 (MQ=255)                                                                  cGATGAACGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGc                                                                        >  4:64937/1‑148 (MQ=255)                                                                  cGATGAACGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGc                                                                        <  3:64937/148‑1 (MQ=255)                                                                   gATGAACGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCt                                                                       <  7:28792/148‑1 (MQ=255)                                                                          gTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTc                                                               >  5:298176/1‑149 (MQ=255)                                                                               ggggCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGAtt                                                           <  5:104059/148‑1 (MQ=255)                                                                                gggCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGCCATGGTCGGGGCGCTCGCGGGCTCCTCTTTAGAttt                                                          >  4:4628/1‑148 (MQ=255)                                                                                gggCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTa                                                         <  7:288425/149‑1 (MQ=255)                                                                                   cGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTAccg                                                      <  5:197691/149‑1 (MQ=255)                                                                                    gAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCggg                                                                         <  1:253927/129‑1 (MQ=255)                                                                                    gAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCggg                                                                         >  2:253927/1‑129 (MQ=255)                                                                                    gAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTAccgc                                                     >  2:191180/1‑149 (MQ=255)                                                                                     aaCAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTAccgcc                                                    >  3:156140/1‑149 (MQ=255)                                                                                       cAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCgg                                                  >  1:309057/1‑149 (MQ=255)                                                                                        aaaaTGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGt                                                 <  1:191180/149‑1 (MQ=255)                                                                                          aaTGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTgg                                               >  1:136247/1‑149 (MQ=255)                                                                                           atccgcgcgcAAGGGCGCCCGGATGGCCCCTGCATCA‑CGCCTTCTTTCATCCGGACTAAGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTggg                                              <  8:27947/146‑1 (MQ=255)                                                                                           aTGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATCACCGTCGGCGCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTggg                                              <  8:32830/149‑1 (MQ=255)                                                                                              cgcgcgcAAGGGCGCCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGa                                            <  8:281597/148‑1 (MQ=255)                                                                                              cgcgcgcAAGGGCGCCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGa                                            <  6:298176/148‑1 (MQ=255)                                                                                               gcgcgcAAGGGCGCCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACt                                          >  7:176601/1‑149 (MQ=255)                                                                                               gcgcgcAAGGGCGCCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACt                                          >  8:51961/1‑149 (MQ=255)                                                                                                        ggCGCCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGAc                                           >  2:298801/1‑139 (MQ=255)                                                                                                        ggCGCCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGAc                                           <  1:298801/139‑1 (MQ=255)                                                                                                         gcgcCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGc                                >  6:281523/1‑149 (MQ=255)                                                                                                          cgcCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGcc                               <  4:156140/149‑1 (MQ=255)                                                                                                          cgcCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTCGGCTATGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGcc                               <  4:47514/149‑1 (MQ=255)                                                                                                           gcCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCTTCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACCCCCCGcc                               >  4:112446/1‑148 (MQ=255)                                                                                                           gcCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGccc                              <  8:176601/149‑1 (MQ=255)                                                                                                            cccGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGAGCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCt                             <  3:4628/149‑1 (MQ=255)                                                                                                               ggATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGAACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGaa                          <  8:112126/149‑1 (MQ=255)                                                                                                               ggATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGAACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGaa                          <  4:197982/149‑1 (MQ=255)                                                                                                                     ccccTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAgcccga                    >  2:1535/1‑149 (MQ=255)                                                                                                                        cTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACAccc                                   <  8:156591/131‑1 (MQ=255)                                                                                                                        cTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACAccc                                   >  7:156591/1‑131 (MQ=255)                                                                                                                        cTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAgaccgacc                  >  3:262508/1‑148 (MQ=255)                                                                                                                            aTGA‑CGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCa             >  1:292514/1‑149 (MQ=255)                                                                                                                                   ccTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGc       <  2:136247/149‑1 (MQ=255)                                                                                                                                    cttcttTCATCCGGATTATGACCGTCGGCTCCGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCttt                                                                <  8:227984/91‑1 (MQ=255)                                                                                                                                    cttcttTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGaa                          <  3:226984/129‑1 (MQ=255)                                                                                                                                    cttcttTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGaa                          >  4:226984/1‑129 (MQ=255)                                                                                                                                    cttcttTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGca      >  2:179204/1‑149 (MQ=255)                                                                                                                                    cttcttTCATCCGGAATATGACCGTCGGCTCAGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCttt                                                                >  7:227984/1‑91 (MQ=255)                                                                                                                                        tttCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGt  >  5:286137/1‑149 (MQ=255)                                                                                                                                          tCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACGCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGt  <  7:70958/147‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                            GCGGCGGAAACATTGGGGAATTGCGCGTGGTGCTGGGTGAACATGAAACGCTCCTCGTGAATGATCGATGAACGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGA‑CGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTTTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGT  >  NC_002947/616258‑616540

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