Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	537,458	T→A	9.0%	noncoding (54/66 nt)	PP_mr07 →	P26 ncRNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	537,458	0	T	A	9.0%	189.6 / 7.8	67	noncoding (54/66 nt)	PP_mr07	P26 ncRNA
Reads supporting (aligned to +/- strand):  ref base T (34/27);  new base A (4/2);  total (38/29)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.09e-01

CTCAGGTTGTAGCTGAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGAC  >  NC_002947/537311‑537605                                                                                                                                                    |                                                                                                                                                    ctcAGGTTGTAGCTAAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCt                                                                                                                                                     <  8:99562/148‑1 (MQ=255)    aGGTTGTAGCTGAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATg                                                                                                                                                 >  3:148973/1‑149 (MQ=255)      gTTGTAGCTGAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGcc                                                                                                                                               >  8:78202/1‑149 (MQ=255)         gTAGCTGAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCAcc                                                                                                                                            >  2:90966/1‑149 (MQ=255)            gCTGAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGc                                                                                                                                         >  3:255530/1‑149 (MQ=255)                aaGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCttt                                                                                                                                     >  4:303633/1‑149 (MQ=255)                aaGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCCGGTGGCATAAGCCACCAGccat                                                                                                                                      >  2:178694/1‑146 (MQ=255)                 aGGCGTTTCGAAAGAAGCCGCTGAAGACGCAAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCCGGTGGCATAAGCCACCAGccatca                                                                                                                                    <  1:178694/149‑5 (MQ=255)                   gCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCCAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGGTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTc                                                                                                                                  <  5:314581/149‑1 (MQ=255)                    cGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTcc                                                                                                                                 <  7:78202/149‑1 (MQ=255)                     gTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCATAAGCCACCAGc                                                                                                                                         >  1:94242/1‑140 (MQ=255)                     gTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCATAAGCCACCAGc                                                                                                                                         <  2:94242/140‑1 (MQ=255)                           aaaGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTAtt                                                                                                                           <  2:277236/148‑1 (MQ=255)                           aaaGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATtg                                                                                                                          <  7:118143/149‑1 (MQ=255)                            aagaagCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATtgg                                                                                                                         <  4:148973/149‑1 (MQ=255)                            aagaagCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATtgg                                                                                                                         >  7:259654/1‑149 (MQ=255)                                agCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTGGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTCCCGTTATTGGTGGc                                                                                                                     <  6:277470/149‑1 (MQ=255)                                   cGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCgg                                                                                                                  >  8:220715/1‑149 (MQ=255)                                        aaGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAgg                                                                                                             >  2:82245/1‑149 (MQ=255)                                        aaGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAgg                                                                                                             >  6:88554/1‑149 (MQ=255)                                        aaGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATCCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAgg                                                                                                             >  7:12784/1‑149 (MQ=255)                                         aGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGc                                                                                                                     >  4:308772/1‑140 (MQ=255)                                         aGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGc                                                                                                                     <  3:308772/140‑1 (MQ=255)                                         aGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGc                                                                                                            <  6:50173/149‑1 (MQ=255)                                             gCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGt                                                                                                        <  3:91438/149‑1 (MQ=255)                                               tAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCAc                                                                                                                                             <  3:302482/110‑1 (MQ=255)                                               tAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCAc                                                                                                                                             >  4:302482/1‑110 (MQ=255)                                                 agaagaagCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTccc                                                                                                     >  8:202056/1‑148 (MQ=255)                                                  gaagaagCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGa                                                                                                   <  1:216915/149‑1 (MQ=255)                                                   aagaagCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGAt                                                                                                  <  5:88554/149‑1 (MQ=255)                                                     gaagCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATaa                                                                                                >  2:304148/1‑149 (MQ=255)                                                       agCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCTAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAAc                                                                                               <  6:57052/148‑1 (MQ=255)                                                                aagCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGc                                                                                        >  8:86234/1‑146 (MQ=255)                                                                   cAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGaca                                                                                  >  1:321689/1‑149 (MQ=255)                                                                          tAACGTCGTGCTGACTTATTCTCGACGTTGCGTTGCCTGCCTGAGCGTTGCGCGAACGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTTTTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCaa                                                                                      <  7:211650/138‑1 (MQ=255)                                                                          tAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCAtgt                                                                           >  3:127503/1‑149 (MQ=255)                                                                          tAAAGTCGAGCTGAAGTAATCTCGACCATGCGAAGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTAAGCCACCGGCCATTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCCAAAACGACCGGCAAGacaacggcag                                                                           >  8:211650/1‑142 (MQ=255)                                                                            aaGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTAt                                                                                                                            <  1:60898/98‑1 (MQ=255)                                                                            aaGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCAtgt                                                                           <  1:90966/147‑1 (MQ=255)                                                                            aaGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATgtgg                                                                         <  1:304148/149‑1 (MQ=255)                                                                            aaGTCGAGCTGAAGTAATCTCGACCTTGCGATGCAAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTAt                                                                                                                            >  2:60898/1‑98 (MQ=255)                                                                               tCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATgtgggtg                                                                      >  1:91888/1‑149 (MQ=255)                                                                                     tGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCa                                                                 >  3:97326/1‑148 (MQ=255)                                                                                       aaGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAAc                                                              <  7:220715/149‑1 (MQ=255)                                                                                            aTCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAg                                                          <  4:127503/148‑1 (MQ=255)                                                                                             tctcGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAggg                                                        >  7:11968/1‑149 (MQ=255)                                                                                                 gACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGttt                                                    <  4:255530/149‑1 (MQ=255)                                                                                                   ccTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGc                                                  <  3:303633/149‑1 (MQ=255)                                                                                                       gCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGa                                              >  6:13830/1‑149 (MQ=255)                                                                                                         gATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGAtt                                            >  6:326136/1‑149 (MQ=255)                                                                                                         gATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGAtt                                            >  6:313157/1‑149 (MQ=255)                                                                                                             ccAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCt                                        >  7:66377/1‑149 (MQ=255)                                                                                                             ccAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCAACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGTTCCCGATAACGAGCTGCaa                                                                                      >  3:152024/1‑103 (MQ=255)                                                                                                              cAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATCCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTg                                       <  8:12784/149‑1 (MQ=255)                                                                                                                 ccTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGgct                                    >  1:244032/1‑149 (MQ=255)                                                                                                                   tGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGgctgc                                  >  2:177957/1‑149 (MQ=255)                                                                                                                   tGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGgctg                                   >  5:105909/1‑148 (MQ=255)                                                                                                                    gAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGgct                                    <  6:133916/146‑1 (MQ=255)                                                                                                                    gAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGgct                                    >  5:133916/1‑146 (MQ=255)                                                                                                                         ttGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTccctc                            >  5:296372/1‑147 (MQ=255)                                                                                                                            aGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCgg                         <  5:313157/149‑1 (MQ=255)                                                                                                                            aGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCgg                         <  8:259654/149‑1 (MQ=255)                                                                                                                              cGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTTCTCCcggcgg                         >  6:20084/1‑147 (MQ=255)                                                                                                                                    gCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGcc                 <  1:82245/149‑1 (MQ=255)                                                                                                                                           gccggtggCATAAGCCACCAGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAg          >  6:155587/4‑149 (MQ=255)                                                                                                                                           gCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAg          >  3:144915/1‑149 (MQ=255)                                                                                                                                            ccggtggCATAAGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAg          >  8:155671/3‑148 (MQ=255)                                                                                                                                            cTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGc         <  7:202056/149‑1 (MQ=255)                                                                                                                                                    ttATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCaaa                                                               <  5:161768/87‑1 (MQ=255)                                                                                                                                                    ttATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCaaa                                                               >  6:161768/1‑87 (MQ=255)                                                                                                                                                    ttATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGAc  >  7:237460/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                    CTCAGGTTGTAGCTGAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGAC  >  NC_002947/537311‑537605

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