Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,684,887	T→A	5.3%	Q79L (CAG→CTG)	PP_1480 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,684,887	0	T	A	5.3%	239.8 / 2.4	75	Q79L (CAG→CTG)	PP_1480	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (36/35);  new base A (2/2);  total (38/37)
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.87e-01

GGAGCTCACTGATTGATTCACCGTGAAAGCCGATCACGTCACGGATGCCGGCAACATGACCGATCAGTAACTGGTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTCCCATTAA‑GCTCGAAGCGCACCCGCGATCCACGGCCTTCGATCGCCTTGGCGCCTATCGATAT  >  NC_002947/1684742‑1685032                                                                                                                                                  |                                                                                                                                                   ggAGCTCACTGATTGATTCACCGTGAAAGCCGATCACGTCACGGATGCCGGCAACATGACCGATCAGTAACTGGTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGtt                                                                                                                                                 >  8:179681/1‑149 (MQ=255)  gAGCTCACTGATTGATTCACCGTGAAAGCCGATCACGTCACGGATGCCGGCAACATGACCGATCAGTAACTGGTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTc                                                                                                                                                <  2:44524/149‑1 (MQ=255)   aGCTCACTGATTGATTCACCGTGAAAGCCGATCACGTCACGGATGCCGGCAACATGACCGATCAGTAACTGGTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTcc                                                                                                                                               <  6:164873/149‑1 (MQ=255)      tCACTGATTGATTCACCGTGAAAGCCGATCACGTCACGGATGCCGGCAACATGACCGATCAGTAACTGGTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAg                                                                                                                                             >  2:48801/1‑148 (MQ=255)       cACTGATTGATTCACCGTGAAAGCCGATCACGTCACGGATGCCGGCAACATGACCGATCAGTAACTGGTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCAGTTCCTgg                                                                                                                                            >  5:224241/1‑148 (MQ=255)         cTGATTGATTCACCGTGAAAGCCGATCACGTCACGGATGCCGGCAACATGACCGATCAGTAACTGGTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAAc                                                                                                                                         <  7:179681/149‑1 (MQ=255)               gattCACCGTGAAAGCCGATCACGTCACGGATGCCGGCAACATGACCGATCAGTAACTGGTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGgcg                                                                                                                                    <  5:221517/148‑1 (MQ=255)               gattCACCGTGAAAGCCGATCACGTCACGGATGCCGGCAACATGACCGATCAGTAACTGGTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCAGTTCCTGGAACAGgca                                                                                                                                    <  6:224241/148‑2 (MQ=255)                attCACCGTGAAAGCCGATCACGTCACGGATGCCGGCAACATGACCGATCAGTAACTGGTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGgcgcg                                                                                                                                  >  4:85878/1‑149 (MQ=255)                    aCCGCGAAAGCCGATCACGTCACGGATGCCGGCAACATGACCGATCAGTAACTGGTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCt                                                                                                                              >  5:64630/1‑149 (MQ=255)                      cgcgAAAGCCGATCACGTCACGGATGCCGGCAACATGACCGATCAGTAACTGGTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTcc                                                                                                                            <  6:64630/146‑1 (MQ=255)                        tGAAAGCCGATCACGTCACGGATGCCGGCAACATGACCGATCAGTAACTGGTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGt                                                                                                                          <  2:78905/149‑1 (MQ=255)                           aaGCCGATCACGTCACGGATGCCGGCAACATGACCGATCAGTAACTGGTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTAcc                                                                                                                       <  2:268564/149‑1 (MQ=255)                            aGCCGATCACGTCACGGATGCCGGCAACATGACCGATCAGTAACTGGTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTa                                                                                                                         >  4:366020/1‑146 (MQ=255)                            aGCCGATCACGTCACGGATGCCGGCAACATGACCGATCAGTAACTGGTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTa                                                                                                                         <  3:366020/146‑1 (MQ=255)                                gATCACGTCACGGATGCCGGCAACATGACCGATCAGTAACTGGTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGa                                                                                                                    >  7:204157/1‑147 (MQ=255)                                      gtcacgGATGCCGGCAACATGACCGATCAGTAACTGGTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAAc                                                                                                                                         >  6:84479/1‑120 (MQ=255)                                       tcacgGATGCCGGCAACATGACCGATCAGTAACTGGTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTg                                                                                                           <  3:85878/149‑1 (MQ=255)                                       tcacgGATGCCGGCAACATGACCGATCAGTAACTGGTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTg                                                                                                           >  8:118793/1‑149 (MQ=255)                                             aTGCCGGCAACATGACCGATCAGTAACTGGTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATcc                                                                                                     <  3:103554/149‑1 (MQ=255)                                                ccGGCAACATGACCGATCAGTAACTGGTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTg                                                                                                  <  8:204157/149‑1 (MQ=255)                                                ccGGCAACATGACCGATCAGTAACTGGTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCAGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTACAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGCCATCCTTg                                                                                                  >  7:355658/1‑149 (MQ=255)                                                 cGGCAACATGACCGATCAGTAACTGGTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCAGTTCCTGGAACAGgcag                                                                                                                                   >  4:186702/1‑113 (MQ=255)                                                 cGGCAACATGACCGATCAGTAACTGGTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCAGTTCCTGGAACAGgcag                                                                                                                                   <  3:186702/115‑3 (MQ=255)                                                          gcccGATCAGTAACTGGTCTTCATCGCTTTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTc                                                                                                <  7:138413/139‑1 (MQ=255)                                                          gACCGATCAGTAACTGGTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCg                                                                                        >  6:34341/1‑149 (MQ=255)                                                           aCCTATCAGTATCTGGTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATTCg                                                                                        <  6:411861/148‑1 (MQ=255)                                                           aCCGATCAGTAACTGGTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCgg                                                                                       >  8:225064/1‑149 (MQ=255)                                                             cGATCAGTAACTGGTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTAcc                                                                                                                       <  3:300447/115‑1 (MQ=255)                                                             cGATCAGTAACTGGTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTAcc                                                                                                                       >  4:300447/1‑115 (MQ=255)                                                                  attAACTTGTCTTCATCTCTGTAATCAATGCGTGTGGCATAGCCGCTGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCATGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATg                                                                                 <  8:355658/146‑1 (MQ=255)                                                                  aGTAACTGGTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATg                                                                                 <  7:225064/148‑1 (MQ=255)                                                                   gTAACTGGTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATc                                                                                                      <  4:218378/126‑1 (MQ=255)                                                                   gTAACTGGTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATc                                                                                                      >  3:218378/1‑126 (MQ=255)                                                                   gTAACTGGTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATgg                                                                                >  6:205003/1‑148 (MQ=255)                                                                   gTAACTGGTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGt                                                                               <  5:210629/149‑1 (MQ=255)                                                                    tAACTGGTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGa                                                                                   <  8:112216/144‑1 (MQ=255)                                                                    tAACTGGTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGa                                                                                   >  7:112216/1‑144 (MQ=255)                                                                          gTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATa                                                                                                                  >  2:46081/1‑107 (MQ=255)                                                                                aTCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTg                                                                          <  8:178746/141‑1 (MQ=255)                                                                                aTCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTg                                                                          >  7:178746/1‑141 (MQ=255)                                                                                       tAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGa                                                                                   >  5:207792/1‑125 (MQ=255)                                                                                       tAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGa                                                                                   <  6:207792/125‑1 (MQ=255)                                                                                         aTCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTCCCATTaa                                                          <  5:205003/148‑1 (MQ=255)                                                                                         aTCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTCCCATTaa                                                          <  7:118793/148‑1 (MQ=255)                                                                                         aTCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTCCCATTAA‑g                                                        >  6:368068/1‑149 (MQ=255)                                                                                          tCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTCCCATTAA‑Gc                                                       >  7:302412/1‑149 (MQ=255)                                                                                          tCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTCCCATTAA‑Gc                                                       >  6:90112/1‑149 (MQ=255)                                                                                            aaTGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTCCCATTAA‑GCTc                                                     <  8:302412/149‑1 (MQ=255)                                                                                                      gCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTCCCa                                                              >  4:388155/1‑131 (MQ=255)                                                                                                      gCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTCCCa                                                              <  3:388155/131‑1 (MQ=255)                                                                                                       cATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAATTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTCCCATTAA‑GCTCGAAGCGCAccc                                          >  3:365317/1‑149 (MQ=255)                                                                                                       cATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTCCCATTAA‑GCTCGAAGCGCAccc                                          <  5:90112/149‑1 (MQ=255)                                                                                                          aGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTCCCATTAA‑GCTCGAAGCGCAc                                            >  3:177892/1‑144 (MQ=255)                                                                                                          aGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTCCCATTAA‑GCTCGAAGCGCAc                                            <  4:177892/144‑1 (MQ=255)                                                                                                                  tAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTCCCATTAA‑GCTCGAAGCGCACCCGCTATCCACgg                               <  5:368068/149‑1 (MQ=255)                                                                                                                    ggTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTCCCATTAA‑GCTCGAAGCGCACCCGCGATCCACGGcc                             <  5:34341/149‑1 (MQ=255)                                                                                                                     gTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTCCCATTAA‑GCTCGAAGCGCACCCGCGATCCACGGCCt                            >  6:195543/1‑149 (MQ=255)                                                                                                                       cATTACATTCATGGTGTTACTCCTGCCTGTTACAGGAACTGGCTCGCGCTCCTTACATGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGAAATGGTCCCATTAA‑GCTCGAAGCTCACCCGCGATCCACCGCCt                            >  4:225222/1‑147 (MQ=255)                                                                                                                        aTTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTc                                                                 >  6:88286/1‑110 (MQ=255)                                                                                                                        aTTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTc                                                                 <  5:88286/110‑1 (MQ=255)                                                                                                                        aTTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTCCCATTAA‑GCTCGAAGCGCACCCGCGATCCACGGCCTTc                          >  8:385231/1‑148 (MQ=255)                                                                                                                        aTTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTCCCATTAA‑GCTCGAAGCGCACCCGCGATCCACGGCCTTc                          >  1:145409/1‑148 (MQ=255)                                                                                                                           aCATTCATGGTGTTACTACTGCCTGTTCCAGGAACTGGGGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGTTGCGGGGGTTGGTACGTGGGCATGGTCCCATTAATGCTCGAAGCGCACCCGCGATCCACGGCCTTCGa                        <  1:48801/148‑1 (MQ=255)                                                                                                                            cATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTCCCATTAA‑GCTCGAAGCGCACCCGCGATCCACGGCCTTCGATc                      <  7:385231/148‑1 (MQ=255)                                                                                                                            cATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATACTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTCCCATTAA‑GCTCGAAGCGCACCCGCGATCCACGGCCTTCGATc                      <  3:265736/148‑1 (MQ=255)                                                                                                                             aTTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTCCCATTAA‑GCTCGAAGCGCACCCGCGATCCACGGCCTTCGATCGc                    >  5:379443/1‑149 (MQ=255)                                                                                                                             aTTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTCCCATTAA‑GCTCGAAGCGCACCCGCGATCCACGGCCTTCGATCGc                    >  4:263649/1‑149 (MQ=255)                                                                                                                             aTTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTCCCATTAA‑GCTCGAAGCGCACCCGCGATCCACGGCCTTCGATCGc                    >  4:197815/1‑149 (MQ=255)                                                                                                                               tCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTCCCATTAA‑GCTCGAAGCGCACCCGCGATCCACGGCCTTCGATCGCCt                  >  3:264063/1‑149 (MQ=255)                                                                                                                               tCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTCGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTCCCATTAA‑GCTCGAAGCGCACCCGCGATCCACGGCCTTCGATCGCCt                  >  6:171743/1‑149 (MQ=255)                                                                                                                                 aTGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTCCCATTAA‑GCTCGAAGCGCACCCGCGATCCACGGCCTTCGATCGCCTTg                >  5:32315/1‑149 (MQ=255)                                                                                                                                       ttACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTCCCATTAA‑GCTCGAAGCGCACCCGCGATCCACGGCCTTCGATCGCCTTGGCGcc           >  3:385710/1‑148 (MQ=255)                                                                                                                                          cTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTCCCATTAA‑GCTCGAAGCGCACCCGCGATCCACGGCCTTCGATCGCCTTGGCGCCTAt        <  5:195543/148‑1 (MQ=255)                                                                                                                                          cTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTCGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTCCCATTAA‑GCTCGAAGCGCACCCGCGATCCACGGCCTTCGATCGCCTTGGCGCCTAt        <  5:171743/148‑1 (MQ=255)                                                                                                                                                cctgTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTACCATTAA‑GCTCGAAGCGCACCCGCGATCCAAGGCCTTCGATCGCCTTGGCGCCTATCGatat  >  7:179534/1‑148 (MQ=255)                                                                                                                                                  |                                                                                                                                                   GGAGCTCACTGATTGATTCACCGTGAAAGCCGATCACGTCACGGATGCCGGCAACATGACCGATCAGTAACTGGTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTCCCATTAA‑GCTCGAAGCGCACCCGCGATCCACGGCCTTCGATCGCCTTGGCGCCTATCGATAT  >  NC_002947/1684742‑1685032

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