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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,684,893	0	A	T	5.1%	254.1 / 2.7	79	L77Q (CTG→CAG)	PP_1480	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (39/36);  new base T (2/2);  total (41/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

TCACTGATTGATTCACCGTGAAAGCCGATCACGTCACGGATGCCGGCAACATGACCGATCAGTAACTGGTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTCCCATTAA‑GCTCGAAGCGCACCCGCGATCCACGGCCTTCGATCGCCTTGGCGCCTATCGATATCAGAAAT  >  NC_002947/1684747‑1685039                                                                                                                                                   |                                                                                                                                                    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                                                                                                                  >  8:118793/1‑149 (MQ=255)                                  tcacgGATGCCGGCAACATGACCGATCAGTAACTGGTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTg                                                                                                                  <  3:85878/149‑1 (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                                                                                                                              >  4:300447/1‑115 (MQ=255)                                                        cGATCAGTAACTGGTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTAcc                                                                                                                              <  3:300447/115‑1 (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                                                                                          >  7:112216/1‑144 (MQ=255)                                                               tAACTGGTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGa                                                                                          <  8:112216/144‑1 (MQ=255)                                                                     gTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATa                                                                                                                         >  2:46081/1‑107 (MQ=255)                                                                           aTCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTg                                                                                 >  7:178746/1‑141 (MQ=255)                                                                           aTCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTg                                                                                 <  8:178746/141‑1 (MQ=255)                                                                                  tAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGa                                                                                          >  5:207792/1‑125 (MQ=255)                                                                                  tAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGa                                                                                          <  6:207792/125‑1 (MQ=255)                                                                                    aTCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTCCCATTaa                                                                 <  7:118793/148‑1 (MQ=255)                                                                                    aTCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTCCCATTaa                                                                 <  5:205003/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                                                   <  4:177892/144‑1 (MQ=255)                                                                                                     aGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTCCCATTAA‑GCTCGAAGCGCAc                                                   >  3:177892/1‑144 (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                                                                        <  5:88286/110‑1 (MQ=255)                                                                                                                   aTTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTc                                                                        >  6:88286/1‑110 (MQ=255)                                                                                                                   aTTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTCCCATTAA‑GCTCGAAGCGCACCCGCGATCCACGGCCTTc                                 >  1:145409/1‑148 (MQ=255)                                                                                                                   aTTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTCCCATTAA‑GCTCGAAGCGCACCCGCGATCCACGGCCTTc                                 >  8:385231/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:197815/1‑149 (MQ=255)                                                                                                                        aTTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTCCCATTAA‑GCTCGAAGCGCACCCGCGATCCACGGCCTTCGATCGc                           >  4:263649/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)                                                                                                                                               ttCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTGGGGATGCGGTCGATGGTACGTGGCAATGGTCCCATTAA‑GCTCGAAGCGCACCCGCGATCCACGGCCTTCGATCGCCtt                        <  5:294349/129‑1 (MQ=255)                                                                                                                                               ttCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTCCCATTAA‑GCTCGAAGCGCACCCGCGATCCACGGCCTTCGATCGCCTTGGCGCCTATCGATATCAGaa    <  3:263649/149‑1 (MQ=255)                                                                                                                                               ttCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTCCCATTAA‑GCTCGAAGCGCACCCGCGATCCACGGCATACGATCCCCtt                        >  6:294349/1‑129 (MQ=255)                                                                                                                                                  cAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTCCCATTAA‑GCTCGAAGCGCACCCGCGATCCACGGCCTTCGATc                             <  7:318840/121‑1 (MQ=255)                                                                                                                                                  cAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTCCCATTAA‑GCTCGAAGCGCACCCGCGATCCACGGCCTTCGATc                             >  8:318840/1‑121 (MQ=255)                                                                                                                                                  cAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTCCCATTAA‑GCTCGAAGCGCACCCGCGATCCACGGCCTTCGATCGCCTTGGCGCCTATCGATATCAGAAAt  >  6:17658/1‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                                    TCACTGATTGATTCACCGTGAAAGCCGATCACGTCACGGATGCCGGCAACATGACCGATCAGTAACTGGTCTTCATCGCTGTAATCAATGCGTGCGGCATAGCCGCGGTAGGTCATTACATTCATGGTGTTACTCCTGCCTGTTCCAGGAACTGGCGCGCGCTCCGTACCTGATAGGATTTGGCATCCTTGTCGGGATGCGGTCGATGGTACGTGGCAATGGTCCCATTAA‑GCTCGAAGCGCACCCGCGATCCACGGCCTTCGATCGCCTTGGCGCCTATCGATATCAGAAAT  >  NC_002947/1684747‑1685039

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